Recent mathematical developments in the Skyrme model 



T. Gisiger and M.B. Paranjape 

Groupe de Physique des Particules, Universite de Montreal, CP. 6128, succ. 
centre-ville, Montreal, Quebec, Canada, H3C 3J7 



Abstract 

In this review we present a pedagogical introduction to recent, more mathematical 
developments in the Skyrme model. Our aim is to render these advances accessible 
to mainstream nuclear and particle physicists. We start with the static sector and 
elaborate on geometrical aspects of the definition of the model. Then we review the 
instanton method which yields an analytical approximation to the minimum energy 
configuration in any sector of fixed baryon number, as well as an approximation 
to the surfaces which join together all the low energy critical points. We present 
some explicit results for B = 2. We then describe the work done on the multibaryon 
minima using rational maps, on the topology of the configuration space and the 
possible implications of Morse theory. Next we turn to recent work on the dynamics 
of Skyrmions. We focus exclusively on the low energy interaction, specifically the 
gradient flow method put forward by Manton. We illustrate the method with some 
expository toy models. We end this review with a presentation of our own work 
on the semi-classical quantization of nucleon states and low energy nucleon-nucleon 
scattering. 



1 Introduction 



The Skyrme model[^ was first proposed by T.H.R. SkyrmeQ in the sixties, as 
a revolutionary idea for incorporating baryons in the non-hnear sigma model 
description of the low-energy interactions of pions. This sigma model consists 
of a unitary matrix valued field U{x, t) of dimension 2 x 2 or 3 x 3 depending 
on the number of hght quark flavours that are considered. The dynamics is 
described by the Lagrange density 

C = ii{U^d^UU^d^U) (1) 

^ For an interesting compilation of the life, work and influence in physics of T.H.R. 
Skyrme, see Selected Papers, with commentary, of Tony Hilton Royle Skyrme (World 
Scientific Series in 20th Century Physics- Vol. 3), G.E. Brown editor 
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where is the pion decay constant. Skyrme noted the existence of topologi- 
cally non-trivial field configurations of finite energy. These were however, un- 
stable against collapse, which can be adduced by simple application of scaling 
arguments. Skyrme then added a higher derivative term to the Lagrange den- 
sity rendering these configurations stable. This term is now called the Skyrme 
term and the Skyrme Lagrange density Csk is given by 

= tiiU^d,UU^d>^U) + tr([f/t9^f/, U^dM?) (2) 



where e is a new, dimensionless coupling constant. Since each derivative cor- 
responds to a momentum, this term is clearly of higher order in the low- 
energy (momentum) approximation. Skyrme proposed the interpretation of 
these topological solitons (stable, localized, finite-energy solutions of the clas- 
sical equations of motion) as the nucleons and identified the topological wind- 
ing number of the soliton with the baryon number. The technology of quan- 
tum field theory in the sixties was not sufficiently advanced to treat solitons 
and it took almost twenty years before the ideas of Skyrme were revived by 
Balachandran et al and Witten[^], and vindicated with surprising accord 
with experiment [|] . Witten[^ described another topological density which 
should be added in the effective action, the celebrated Wess-Zumino-Novikov- 
Witten|J|,| term 



where U G SU(3), N & TZi {D^ is in fact a 5-dimensional manifold with only 
its boundary giving the usual 4-dimensional space-time). Witten^ showed its 
relation to the underlying microscopic theory of the strong interactions QCD, 
with the number of colours giving the quantized coefficient in Ty/zNW- Since 
these seminal papers there has been an enormous amount of work relating the 
Skyrme model to phenomena in nuclear and particle physics, for instance 
targeting the spectrum of excitations of baryons0, the inclusion of strange 
degrees of freedom in the model P], the nucleon-nucleon potential P,p!0|, scat- 
tering 7C — N states[|ll[], high density baryon matter as a Skyrme crystal p!2 
and the nucleon-anti-nucleon annihilation [pr^, p!^JT5|Jl6|Jl7| , p!8[| to name a few. 
We will not consider these developments in detail here and refer the in- 
terested reader to the literature and to the many excellent reviews on the 
subject 



Concurrently, there were certain mathematical advances in the Skyrme model 
which were not strongly based on making any contact with phenomenology, 
but moreover to understand the mathematical content of the model. These 
concerned two main areas, that of the exact nature of the minima or critical 



2 



points of the static energy functional and secondly a description of the scat- 
tering of the corresponding solitons. It turns out in a certain approximation, 
that of low energy scattering, these two aspects are not unrelated. There were 
also some auxiliary mathematical and physical results, concerning geometri- 
cal insights into the model. Our review will primarily focus on these relatively 
recent mathematical developments. We start with the static sector and intro- 
duce the Skyrme model on a general Riemannian manifold. Then we present 
the methods using instantons and rational maps to obtain useful Skyrme con- 
figurations, followed by a short description of Morse theory and its application 
in the model. Next we move to the dynamics where we treat the gradient flow 
method put forward for studying soliton scattering. We terminate with it ap- 
plications to the baryon number 2 sector of the model. Most of the advances 
which we will consider were made by N.S. Manton, among others (our original 
contributions to this subject are secondary). Since these advances generally 
use the language and formalism of differential geometry, the reader should be 
familiar with these notions (any standard course on tensor analysis/differential 
geometry /general relativity should be adequate 



ga,24,25,2f 



We start the study of the statics of the model in section 3 with the geometrical 
aspects of the model as first discussed by Manton and Ruback [p7| and then by 
Mantonp8[ and Lossp9[. We show in detail how the model can be understood 
as a theory of elasticity in curved space. This fruitful approach allows us to 
explain, for example, why the Skyrmion (the lowest energy solution of the 
model with baryon number 1) does not saturate the Bogomolnyi bound in 
ordinary space, but does so for a space with great curvature. This situation is 
connected to chiral symmetry restoration and deconfinement. In section 3.1 we 
first review the theory of non-linear elasticity, both in flat and curved space, 
as well as some tools of differential geometry. We then move on to show how 
those tools can allow us to connect elasticity and the Skyrme model. This 
shows how it is possible to formally interpret a classical field theory as an 
elasticity theory. 



Section 3.2 describes the instanton method put forward by Atiyah and Manton |[30 
to approximate critical points of the Skyrme energy functional. Indeed, the 
Skyrme model has always been plagued (as T.H.R. Skyrme realized already 
in the sixties) by the absence of analytical solutions: to obtain solutions of 
the model, one always has to solve partial differential equations numerically, 
except in the simplest cases like the B = 1 Skyrmion, where one has only 
to solve an ordinary differential equation numerically. The instanton method 
permits one to obtain analytical expressions (in the simplest cases) or at least 
limit the difficulties to solving ordinary differential equations. The tradeoff is 
that in this method (which is by nature approximate) the error introduced is 
hard to estimate, except by comparing the resulting approximation with the 
full "exact" numerical solution of the problem. The method seems neverthe- 
less to work reasonably well, and has been helpful in investigating the bound 
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states in the sector B = 2, 3 and 4, and Skyrme crystals as well 
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,32 



,33 



The next subsection 3.3 presents the use of rational maps|3^ in finding ap- 
proximations to minima with particular symmetries for sectors comprising of 
a range of baryon numbers. Rational maps are mappings of the complex plane 
to itself comprised of ratios of polynomials. The choice of the coefficients of 
the polynomials can encode the mappings with complicated symmetries, es- 
pecially when viewed via stereographic projection as mappings of 5*^ —>■ S"^. 
The former S"^ corresponds to the angular degrees of freedom about a Skyrme 
configuration while the latter corresponds to a selected S"^ in the group SU(2) 
obtained via the Hopf projection as will be explained in subsection 3.3. 



We give in subsection 3.4 a short introduction to Morse theory and how it could 
be used to find new solutions to the equations of motion Morse theory 
relates the existence of critical points of a function defined on a manifold to 
non-trivial topological aspects of the manifold. Here the function in question 
is the energy functional defined on the manifold of the space of all static 
field configurations. Although the method has not borne fruit in the Skyrme 
model, it has already been useful in the analysis leading to the sphaleron of 



the standard electroweak model[36 



In the next section (section 4), we present some recent developments in the 
study of the dynamics of the model, namely Skyrmion-Skyrmion and nucleon- 
nucleon scattering. 



We start by presenting in section 4.1 Manton's method[0] for truncating 
the degrees of freedom of a system, thereby possibly rendering it tractable. 
This formalism was published by Manton after his work on BPS monopole 
(topological solitons of the massless SU(2) Higgs model) scattering |3^. It is 
a much more general method in the sense that it is applicable not only to 
soliton problems (namely scattering thereof) but to many systems with a great 
number of degrees of freedom under the right circumstances. Even though the 
underlying ideas are intuitively simple, the formalism needs to be expressed 
using differential geometry. This is why in section 4.2 we introduce this method 
using a number of simple examples to help us emphasize the ideas behind the 
method and to get a feel for how it works. Results obtained using this method 
are then compared with "exact" numerical results. In section 4.3 we present 
the application of the method to the case S = 2 in the Skyrme model. 



We close this review with a discussion of our work on Skyrmion-Skyrmion ||39|| 
and nucleon-nucleon|^0|] scattering using Manton's method but in a semi- 
classical, rather than purely quantum mechanical perspective. 
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2 The Skyrme model: Lagrangian and notation 



The Skyrme model including mass is described by the Lagrange density, 



-]-mlii{U + U^) (4) 



where U{x) is a unitary matrix valued field. The Lagrangian is also often 
written in terms of F^^ = 2/,^ but this is only a matter of convention. The 
quantum fiuctuations of U{x) represent the low energy mesons made up of 
quark-anti-quark pairs. In the chiral limit, all light quarks {u, d, s) are massless 
and degenerate and the corresponding flavour symmetry dictates that U{x) is 
a 3 X 3 matrix. In this review, however, we will be more focused on the fact 
that U{x) is a unitary matrix; the simplest example of this is a 2 x 2 matrix, 
the only case that we will consider here. Phenomenologically, this means we 
consider the explicit breaking of the SU(3)f symmetry to be large. The Skyrme 
Lagrangian (^) corresponds to the first two terms of a systematic expansion 
in derivatives of the effective Lagrangian describing low energy interactions of 
pions plus the mass term. It should be derivable from QCD hence /jr, e and 
mjr, the pion mass, are in principle calculable parameters. These calculations 
are actually unfeasable and we take /,r, e and m^r from phenomenological fits. 
We then find /,r to be in the range of 130-190 MeV and e ~ 5. In this article, 
we will always take m,r for simplicity. 

The energy functional, coming from the static part of the Lagrangian, can 
be written in a more elegant way by adopting a convenient choice of units 
(so-called natural units) p7|] : 



^sk = / d'x 



-- tT{U^diUU^diU) - — tT{[WdiU,U^djU]^] 
2 16 



(5) 



where the unit of energy is now f-j^/^e ~ 6 MeV and the unit of length 2/e/^ ^ 
0.6 fm. 

The baryons arise as topological solitonic solutions of the equations of motion . 
These topological solitons correspond to non-trivial mappings of IR^ plus the 
point at infinity into SU(2): 

U{x) : R3 + cx) ^ SU(2) = Si (6) 
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Fig. 1. Artistic impression (by T. Gisiger) of the field of the Skyrmion or hedgehog. 
But topologically 



+ oo = S^; (7) 

thus the homotopy classes of mappings 

U{x) : 5^ (8) 



which define 113(5'^) = 7Zi, characterize the space of configurations. The topo- 
logical charge of each sector is given by 

N = f ci^f gijk tT{U^diUU^djUU^ch,U), (9) 



which is an integer and is identified with the baryon number 

The solution with baryon number = 1 with lowest energy is the Skyrmion. 
It is parametrized presumably (since this has still not been rigorously, math- 
ematically proven) by the following expression 

U{x) = e'^"-"^(") (10) 



where r are the Pauli matrices, f is the unit position vector and r is its length. 
The field U{x) therefore points radially, and is sometimes called "hedgehog" 
(see figure |I]). /(r) is a function to be obtained from the equations of motion 
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of the theory by replacing the ansatz ([T0|) in the static equations of the theory: 



diUU\ 



dyUU\d;UU^ 



which give for (|T0|) 



+ 2 sin^ F)F" + -pF' + sin 2fF'^ - - sin 2F 



sin^ F sin 2F 
? 



(12) 



where p = 2efT^r and /(r) = F{p). This equation has to be solved numerically 
using the boundary conditions imposed on the ansatz by the requirements 
that it have baryon number unity and finite energy, namely /(O) = vr and 
/(+oo) = 0. The asympotic behaviour of /(r) for large values of r is obtainable 
from CO): 



fir) 



K 

^2' 



(13) 



Using numerical integration we find that k ~ H^Ell (The equation 

for /(r) can equally well be obtained by replacing the ansatz (|TUp in the energy 
functional (|^), extremizing relative to f{r), and integrating over 6 and (p). 



Using the definitions 

Cl{U) = tT[r^U^d,U] 



(14) 
(15) 
(16) 



the Skyrme Lagrange density can be written in the following fashion 



n{u)^-n{uYn{u),-n{uY 
-n{u)^-n{uyn{u),-n{uY 



(17) 



where the "■" represents, for example, TZ^-TZ'^ = TZ^Tl^'^. We can separate this 
into a kinetic energy T which is the part quadratic in time derivatives and a 
potential energy V without any time derivatives: 



f 1 

r = ^7^o • 7^o — ^ 

2 2e2 



7^o ■ 7^i 7^o ■ 7^i - 7^o • 7^o n-, ■ 
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f2 1 



7^i ■ 7^i 7^j • 7^j - 7^i • 7^j 7^j • 7^i 



(19) 



The parametrization of a Skyrmion with (iso) orientation defined by a time 

— * 

dependent SU(2) matrix A{t) and position R{t) is 



[/(f , t) = f/(f - R{t)) A{t)\ 



(20) 



This gives a Skyrmion 6 degrees of freedom. After replacing this ansatz in the 
Skyrme Lagrangian and integrating over all space, we findH]: 



L = -2M + -MR 2 + 2A (7^^(A) 7^^(A)) 



(21) 



where 



+ 00 

M= An J r'^dr 





X 



dr J 



+ 



1 sinV r sinV ^ ^/df 

2g2 y,2 y y,2 y Qy, 



(22) 



is the mass of a Skyrmion and (see for instance ||20| 



A = (c/tt)^ j r'^dr sin^ / 



1 + 



4 /^.^ftoV 



(23) 



is its moment of inertia. This roughly gives the Skyrmion a mass of 850 MeV 
(quite close to the nucleon mass), or in natural units of energy ~ 1.23 x 
127r2[^. The moment of inertia is roughly equal to (1/195) MeV^^[^. 



3 Study of the statics of the Skyrme model 



We now start the first part of this review which is devoted to the study of the 
static sector of the model. Most of the material concerns geometric aspects 
of the model, but a few pages are devoted to the instanton method, rational 
maps and Morse theory. 
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3.1 Geometry of the Skyrme model 



In this subsection we will first briefly review the theory of non-linear defor- 
mations of a body. Most of the material comes from Ogden|]42| and readily 
generalizes to the case of a fleld theory. We then take up the case of the 
Skyrme model deflned on various spatial manifolds, presenting the flndings of 
Manton[gB| and Manton and RubackfS^]. 



3.1.1 Non-linear deformation of a body 

Let us consider a body B at rest, i.e. , free of oscillations or interior motion: 
the action of exterior forces (like gravity) if any, and that of the interior forces 
caused by the nature of the body (the interactions between the atoms of the 
body, for instance) exactly cancel each other. We will call this state of the 
body its initial conflguration. Let us now change the shape of the body by 
applying forces to it, until it reaches a new static conflguration. Then it is 
possible to characterize this action by a function x which we now deflne. Let 
X be the initial position of a given point P of the body, and x its position 
in the flnal conflguration. Then we can deflne the deformation function x 
mapping X to x for every point of the body: 



X^x = x{X). (24) 

For our needs, we will restrict ourselves to x being continuous and twice 
differentiable. In the case of a real dynamical system, x could depend on 
time but here we will only consider static conflgurations. See Figure for an 
example of deformation parametrized by a function x- 

To make things more concrete, let us introduce coordinate systems in the 
initial and flnal conflgurations. Using two different coordinate systems to de- 
scribe the deformation of a body can seem like overkill but we follow this 
course for two reasons. The flrst is that it enables one to choose coordinate 
systems which best suit each conflguration, simplifying the computations that 
follow (which are usually quite complex in real life problems). The second, 
more important reason is that it will make the jump to Manton's elasticity 
formalism in the frame of fleld theory easier and more natural. We note that 
we will only be working with bases of vectors which are locally orthogonal, and 
that most of the following equations are written in cartesian coordinates. All 
tensorial equations are readily generalized to arbitrary curvilinear systems. 

Let O be the origin of the system in the initial conflguration and {-Em} its 
base. Following the convention used by Manton[^], we will use the indices 
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Fig. 2. Example of deformation of a solid (left) before the deformation and (right) 
after a simple stretching along the z axis. 

m,n,o,p to denote basis vectors. The position X of a point of the body is 
defined by 



= X^E^ (25) 

where summation over m is imphcit. X then has coordinates X™ in this co- 
ordinate system. 

Let o and {^^} be the origin and base of the coordinate system we choose to 
describe the final configuration, respectively (We will use the indices ii,v p,a 
here). Then we express x in the final basis as 

^x^i,. (26) 



Then equation ( ^4]) can be written in terms of the components 

x^ = x^(X-). (27) 

The deformation of the body is completely contained in the map x- How- 
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ever the map contains much redundancy which we will isolate next. Under 
deformation an element of length dX transforms according to 

df = ^ dX (28) 

dx 



or in coordinates 

dx^" = dX"^. (29) 



Let us introduce the second order tensor 

A = A^^ ® (30) 



with 



A 



fira 



dX^ 



(31) 



which is called the deformation gradient (relative to the reference configura- 
tion). This is just the local Jacobian of the deformation defined by (p^. 

Intuitively it is clear that A completely represents the deformation, but as we 
shall demonstrate, it is also sensitive to the choices of the bases {-Em} and 
{■Cm} which have no physical content. For instance, a physical quantity such 
as the amount of energy stored in the body by the deformation should not 
depend on these choices. 

Let us consider new bases {E'^} and {^'^} defined by 

-E'm = ^mn-En (32) 

S = W,,C (33) 



where Vt and W are orthogonal matrices. The new bases correspond to a new 
choice of orientation for the bases in the initial and final configuration. The 
tensor A expresses itself in each basis as follows: 



A — A^^ -Em 

= Al,C(^K (34) 

where 

A',^ = W^^A^rain^Un (35) 



11 



and ff^ is the transpose of Q. Therefore A depends strongly on the choice 
of the bases. In the normal theory of elasticity, A is considered to be non 
singular: no dX can be deformed into a line element dx with zero length. 
Such an annihilation of a line element would imply an infinite force acting 
on the body, which is unrealistic. But in the next subsections, devoted to the 
study of the Skyrme model as a theory of elasticity, we will see that there is 
physical meaning to a singular or in fact zero tensor A. 

A step toward a better description of the deformation is to use the right 
Cauchy-Green deformation tensor D = AA"^. It is symmetric and positive 
definite by construction and under the change of coordinate axes mentioned 
above, D ^ D' = WDW"^ which is physically more sound but still sensitive 
to the choice of the coordinate system of the final configuration. 

The right CaTichy- Green tensor is related in a simple way to the usual strain 
tensor of elasticity. To show this let us define the displacement vector u: 

u{X) = f - X = x{X) - X (36) 



which measures how much a point moves during the deformation. It is usually 
assumed to be very small in the ordinary, i.e. linear, theory of elasticity. Define 
G the displacement gradient 

^,^ = (^-1^ = ^. (37) 



Then one can define two tensors of strain: 



E = -[G + G^ + G^G) 
= i(^AT-I) = i(D-I) (38) 

and 



F^-{G + G'^ + GG'^) 

^i(A^^-I). (39) 



In component formulation this gives : 



E — — \ 



^ aX^^^^. ^IX^^^IX^) ^^^^ 
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and 



F, 



E is the ordinary strain tensor considered in Ogdenp2| while the reader will 
recognize F to be the strain tensor defined in the theory of elasticity, see for 
example Landau[0. Notice that the physical content of D and E coincide 



since they differ by a translation and a factor, as seen from equation (B^). 

The non-redundant description of the deformation is furnished by functions of 
the matrix D which are invariant under conjugation by orthogonal matrices. 
These functions are given by the secular (determinant) equation of the matrix 
D: 

det(D - 61) = 5^ + h{D)6^ + h{D)5 + h{D) = 0. (42) 



One finds 

h=iiD (43) 

h = \[{tvDf-tr{D')] (44) 

/3 = detD. (45) 

The set of eigenvalues of D is invariant under orthogonal conjugation, as are 
the coefficients Ji, which can be verified easily. The set of invariants Ji is 
complete since the eingenvalues of D are uniquely determined by equation 



These eigenvalues can in fact be expressed in terms of the eigenvalues of a 
matrix obtained from A. This can be seen via the following analysis. Any 
non-singular, real matrix A admits a polar decomposition 



A = RU (46) 
= VR (47) 

where R is an orthogonal matrix and U and V are symmetric, non-singular, 
real matrices . Let -U; and Ai denote the eigenvectors and corresponding 



eigenvalues of U 

Uui = X-A i = 1, 2, 3 (4J 
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where the Ui can be chosen orthonormal since U is a symmetric, real (hermi- 
tian) matrix. Then 



Aui = RUui = Xi{Rui) = V{Rui), (49) 

hence Rui is an eigenvector of V with the same eigenvalue Ai. 
Now 

D = AA^ = V'^ = RU'^R^, (50) 

thus 

D{Rui) = X-{Rui) (51) 



and the three invariants (^31) , (^if) and (|45|) are 

h = \l + \l + Xl (52) 
I2 = A^A2 ~l~ A2A3 -|- A3A^ (53) 
/3 = A?A^A^. (54) 

These invariants will be useful in writing the energy stored in the body due 
to the deformation. The same applies to the Skyrme model, as we will show 
below. We can give a physical interpretation of the three invariants in terms 
of stretching of a set of vectors and various associated geometrical quantities. 

The simplest invariant is which is related to the change in the volume 

defined by three non-coplanar vectors dX'-^-', dX*^^) and dX'-^-', and their 
images dx^^\ dx*^^^ and dx*-^-* under the deformation. By equation (pH]), we 
have 

dx = A-dX (55) 
which gives 

dv = det A dV = A1A2A3 dV (56) 



where 



dV= dX(^) ■ dX(^) X dX(2) (57) 
dt. = df ■ df X df (58) 
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assuming that dx^^^ and dX^'^ form right handed triads. So the third invariant 
is just the square of the change of a volume under the deformation. 



Let us define the stretch of a hne element of the body as measured by the 
ratio between the length of the initial line element and that of its image under 



the deformation. From (pSf ) and (|3T|) we have 



m\ dx\ = AM\ dX\ (59) 

where m, M are unit vectors along the direction of dx and dX and | . . . | 
indicates the length. Taking the norm of both sides of (^) we define 

A(M) = = ^M-A^^A-M (60) 

IdXl 



which gives the stretch in the direction M at X. Taking M to be colinear to 
the n-th eigenvector of A^A with eigenvalue A^, we get 

X\Un) = Xl (61) 



The first invariant Ji is then just the sum of the squares of the stretch along the 
three eigenvectors of U. These eigenvectors correspond exactly to the usual 
principal directions of strain. Indeed from ( PU] ) we see that the matrix F 
corresponds to the usual definition of the strain tensor. However 

F = i(f/2-l) (62) 



hence 



Fu, = ^{X^-l)u, (63) 



thus showing that the Ui correspond to the principal directions of strain. 

The second invariant is the most subtle and is related to the change under 
the deformation of the area elements defined by an orthonormal triad. An 
orthonormal triad -Oi defines the three area elements 

^eijki)j X V], i,j, k = l,2, 3. (64) 



Under the deformation these are transformed to 

ieijk(A{)j) X (Avi,). (65) 
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Computing the norm and summing over i yields 



eijkeiim(^t^j X Av],) ■ {Avi x Avm) 
{Avj X Av]^) ■ {Avj X Avi,) 
{Avj ■ Av^){Av^ ■ Avi,) - {Avj ■ Av^){Avi ■ Av^] 



With 



Vi = Mk 



(66) 
(67) 
(68) 

(69) 



where vf are the elements of an orthogonal matrix (since are orthonormal) , 
we obtain 



AVi — V^AUi: — V^X]f:RU]f: 



(70) 



so 



Ai}^ ■ AV]<^ = wjf ™AlAm(-R'Ul • RUm) 



Hence the squared norm is now equal to 



(71) 
(72) 



>HA?)(«A^) - {v]viX^){vf^v^Xl 



1 m 1 



(73) 

(74) 

(75) 
(76) 



3.1.2 Geometrical framework for the Skyrme model 

A field corresponds to a mapping tt from the manifold S of ordinary space 
to the manifold E of the target space. A field theory further specifies the 
dynamics obeyed by the field via the Euler- Lagrange equations obtained from 
the Lagrangian of the theory. We will limit ourselves to the case where both S 
and S are 3-spheres, the former being a 3-sphere of radius L, while the latter 
is the 3-sphere of the SU (2) group corresponding to isospin. The case of the 
ordinary Skyrmion defined on is obtained by taking L +00. 
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This modest generalization of the Skyrme model allows for a non-trivial ap- 
plication of the geometrical formalism of non-linear elasticity theory which we 
have just elaborated. Varying the radius L allows us to cover the cases from 
extreme to zero curvature. 

The map vr describes the Skyrme field, whether it is a group of waves with zero 
baryon number, a Skyrmion or a heavy nucleus of high baryon number. Even if 
TT bears similarity with the deformation map x of the above elasticity theory, it 
is different in a fundamental way: it maps a curved space onto another curved 
space. X only maps a set of points in IR^ onto another set of points in IR^. 
Apart from this fundamental difference, there are many similarities between 
the treatment of elastic bodies and field theories. 

The initial spatial manifold with a given metric comes equipped with a tangent 
space at each point. The tangent space Tp(S') of S at the point p has a natural 
basis {d/dp^}. Although linearly independent, these vectors are not necessarily 
orthonormal. By a linear transformation, we can construct an orthonormal 
basis {cm} (m = 1, 2, 3) of Tp(S') given by 

■ d 

em = e'^Q- (77) 

where the coefficient of the linear transformation are called the dreibein. Doing 
the same construction at every point p oi S defines an orthonormal frame at 
each point of S. We will follow the same convention as Manton and use indices 
i,j,k,l with reference to the coordinate basis and m,n,o,p with reference 
to the orthonormal basis. If the reader has difficulty with these notions, we 
recommend the references already mentioned [p^, p^ , p5| ,p^. 

With coordinates vr" on S and its given metric r"^ {a, (3, 6, 7 coordinate indices 
and n, u, p, a orthonormal basis indices), consider the image of the orthonormal 
frame field {cm} under 71. We let tt" = 7r"(p') be the coordinates of the image 
of for efficiency of notation. According to the transformation law of a vector 

e--^.^- (78) 



It is evident that the lengths and directions of the orthonormal triad are 
changed: it is generally no longer orthonormal. As long as the Jacobian of the 
transformation 



is non-singular, the image triad defines a basis of the image tangent space. This 
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is the generic case but not at all the relevant one in many physical situations, 
as we will see below. 

The (inverse) metric of space S is given by 

tiJ = S^^e^ei (80) 



where 6"^^ is the usual Kronecker delta. Its image under vr is 
dp^ dp} 



As we will see later the degree to which t'"^ differs from the intrinsic metric 
r"^ already existant on E, is a measure of the lack of isometricity of the map vr 
and the general energy functional for Skyrme type models measures this non- 
isometricity (an isometry is a map which preserves the metric i.e. distances 
are left unchanged by the map). Intuitively the mapping n produces a strain 
and the energy functional is a measure of the energy attributed to this strain. 
If 



at every point of S then the map is an isometry. 

There are in fact four ways of testing whether or not the map is an isometry, 
in the event vr is invertible. In the following paragraph we will make precise 
the notions of push-forward, pull-back of the inverse metrics and metrics re- 
spectively of the spaces S and S. 

Any map vr between manifolds S and S defines a map vr* called the push- 
forward between the corresponding tangent spaces and vr* called the pull-back 
between their dual spaces (see figure |^). We remind the reader that the dual 
space, sometimes called the space of differential forms, is simply the space of 
real valued linear functions of the tangent space. This means if f G Tp and 
a; e Tp (Tp is just notation for the dual space), u; is a linear function taking v 
to IR, and we write {lj, v) E which is also called the contraction of u) with 
V. The mapping tt is defined as follows: 



p f— > 7r(p). (83) 
The push-forward it* is defined by 

7r*:Tp(^)^T.(p)(S) (84) 
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Fig. 3. Diagram of the push-forward and pull-back for the mapping tt : 5" — > S. The 
same construction is possible for the inverse mapping tt^^ : H ^ S . 

where Tp{S) is the tangent space at the point p to the manifold S while 
Tt^(j,-){Y1) is the tangent space at the point 7i{p) to the manifold E. A vector 
V e Tp{S) with components with respect to the coordinate basis {d/dp'} 

- = (85) 



is mapped to the vector ■n*{v) on S given in terms of the coordinate basis 
{a/^Tr"}, 

7r*(^;) = v'"^ (86) 

^ ' dp' ^TT" ^ ' 



hence the components transform as 



The pull-back tt* is defined as follows: 

7r.:T;(p)(S)^T;(5) (88) 
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where r*(p)(E) is the dual space to T^(p)(E) i.e. the space of differential 1- 
forms, while T*{S) is dual to Tp(iS'). We remind the reader that the coordinate 
basis of the cotangent space is defined with respect to the coordinate basis of 
the tangent space by the condition ( dvr", d/dn^) = 6p. A 1-form lo G T^(p)(S) 
with components cua in the dual coordinate basis {d7r"}, 

uj^uja dvr" (89) 



is mapped to the 1-form 7r*(a;) of T*{S) with dual coordinate basis { dp^} 
n^{uj) ^uja dp' (90) 



hence the components transform as 



These transformation rules generalize tensorially on tensor products of the 
tangent and dual tangent spaces. Evidently, if tt is invertible then 

TT-i : S ^ 5 (92) 



defines a push-forward and a pull-back in the opposite sense. We now come 
to the point of computing images of the metric and the inverse metric under 
TT and under its inverse (if it exists). 

Starting with the metric t on S and r on E we have 



t = iij dp' (g) dpi (93) 
r = Tap d7r" ® dTT^ (94) 

the corresponding pull-backs via tt = 7r(p) for r and 7r~^ = p(7r) for t are 

t' = dp' ® dpi e t;{s) ® t;{s) (95) 



and 



e 7! 



;*(p)(E)®r;(p)(E). 



(96) 
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Also for the inverse metrics 

= r"^ A ^ ^ (98) 



we have the corresponding push- forwards via vr for t ^ and vr ^ for r ^, 



and 



--ly = ^ ^ ^ GTp(5)®Tp(5). (100) 



Now the condition that the map is an isometry is given by any of the following 
statements: 



t' = r 






(101) 


t' = t 




dp^ dp^ 


(102) 


it'')' = r-' 




fi.l ^"P 

dpi dpi 


(103) 






dpi _ 


(104) 



If any one of these equations is true, they are all true. Equations (|101| ) and 
(|104|) require that tt is invertible to make sense. As we will see, it is not 



necessary for vr to be invertible, hence ( |102| ) and ( p.03|) are more fundamental. 
They are all algebraically identical when vr is invertible. We take (|102|) as the 
fundamental relation imposing (locally at the point p) isometry, since ( |102| ) 
does not require the inverse mapping to exist. (|103|) is equally suitable. 

We give an explicit example of the preceding formalism with the SU(2) Skyrme 
field defined on 5*^. We take 5, the initial spatial manifold, to be a topological 
and metrical of radius L. The target manifold is the manifold of the group 
SU(2) which also happens to be an S^. Group manifolds come equipped with 
a natural metric, the so-called Haar measure, and this gives a natural radius 
of one to the target 3-sphere. 
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With the Cartesian coordinates , X'^ , , X'^ on H^, we define a 3-sphere 
of radius L embedded in by the constraint 

3 

Y,X'X' + {Xy = L^. (105) 

i=l 



With the relation 

= ^ X' (106) 



where i = 1, 2, 3 we effect the stereographic projection to IR'^ and obtain the 
following metric 



— — y dx' dx' (107) 

(r74L2 + l)2^ ^ ' 



i=l 



where = I]f=i x^x^. Thus we see that a stereographic projection is simply 
a conformal transformation of flat space {i.e. the metric only changes by an 
overall, space dependent, scaling). 

The natural metric on the target manifold is best expressed in terms of the 
left-invariant 1-forms. These are a natural basis of the co-tangent space of 
SU(2). If tt" [a = 1,2,3) is any set of local coordinates on SU(2), they are 
deflned by the following generalisation of (|T3p: 

jr, = _^tr(r'^f/t(7r)^f/(7r))rf7r<^ (108) 



where t'^ are the Pauli matrices and U taken for convenience to be in the 
fundamental representation of SU(2) (our notation is consistent for this sub- 
section, later we will revert to the original notation of section 2). These are 
left invariants since under the transformation 

U ^VU (109) 



where is a constant element of SU(2), are invariant. (One can also deflne 
right invariant 1-forms by exchanging the role of U and W as in (p!4D ). 

For example let us take the coordinates {vr} deflned by 

U = Vl — vr^ + i TT ■ r = ttq + i vf ■ r. (110) 
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Note that U is a function of vr^, vr^, vr^ and as such, is not a covariant ex- 
pression, in the tensorial sense. Indeed, tt" are also just coordinates and do 
not transform tensorially either. Hence the expression for in the specific 
coordinates chosen does not appear as a tensorial expression. We find: 



1 



3 



2 ^ 



3 

/3= 



(111) 



which are the well known Maurer-Cartan forms p6|] written in this coordinate 
system. 



This set of 1-forms is natural since, first of all, at the identity where tt = 0, 
£^ = dnf". (112) 



Secondly we can obtain at any other point in the group via the pull-back 
of an appropriate mapping defined using the group multiplication. Consider 
a general element Vq in the group with corresponding coordinates vfvo- The 
mapping of a neighborhood Vvo of vfvo to a neighborhood Vi of the identity 
I is furnished by multiplication by VJ^^. If a general element of Vvo is noted 
by V with coordinates vfy, and a general element of Vi is noted by U with 
coordinates vf the map 



: Vvo ^ Vi 



1 - 7?^^, - i TTvo • r)(V 1 - Tfv - i TTv ■ r) 



= \/l - vf(7fvo)2 + i7r(7rvo) ■ r 



(113) 



which gives 



TTlTTv) 



1 - TT^ TTvo + Vl - ^Vn + TTVo X TTy 



(114) 



This induces the pull-back of dvr^ = (C'^I^^q) dvr", clearly where = S^, 



dvr^ 

7f=0 



'-'a 



7f=0 dlTy 



d4 



a4 



7rv=7rvQ 

/3 



dvr'^ 
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which is exactly as we had found before in equation ( |1 1 ID , with vr = vrvo- 
The metric on SU{2) is given by 

T = 6^,C>'(g)C'' (116) 

which clearly indicates the orthonormality of knowing that £^ are invert- 
ible as matrices at each point in the group. A short calculation then shows 
that this implies the metric 



r = -tr 



(t"W^£"wf?)d."sd.". (117) 



Now we return to our setting where tt'^ the coordinates on the group 3-sphere 
are functions of the coordinates on the spatial 3-sphere, since we consider 
a mapping (that we call vr) between these two spaces. The pull-back of the 
metric to the spatial via the mapping is given by 



r' = 7r,(r) = r,^ dx' dx^ 

= -tr UUtt) ^ ^ ^ UUtt) ^ \ K . ^ . ds'(g)dxJ 



V Svr" diff^ J dx^ dx^ 

= -ti[W—W—] dx'®dx^ 
\ ax' axJ / 

= r-^?^ f ^'"^ + (^"^"^ - ^'^"^)') ^ ® (118) 
1 — TT^ ax' V / axJ 

where the last line is relevant to the coordinate system chosen on the group. 

The kinetic term of the Skyrme Lagrangian is obtained by contracting r' with 
the inverse metric on the spatial manifold 

r'^(i + .V4rV*'>^«|j. (119) 



[X^,'IM:'fY^,,U^UU^Uj]. (120) 



5x' dx 

Including the volume measure ^Jg = 1/(1 + j^L?^ 



i=l 

2\3 
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which has the correct hmit as L — > +00. 



3.1.3 Non-linear elasticity theory on a curved space and the Skyrme model 

Now we make the connection with the non-hnear elasticity theory that we have 
treated previously, except generalized to a curved space. Hence we suppose 
that are coordinates on a curved space and that the particle p of the body 
at the point is mapped to the point — The metric at the initial 

point is tij(Xp) while at the image point it is rij(a;p) = tij{x^). The initial triad 
is ejjj(Xp) while the final triad is Cmi^p) = ^mi^p)- The initial orthonormal 
triad is mapped to the final triad in the following way: 

era = ein(^')^ - "^^^'^l^^^ = (122) 



where ^i^{x^) is the inverse dreibein and = d/dx^. Hence the analog of 
the tensor A^^ is 



which is the local Jacobian of the deformation. 
The right-Cauchy-Green deformation tensor is 

i^mm' = A^Al^, (124) 

and the ordinary strain tensor is given by 

E^l{D-I). (125) 

The invariants are defined in the same way as before. 

Now we go to the completely general situation where we are mapping between 
two different curved spaces. The deformation matrix now generalizes to 

- eL^Q (126) 



where the inverse dreibein is defined by the orthonormal basis = d/ Ott" 
in the space tangent to E, while the equivalent strain tensor is 



25 



3 



^mir-i^eJ^r,^. (127) 



Finally, taking the first invariant of D by tracing over m and n gives 



-(l + rV4L^)^Etr(t/tAf/f/t^f;) (128) 
Xl + Xl + \l (129) 



Equation ( |128|) is clearly the usual kinetic term of the Skyrme model La- 
grangian (in the limit where L +oo). 

The Skyrme term is obtained from the curvature tensor defined on the group 
manifold, which is pulled back to the space manifold and then contracted twice 
with the inverse metric there. The curvature tensor is most efficiently defined 
via the machinery of the exterior algebra and the spin connection. The spin 
connection is a 1-form 

ul = u:,drr-^u;,C^, (130) 

which satisfies the structure equation 

dC^ + uj';^ A = (131) 

and the "metricity" condition 

uj^x = -^xv (132) 



where A is the wedge or exterior product (which is simply the antisymmetrized 
tensor product of the forms in question), and uj„\ is the spin connection of 
equation (|130|) with index lowered by 5^^,^. The conditions ( |130|) and (|131|) are 



exactly equivalent to the conditions in the usual formulation of differential 



26 



geometry that there is no torsion (the Christoffel symbol is symmetric in its 
lower two indices) and the metric is covariantly conserved (metricity). 



The curvature is then given by the 2-form 

K = Rt,^^C^ACP (133) 
= dcu;; + cu^ A uj;. (134) 

The spin connection on the group manifold of the target space SU(2) is well 
known |]26| and is given by 



i^^lv = ef,uxC . (135) 



Then using the relation 

d£^ = ej^ C^" A (136) 



which is a little tedious to verify, a short calculation shows 

Kxp = ^'u. (137) 

The pull-back of the tensor R^vXp (pull-backs are only defined for co-tangent 
space tensors) is then 



_ Stt^ Svr^ STr-^ ^ ^ ^ 



We remark that the derivatives d'n'^jd'p^ etc. serve only to change the vari- 
ables from group manifold coordinates vr" to spatial coordinates p\ while the 
Maurer-Cartan forms contain the essential structure. Application of the fol- 
lowing Fierz identities 



r^kV'^ = -^(r^j5ii + r^4j) (139) 
r^vx = -\rlr^,'r\bAv (140) 

yields 

= 2 ix(^^d-^, U^d-p\^'^d^U, U^dxV'\^ (141) 
and contracting twice with the inverse metric on space yields 
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ijki 

(I 



[l + r^lAL''Y tT[[U^dJJ,U^d;Uf] (142) 

i,j=l 



which is obviously the Skyrme term. 

A more geometric and generally valid interpretation is obtained by (for any 
manifold S and S) considering the squared norm of the pull-back of the area 
element defined by two dual basis vectors in the target space The area 
element defined by the two dual basis vectors is 



A^" = (^'AC = C^Q dTr" A dvr^ 

= l{C!iQ-CK) dvr-Advr'^. (143) 
Its pull-back is given by 



A'"'' = A'f dp' A dp) 

1 / \ Stt" Stt^ 

Its squared norm is 



7^ ( CQ - C,;C ) — — dp' A dpi. (144) 



\An^- = it'4^'(ca - ac) (c?a - tfc)|^^|^|^ ("5) 

and expressing = ej^el^^™", summing over /i and u and using the definition 



of the deformation matrix Jj^ = e[^{d7r'^ / dp')(l^ (see equation (|126|) gives 



^|A'^n' = (tr[JJ^])'-tr [JJ^ 



tTi2 

I =yiL[jj J J — iL [J ^ 
= iti[D])^ - ti[D] 

i2\2 I \2\2 I \2\2 



= XiXt + XiXi + XiXi. (146) 

This expression is completely general, allowing for any spatial and target man- 
ifold. Specializing again to the case of and SU(2) we obtain 

A^"" = e^'^'A^ (147) 



where 

= e.^,£^£^ dvr" A dvr'^. (148) 
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Then 

dpi dpi 



A'^ = e,^,Xn^^ dp' A dp^ (149) 



and 



\A.^ I — ^a^u^aXpi CJ^a^lir-^Cf^- 



/3 7 <5 Qp Qpk dpi 

= R' (150) 
by ( |IH7| ) and ( CT - We easily verify that 

This expression for R' and the Skyrme term is in fact identical to that given in 
MantonpS), however there is a slight formal difference. We have pulled-back 



the area elements from the target space (SU(2)) to the spatial manifold and 
computed the sum of their squared norms there. Manton||2^ takes the area 
elements in the spatial manifold and pushes forward their dual area elements 
(tangent space tensors) to the target space and computes their squared norm 
in the target space. This gives the same energy functional. 

To complete our treatment of this example of mapped to SU(2), we show 
the interpretation of the third invariant. Consider the integral coming from 
the third invariant 



J ^3 Vdett d^p = J VdetD Vdett d^p 
s s 

= J A1A2A3 Vdett d 



: J det J Vdett d 



s 

= J detf-^j detC d^p 
s 

= (deg7r) J VdetT d\ (152) 
s 

where the factor deg vr counts the number of times that the mapping tt wraps 
the initial manifold over the target manifold. In reality the last equation is only 
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valid locally on S, the integral gives the volume of the region covered in S. This 
volume must be counted with the appropriate sign depending on whether the 
relative orientation is preserved. Between regions where the relative orientation 
changes sign is a zero of at least one of the A's. This gives rise to natural 
boundaries which should be considered since the sign does not change within 
these regions. Then summing up the volumes of the regions of S with the 
corresponding sign gives exactly the degree of the mapping vr, i.e. the number 
of complete covering of S that the mapping tt provides, multiplied by the 
volume of S. We assume that S* is a manifold without a boundary hence the 
mapping vr must cover S an integral number of times. 



From the previous equation, we have the integral 
1 



degTT 



Vol S 



det(^-^] detC d^p 



(153) 



which can be expressed in terms of U{p) to re-obtain the usual form of the 
baryon number (pf) 



deg. = ^/det(£)det(|) d 



27r 



dp, 



p 



27r2 V 2 



det tr 



1 



^ , det ( tr 
167r2 J V 
s 



1 



167r2 



^fiux tr 



dp\ 
^ dp\ 



dU_ 

d^p 



tr 



dp^ 



^ dp 



tr 



r' 



dU 



,ijk ^3 



p 



247r2 



tTp^diU U^djU U^di,U] e'j"^ d^p 



(154) 



which is the familiar form of the baryon number in the Skyrme model. 



The Skyrme energy function, according to Manton's formalism [pq] , can be 
expressed as 



E = J Vd^ d^p {Xl + + + A^A^ + A^A^ + A^A^; 



(155) 



which is easily re-expressed as 
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E = J Vd^t d^p [(Ai ± A2A3)' + (A2 ± AsAi)^ + (A3 ± A1A2)'] 
s 

TqJ Vdett d^pAiA2A3. (156) 
s 

This shows a novel way of demonstrating the Bogomolnyi bound: evidently 



E>6 



\1\2hVdett d^p 



(157) 



The equality is attained only if (for winding number +1) 
Ai = A2A3 A2 = A3A1 A3 = AiA2. 



(158) 



This system has only three distinct solutions: 

(Ai,A2, A3) : (0,0,0), (1,1,1), (1, -1, -1) + 2 permutations. (159) 



The trivial solution corresponds to mapping 5* to a single point in S and is the 
usual vacuum solution. The third set of solutions is equivalent to the second 
solution after a rotation by 180° in the cotangent space of the target manifold 
about a fixed axis. The second solution implies that the map is everywhere 
an isometry, i.e. , the two 5*^ have the same radius, 1. This shows that for an 
infinite initial sphere, which corresponds to the case of IR'^, the Bogomolnyi 
bound is not saturated and, as is well known, the map is far from the identity 
map. Manton and Ruback|27] and Manton[28| show that as the radius of the 



initial 3-sphere decreases, the map attains the form of the identity for a radius 
of V2. For more details and further applications we refer the reader to the 
hteraturefB7|,BB,B9|,H5|g6|l . 



We close this section with a few general words on this formalism. First of 
all the A's are not independent dynamical variables. Infinitesimal arbitrary 
perturbations are allowed, however integrating to finite deformations is subject 
to consistency conditions. For example there is no deformation of a given 
configuration which can yield Ai = 1 over a finite region, if the region is not 
iso-metric to a part of the target manifold. Actually a smooth mapping tt will 
always give rise to a smooth set of A's. A configuration with a discontinuous 
set of A's is not attainable even though the corresponding energy integral is 
finite. 



We have also made the intuitive paradigm that the Jacobian matrix is a mea- 
sure of the deformation, and hence of the energy. However this is somewhat 
misleading since A; = is clearly a very deformed situation all the same cor- 
responding to zero energy density. For a physical elastic body, Ai = is in 
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fact an infinite-energy deformation hence the corresponding energy functional 
is not at all like the Skyrme energy functional. Very schematically, the energy 
density of an elastic body is 

eel.body ~ P - /)' (160) 

while in the Skyrme model it is like 

esk. -.D + D^. (161) 



Thus the Skyrme ground state is around D = which is quite unlike the case 
of the elastic body, where D = 1. 

This completes our exposition of the interpretation of a field theory, specifically 
the Skyrme model, as a non-linear elasticity theory. In the next two subsections 
we will look in more detail at the Skyrme model and its static, low energy 
configurations. First we will elaborate on the instanton method for obtaining 
an analytical ansatz for the set of relevant low energy configurations, and 
second we will describe the use of rational maps to obtain reasonable ansatze 
for multi-baryonic minima. 



3.2 The instanton method 



The instanton method uses the known solutions of 4-dimensional Euclidian 
Yang- Mills theory called instantons|47] and their moduli spaces to obtain 



Skyrme field configurations. The relation between the two seems tenuous at 
first, however, the known global topology and symmetries of the instanton 
moduli space and its similarities to expected properties of low-energy Skyrme 
field configurations seems to point in that direction. Consider the case of 
B = 2, here know for two widely separated Skyrmions, there are 12 indepen- 
dent degrees of freedom. We expect the relevant low-energy space of configura- 
tions to also have 12 dimensions. Mantonj^ proposed that this sub-manifold 



could be obtained as the union of all gradient flow curves linking together all 
the low-energy critical points. We will return to this subject in much detail in 
section 4. 

This idea to obtain the 12 dimensional sub- manifold, Aiu which should serve 
as the correct truncation of the full field theory description of the interactions 
and dynamics of two Skyrmions, from gradient flow curves although in prin- 
ciple sound is in practice only numerically, approximately implementable. A 
consideration of the symmetries involved, led Atiyah and Mantonl^O] to sug- 



gest that an analytical construction of a manifold, which might be a reasonable 
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approximation to the true sub-manifold, could be obtained from certain in- 
stanton configurations and their holonomies. 



The observation consists of the following two steps. First of all, from any SU(2) 
instanton configuration in IR^, it is possible to obtain a unitary matrix valued 
field defined on IR^ by 



where P denotes the path ordered integral. We will show later, how the baryon 
number of U{x) is equal to the instanton number of r). Secondly it has 
been known that the space of configurations of two instantons interpolates 
continuously and smoothly from an axially symmetric, localized configura- 
tion to two, individual, "spherically" symmetric, well (infinitely) separated 
instanton configurations. It is a matter of verification that the corresponding 
Skyrmion holonomies interpolate smoothly between the toroidal lowest energy 
deuteron to two infinitely separated single Skyrmions. It is also possible to ob- 
tain configurations which correspond to the spherically symmetric dibaryon 
type configuration. 

The most vexxing problem is that the two instanton configurations are 16 
dimensional, that is they have 16 independent parameters. One of these cor- 
responds to a global time translation, the integral over the time direction 
removes this degree of freedom, leaving 15. This is larger than the 12 dimen- 
sional manifold which is being sought. It is not evident what is the proper 
way to reduce the number of parameters by three. Nominally one should re- 
implement the gradient flow method on this sub-manifold of the configuration 
space starting from the highest energy critical point of the Skyrme energy 
functional restricted to the sub-manifold. Such a calculation has not been ef- 
fected, numerically it is just as difficult to work with a discretized version of 
the full problem rather than the one defined on the sub-manifold, hence there 
does not seem to be a compelling motivation to study the gradient flow here. 
The problem has been studied in detail for the case of the most attractive 
channel, which we will return to a little later. The instanton method for three 
and higher baryons is not very efficient. 

3.2.1 Topological numbers 

First let us solidify the connection between baryon number and instanton num- 
ber. An instanton configuration with instanton number k is technically defined 
as a connection on an SU(2) principal bundle over the four-sphere S"^ with sec- 
ond Chern number C2 = k^^^. The relationship between a gauge field defined 




—00 



(162) 
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on and S"^ is obtained via stereographic projection. The Tisiial instanton 
configurations which satisfy the Yang-Mills equations of motion defined on 
have non zero field strength F^^, in a localized region of space-time, and 
achieve a pure gauge field type configuration towards Euclidian infinity. This 
means that on the manifold at infinity of H^, which is topologically S^, an 
instanton configuration defines a smooth group element valued configuration 
U {x, r) I (s,t)-»oo and the gauge field is given by 

= Ud^Ul (163) 
The second Chern number corresponds to the integral 



A; = ^/d^a; e'^'^^^tr(F^.F,,) 



-1 
1 

32^ J " ■""'^ 
1 



327r2 

oo 



da^Ki' (164) 



where 



^ _2 e/^-^- tr(^,FA, - '^-A^A^A,) (165) 



and 

F^u = d^A, - d,A^ + [A^, A,]. (166) 



Replacing A^ = Ud^W gives 
4 



K'* = - e'^'^^^tr(^[/a,[/^ UdxU^ Ud^U^^ , (167) 

and k is given by the integral 

k = W{U) = j da^ e"""^" ii[ud^U^ UdxU^ Ud^U^y (168) 



This is exactly the measure of the winding number of the mapping of — > 
SU(2) defined by the group element at infinity. These mappings define elements 
of the homotopy group 

n3(SU(2)) = 7L, (169) 
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the integer corresponding to the Chern number. The configuration on does 
not attain a constant value at infinity, hence it is difficult to interpret this field 
as corresponding to a field on after stereographic projection. (To be precise, 
the stereographic projection involved is actually the conformal mapping taking 
5*^ — > IR''. Qfj^u = l/(|a;p/4i?^ + 1)^ r/^i, is the conformal transformation taking 
IR^ to S*^ with the coordinate transformation exactly as in equation (|106|) . This 



mapping has the advantage of mapping solutions of the equations of motion 
on one manifold to solutions on the other because the Lagrangian of Yang- 
Mills theories is conformally invariant ||48||.) However the configuration on S*^ 
is also somewhat subtle, is actually a connection on a non-trivial SU(2) 
principal bundle defined over the base manifold S^. These bundles are specified 
by fixing the transition function which maps the fibre SU(2) = over the 
"northern hemisphere" of S"^ to the over the "southern hemisphere" at the 
"equator". The equator of S"^ is simply an S^. Hence the transition functions 
are tantamount to defining a group element over the equatorial S^. This means 
that one defines a mapping of the equatorial 5*^ to the group SU(2) = 5''^. Such 
mappings fall into the disjoint homotopy classes labelled exactly as in (|169| ). 



The second Chern number of the bundle ( |164D is exactly equal to the integer 



characterizing the homotopy class of the transition function. 

Geometrically there is no constraint on the size of the coordinate charts; one 
need not restrict oneself to equal hemispheres. There is no hindrance to ex- 
tending the southern hemisphere to include the whole of S"^ except for one 
point, the north pole. Indeed, in this way we will extend the solution of the 
equations of motion to almost everywhere on S"^. The solution will be singular 
at the north pole however only for the connection. The integral (|164|) will be 



an integral over the whole sphere of only the field strengths associated with the 
connection and the density ii{F^jjF^^) is non-singular over the whole sphere. 
The result must still give the second Chern number. Now the conformal pro- 
jection of this field configuration to IR^ will give the field configuration that 
satisfies the equations of motion on IR^ and furthermore whose integral cor- 
responding to ( |164| ) is also C2{A^) (since ( |164| ) is independent of the metric). 



The values achieved at the north pole by defined by the limiting value 
of the configuration along any path leading to the north pole are all equal 
modulo gauge transformations. The conformal transformation maps this field 
on to a configuration on which becomes a pure gauge configuration at 
infinity since the field strength at the north pole is diluted over the entire 5*^ 
manifold at infinity. 



The path ordered integral ( [L62| ) which defines U{x) starts at r = — oo and 
follows a straight line to r = +oo. This corresponds to a curve on 5*^ which 
starts at the north pole, follows a particular path on the S"^ and returns back to 
the north pole (see figure The path on is simply given by the intersection 
of the S"^ with a 2-plane, that which is defined by the line of integration in IR^ 
and the north pole. This intersection is actually just a circle. The set of such 
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Fig. 4. Diagram of the stereographic projection from S'^ to M^. 

curves on S'^ leave the north pole on one side of a 3-dimensional hyper surface, 
circle around the S^, and return to the north pole from the other side of the 
hypersurface. The hypersurface is in fact just a "great" 3-sphere, exactly like 
the equatorial except that it passes from north pole to south pole, and 
back. It is exactly the inverse (stereographic) projection of the spatial IR^ (of 
the R"^) onto the 5*^. The curves leave the north pole, intersect this great 
exactly once and come back to the north pole in a symmetric fashion. Hence 
each curve on the 5^ defines a unitary matrix valued configuration on the great 
at the point where the curve intersects this great S^. This configuration 
is by construction continuous. Consequently we manage to define a winding 
number ( |169| ). 

The winding number is invariant under any continuous deformation which 
keeps a one to one relation between each curve and the points of the hyper- 
surface. Envisage the following deformations. The lines of integration are well 
represented by the lines of forces emanating from and returning to an ideal 
pointlike "dipole" situated at the north pole. The great sphere separates these 
lines of force at the north pole into outgoing lines on one side, and incoming 
lines on the other side. We simply imagine moving the two charges comprising 
the dipole apart, keeping the lines of integration the same as the lines of forces 
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leaving the positive charge on one side and arriving at the negative charge on 
the other. Such a modification of the fines of integration wifi result in a fio- 
motopy of tfie original U{x), and fience will not cfiange tfie fiomotopy type. 
Finally we will arrive at tfie situation wfiere tfie two cfiarges occupy antipo- 
dal points of tfie four spfiere (actually on tfie equator) and tfie lines of force 
emanate symmetrically from one cfiarge, cross tfie great 3- spfiere and finally 
terminate on tfie opposite cfiarge at tfie antipodal point. We make one furtfier 
fiomotopy, we rigidly rotate tfie system of cfiarges, lines of integration and tfie 
great 5*^ until tfiey are vertical, sucfi tfiat tfie positive cfiarge is at tfie soutfi 
pole and tfie negative cfiarge is at tfie nortfi pole, and tfie usual equator now 
corresponds to tfie great S^. Sucfi a deformation requires a simultaneous redef- 
inition of tfie stereograpfiic projection, and a deformation of tfie U{x) (because 
tfie integration lines are cfianging) but it is clearly a continuous deformation 
keeping tfie fiomotopy type invariant. Finally to be complete, we fiad started 
witfi closed line integrals originally (leaving and returning to tfie nortfi pole) 
but now we fiave open line integrals starting at tfie soutfi pole and terminating 
at tfie nortfi pole. We can easily remedy tfiis by adding one patfi to all of tfie 
otfiers, starting at tfie nortfi pole and descending down a fixed meridian to tfie 
soutfi pole for all of tfie line integrals. Tfiis simply left-multiplies eacfi U{x) by 
a constant unitary matrix, wfiicfi again does not modify tfie fiomotopy type. 
Closed contour patfi ordered exponential integrals are gauge covariant fience 
we can return to tfie original description of tfie instanton on S"^ witfi more tfian 
two patcfies. Now we fiave arrived at tfie starting point of tfie demonstration 
given in Manton and Atiyafi [pl|] wfiere tfiey sfiow tfiat tfie winding number 
of tfiis configuration is tfie same as tfie instanton number (tfie second Cfiern 
number of tfie instanton bundle in question). Tfiis demonstration proceeds as 
follows. 

Tfie first step is to use tfie gauge freedom to put tfie gauge field in tfie specific 
gauge wfiere tfie component of tfie gauge field along tfie meridional directions 
vanisfies. Tfiis can be establisfied in eacfi patcfi separately. Tfien tfie integral 

Pexp(- / A^dx^") = 1 (170) 



since tfie inner product A^dx^ is zero along tfie meridional patfi. However 
tfie definition of tfie patfi ordered exponential, wfien tfie patfi of integration 
crosses a boundary between patcfies is sucfi tfiat one must multiply tfie contri- 
bution coming from tfie first leg of tfie patfi by tfie transition function before 
continuing witfi tfie integral in tfie second patcfi (see figure |^) . 

1 1 TO 

Pexp ( - y A^dx^") = Pexp{- J A^dx>')G{x)P exp ( - J A^dx^"). (171) 

TO 



37 



Region I 



Region // 



r 



r = ro 



r = 1 



transition function 

G{x) 



Fig. 5. Diagram of the role of the transition function G{x) between regions / and 
II. 

Such an expansion is gauge covariant under simultaneous independent gauge 
transformations in the patches / and // since the transition function is defined 
to transform as 



G{x) UuG{x)ul 



(172) 



while 



Pexp(- J A^dxn^Ui{To)Pexp{- J A^dx'')Ul{0) 



1 1 
Pexp ( - y A^dx^) ^ C/„(l)Pexp ( - J A^dx^)Ul,{ro) 



yielding 



TO 



TO 



(173) 
(174) 



1 1 

Pexp{- J Ai.dx") t/ii(l)Pexp ( - J ^^dx'^)t//(0) (175) 
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Now since A^dx^ = in our gauge, we obtain 

U{x) = UoG{x) (176) 

where G{x) is the transition function at the equator. Then, the baryon number 
B = W{U{x)) = W{G{x)) = C2{A^) = k- (177) 

hence the baryon number and the instanton number are identical 

3.2.2 The sector B = l 

For k = 1 the instanton profile is given by 

A^ = ia^^dy\\ip (178) 

where the definition of a can be found in the article of Jackiw, Nohl and Rebbi 
(JNR)[^], however the time component is explicitly 

^4 = i^-r (179) 



2 p 



with 



This is in the 't Hooft gaugeH- The singularities in p are gauge artefacts and 



hence do not contribute to gauge covariant quantities such as the definition 
of the Skyrme field. A local gauge transformation moves the singularity to 
wherever we want, but of course does not affect the U{x). We will not show 
this here. The function p ( |180| ) has 5 parameters, but three others are to 
be added in because of global gauge transformations, which were factored 
out in the definition of instantons. This yields 8 parameters. Evidently the 
4 translation parameters X'^ are the center of mass coordinates, the spatial 
ones determine the spatial center of mass of the corresponding Skyrmion, the 
temporal one is absorbed by the integration in the time direction. Factoring 
out these four leaves 4 parameters, A governing the overall scale, and three 
corresponding to global gauge transformations. Hence 

P=1 + TT^^ (181) 
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where r'^ = x ■ x. Then 

i X 2x ■ 1 

2p(r2 + r2)2 

and 



i X 2x ■ T , ^ 

A^ = -i;- ,.., . (182) 



C/(x)=Pexp -/-i- 

^ — OO 

oo 

= exp (^iXx ■ T J 



oo ^ 

I X 2x ■ T 



p (r^ + r2)2 
(ir 



r2 + r2)2 + A(r2 + r2). 



Thus 



/(r) = 7r(l-(l + A)-i) (184) 

which satisfies /(O) = vr. This point is actually determined by the limit from 
non-zero values of r, since the integral above is not well defined for r = 0. 
The singularity is gauge dependent, as mentioned above, hence by a local 
gauge transformation we can move the singularity away from r = without 
affecting the value for the Skyrme field. The Skyrme field so obtained will be 
continuous at r = 0, hence the value at r = can equally well be defined 
as the value obtained from the limit of non-zero r. Furthermore, f{oo) = 0. 
For minimum energy one finds A = 2.109, and the corresponding energy is 
E = 1.2432 X 127r^. This exceeds the numerically obtained minimum energy 
solution by only 1%. 

Adding in gauge transformations, i.e. global iso-rotations simply combs the 
Skyrmion profile without affecting the energy. The iso-rotation parameters 
are already evident in the JNR[^ parametrization of the instantons. Here we 
have 



which has 10 parameters. It is known that several of these are local gauge 
artefacts, indeed, the instanton obtained is gauge equivalent to the 't Hooft 
instanton ( |180| ) with 



-^1 + -^2 
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-^1-^2 I |2 



(Ai + A2)^ 



Xi-X2\\ (187) 



The JNR instanton generates, however, a Skyrme field differing from that 
obtained from the 't Hooft instanton by a global gauge transformation given 
by 

(X2-Xif (X2-X1) ^ , , 

Uo = \^ ^ + i ' ' -r. (188) 



This completes the case k = 1{B = 1). 



3.2.3 The sector B = 2 



In this subsection we present the findings of Atiyah and Manton[^, where 
they relate the parameters (moduli) of /c = 2 instanton configurations intro- 
duced analytically by JNR and geometrically by Hartshorne[Q, to the various 
parameters of the corresponding B = 2 Skyrme configurations. 



For k = 2{B = 2), the JNR parametrization is 

"^1 I -^2 , A3 /1cn^ 



a 15 parameter solution. It is clear that the overall scale of the A's is never a 
parameter, yielding 14 parameters. In addition there is an explicit one param- 
eter family of local gauge transformations included in p, reducing the number 
of true parameters to 13. We will return to this redundancy later. Integrating 
over Euclidean time to obtain the Skyrme field reduces the number of pa- 
rameters to 12. Finally putting in the 3 iso-rotational degrees of freedom, as 
they are not included in the solution, implies that the corresponding Skyrmion 
fields will have a total of 15 parameters. 

In general, for higher k, there is no local gauge transformation in the JNR 
expression for p, which thus has 5fc + 4 parameters. Integrating over time re- 
moves 1 but adding in three for global iso-rotations yields in general 5k + Q 
parameter Skyrme fields. This is obviously not the full complement of 6k that 
we expect for B = k Skyrmions. The full instanton moduli space is actually 
supposed to be 8k dimensional. These include 4 positions, 1 scale and 3 iso- 
rotations per instanton, not removing the 3 overall iso-rotation parameters. 
Integrating to get the Skyrmions removes one parameter, implying an 8/c — 1 
dimensional manifold. This manifold would correspond to 3 positions, 3 isoro- 
tations and 1 scale per Skyrmion, and k — 1 relative "time" coordinates. These 
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5* 



Fig. 6. Diagram of the cut by a plane, exhibiting a triangle and ellipse. 

time coordinates serve simply to fix the order of tlie individual Skyrmion fields 
in the corresponding product ansatz type configuration, if the Skyrmions are 
well separated in these time coordinates. This ordering is of course quite irrel- 
evant if the Skyrmions are well separated spatially. However the interpolation 
between different orderings is quite important when they are spatially close 
together. 



The full 8k dimensional manifold of instanton solutions is well understood 
algebraically but not analytically. The largest manifold of analytically explicit 
solutions corresponds to the 5fc + 4 dimensional manifold of JNR. For k = 2 
these have 13 parameters and for which there is an algebraic characterization 
given by Hartshorne[^. Hartshorne proves that there is a 1-1 correspondence 
between the instanton solutions (on S"^) and a set of ellipses that are interior to 
the S"^ (we imagine that the S*^ in question is embedded in IR^). The ellipses lie 
in a 2-dimensional plane that intersects the S"^ in a (coplanar) circle, subject 
to one condition. The condition insists that the ellipse can be circumscribed 
by a triangle whose vertices lie on the circle (see figure |^). This condition is 
called the Poncelet condition, and he proved that if such a triangle exists. 



then there is a one parameter family of such triangles, called a porism|31 



The Poncelet condition is effectively one condition on the parameters of the 
ellipse, for a fixed circle. An ellipse is parametrized in general by 5 parameters: 
2 give the semi-major and semi-minor axes, 2 fix the position of the center and 
one fixes its orientation in the plane relative to a fixed set of coordinate axes 
(rotations). The Poncelet condition leaves 4 degrees of freedom. Now a 2-plane 
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passing through the origin in is specified by a division of the into the 
2-dimensional space of the plane and the 3-dimensional space orthogonal to 
the plane. The action of the orthogonal group 0(5) generates all the different 
possibilities from any given initial one. However, if the orthogonal group acts 
only in the 2-plane or in the 3-dimensional orthogonal space, then we obtain 
nothing new. Hence the dimension of the space of 2-planes in passing 
through the origin is dim 0(5) — (dim 0(3) + dim 0(2)) = 10 — (3 + 1) = 6, 
where we have used that the dimension of O(A^) is A^(A^ — l)/2. Removing 
the condition that the plane passes through the origin allows 3 translations, 
one for each independent orthogonal direction, yielding 9 degrees of freedom. 
Adding in the four degrees of freedom of the ellipse yields 13 in total. 

Given the circle and ellipse, the corresponding family of triangles is neatly 
described by a cubic equation. Suppose that s = tan(|) is a variable along the 
circle {6 G [— tt, tt] an angular variable around the circle). Then the vertices of 
the triangle correspond to si, 82,83, which are without loss of generality the 
roots of a cubic polynomial equation, 

P08^ + Pl8'^ + P2S + P3 = 0. (190) 



Clearly there are some global constraints that they the PiS must satisfy so that 



the cubic (|190|) has three real roots, but these constraints do not remove any 
degrees of freedom. However, the piS are real, their overall scale is irrelevant, 
and they are not all zero, hence they define a ray in IR^. The space of all rays 
in IR^ is called the real projective space of dimension three, denoted IRP3. The 
Poncelet condition is expressed in this light by requiring that the coefficients, 
of the cubic equations for all the triangles in the porism, must lie along a 
straight (projective) line in IRP3. Hence if 

qos^ + qi8^ + q2S + gs = 0. (191) 



is a cubic whose roots correspond to the vertices of another solution of the 
Poncelet condition, then all solutions are obtained by the interpolation 

(/ipo - i^qojs'^ + (/upi - yqi)s^ + {^p2 - yq2)s + (yups - uq-i) = 0. (192) 

Now if for some /2 and u the coefficient ftpo + uqo = 0, then the surviving 
quadratic equation gives two of the roots, while the last root is pushed off 
to s = ±cxD (which is the same point for either sign). The requirement that 
the ellipse be interior to the circle implies that the roots of ( |192D are real and 
distinct for all /i and u. This last constraint precludes the possibility that also 
jlpi + vqi = since then we get one finite root but two roots get pushed off 
to s = ±00, which means that they are the same point on the circle. This is 
not permitted for an ellipse that is interior to the circle. 
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Fig. 7. Diagram of a triangle showing the triangle and the ellipse it defines, 

touching it at the points Ai, A2 and ^3. 

The JNR parametrization corresponds to instantons in H^. However these are 
related to those defined on S'^ by a conformal transformation. This corresponds 
to a specific stereographic projection of 5''^ to H^. This projection takes circles 
to circles and triangles and ellipses in the interior of the also to triangles 
and ellipses, respectively. The JNR parametrization immediately gives us the 
ellipse and the circle. The points Xi determine a circle and the vertices of a 
triangle. Then we use 

Ai _ XiA^ X2 _ X2A1 A3 _ X3A2 (193) 
A2 ^3^2 -^3 ^1X3 Ai A2X1 



to determine the points A[, which are defined to be on the line joining Xj with 
Xk (with i,j,k all distinct, and where ^1^43 signifies the distance between 
Xi and A^ along the line joining Xi with X2 for example). The Ai's and the 
XiS are part of the JNR parameters, the AiS are uniquely determined by the 
ratios of the Ai's. The Xi's determine the circle uniquely, and the ellipse is 
determined also uniquely by the points Ai along with the condition that the 
ellipse be circumscribed by the triangle (see figure 

The description of the instanton in terms of the triangle and ellipse allows us 
to readily understand the symmetries of the instanton and hence the resulting 
Skyrmion. The Ai's seem to correspond to one triangle in the porism, however 
they in fact afford an interpretation in terms of an infinitesimal variation of 
the triangle within the porism. If X/ are the vertices of an infinitesimally 
close triangle to the one determined by the Xi's, then the line X(X2 is an 
infinitesimal rotation of the line X1X2 about the point of tangency (see figure 
By elementary geometry 
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Fig. 8. Diagram of the infinitesimal rotation of X1X2 to XjXg. Similar rotations 
occur to the other two sides of the triangle. 

XiX( = (XiA3)^ (194) 
X,X'^ = {A,X,)^ (195) 

hence 

XiX[ XiA-i Ai 



X2X2 A^X2 A2 



(196) 



using (|193|) . Thus for an infinitesimal change of the triangle, the vertices move 
along the circle by a distance proportional to the Ai's. This means that the 
JNR data determines a 1 parameter family of cubic equations. With X^X^ = 
dOJ siny^i, and Si = tan ^, then 

dsi = ^(l + sf)d^i. (197) 



The family of cubic equations 

fl (s- Si){s - S2){S - S3) + + sl){s - S2){S - S3) + 

A2(l + sDis - S3)is - Si) + A3(l + sDis - si)is - S2)) = (198) 

has three real roots which correspond to the three angles 6i. Indeed if z/ = 
the roots are Si, S2 and S3. 

For an infinitesimal u, Si = S[ + ds[, solving to first non-trivial order yields 
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PL dsids2ds3 + dp (^Ai(l + s^) ds2ds3 + 

A2(l + si) dsidss + Xs{l + S3) dsids2 ] = (199) 



implying dsi = — ^^Ai(l + sf) which corresponds to the desired variation 

Now we can address the case where the circle on passes through the "north 
pole", the point from which we do the stereographic projection. In this case 
the circle projects to a straight line in IR^ and the ellipse also projects to the 
same line. The projections of the porism of triangles gives a triplet of points 
along the line, coming from the vertices. These triples are again the roots of 
a cubic equation as before, however the parameter is just an affine parameter 
along the line. If Xi and are two infinitesimally separated triplets, we can 
define weights (up to an overall constant) by 

Al ^ XiX[ A2 ^ X2X2 A3 ^ X3X3 ^^gg-j 

A2 X2X2 A3 X3X3 Al XiX[ 



then the JNR potential is as before. Conversely, given the X[ along a line in 
IR"^, we may invert the stereographic projection and reconstruct the instanton 
and its associated circle and elhpse on S"^ equally well. The inverse images 
of the Xi's determine the circle and the vertices of the triangle, the weights 
determine the ellipse. The corresponding line of cubics is defined directly for 
the affine parameter s on the line in IR^. An infinitesimal variation (|200|) 
implies the variation 

dsi = Aids. (201) 



Then the projective line of cubics is 



fj,{s - Si) (S - S2){S - S3) + Z/[Ai(s - S2)(S - S3) + 

A2 {s - ss){s - Si) + A3(s - si)(s - S2)] = (202) 

and an infinitesimal variation du yields dsi = — ^Ai in concord with ( |201D . 
The corresponding points on S*^ projecting to these roots gives the porism of 
triangles there. 

We will next consider two special cases which give rise to interesting B = 2 
Skyrmions. First consider the case of an ellipse of very high eccentricity, such 
that it almost touches the surface of the sphere 5"^ at two points (see figure ^. 
The instanton degenerates to two k = 1 instantons near these points. Stereo- 
graphic projection gives 2 well separated and localized instantons in R'^. For 
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Fig. 9. Triangle corresponding to a pair of well separated instantons or, correpond- 
ingly, to two well separated Skyrmions in the product ansatz. 

simplicity we take Xi not be collinear. With Ai = 1 (without loss of generality) 
and A2, A3 << 1 then for x near X2 the JNR potential is approximately 

1 A2 A3 1 A2 / X 

P = I ^ + I W + I W - I ^ + I W 203 



which is exactly the JNR potential for a single instanton centered at ~ 
X2 with scale parameter A2IX2 — Xip. The scale being proportional to A2, this 
allows for the possibility that the instanton here has little overlap with the 
corresponding one located at X3. The integration along the time lines gives two 
well localized Skyrmions. The minimum energy single Skyrmions are obtained 
for A ~ 2.109, which is not particularly small. However given that the scale 
size of the instantons is quadratically dependent on the separation, it is always 
possible to choose this so that the resulting single Skyrmions are well separated 
and of minimal energy. 

If Xi is relatively near X2 and it changes while keeping everything else fixed, 
the Skyrmion at position X2 will vary over all possible orientations while the 
one at X3 will remain essentially unchanged. (|203|) implies a gauge transfor- 
mation ( [TH5D 



X2 — Xi X2 — Xi 



which covers SO (3), hence all possible orientations, twice as Xi varies over a 3 
sphere centered on X2. Hence all possible relative orientations are permitted 
and this case contains all well separated Skyrmions which gives rise to the 
product ansatz. 
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The second case that we will consider is the spherically symmetric situation, 
which gives rise to hedgehog fields for B = 2. Consider the potential 



Ai A2 At 
\ = \ 

+ (r — Ti)^ + (r — r2)^ + (r — r^y 

N(r,T) , , 

^ ' ' (205) 



nf=i(r2 + (r-r02) 

with = X ■ X. The corresponding JNR instantons are situated along the 
time axis. A^(r, r) is a quartic polynomial in r with r dependent coefficients 
which is positive since the denominators can never simultaneously vanish. 
is actually gauge invariant. has r dependent roots a(r), a*(r), /3(r), /3*(r), 
with the imaginary parts of a and /9, 3(a), Q'(/3) > 0. Then 



In p = ln(r - a) + ln(r - a*) + ln(r - /3) + ln(r - /?*) 

- An (r - Tj - ir) + In (r - tj + \r)\ (206) 

which in turns yields 



1 dp 1 da 1 da* 1 d/3 1 d/5* 



p dr (r — a) dr (r — a*) dr (r — /3) dr (r — /?*) dr 



j^^xT-rj-ir r-rj + ir 



+ -^^-^)- (207) 

Then 



^4 = ^-4^r.r (208) 
2 p dr ^ ' 



so the Skyrme field from the integration along time lines would give 
1 °° 1 r1 



The TT comes from the difference between instantons on S"^ and IR^, the in- 
tegration along time lines in IR^ must be closed with a semi-circle at inffnity 
to give a truly closed integration contour. With this addition the boundary 
condition \J\ao 1 is satisffed. The integral (|209|) may be computed using 
standard techniques of contour integration, yielding 
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/(r) = 7ri-^(a + /3) + 27r 
dr 

= 7ri-^(^^{a) + ^{/3)^ +2tt (210) 

using a + a* + P + P* is independent of r. The roots of a quartic can be 
found in closed form, however this does not elucidate the properties of the 
configuration. 

The JNR potential has 5 parameters, ti and the 2 independent ratios of the Ai. 

can be varied arbitrarily by a transformation moving along the line of cubics 
(changing the porism of triangles), leaving 4 parameters. We will call these Ti, 
T2, Ai and A2 in the equivalent 't Hooft parametrization. (Take A3 = 1, ra = 
and go to the 't Hooft parametrization.) A rigid translation of the Ti's does 
nothing hence we take Ti = — T and T2 = T, leaving 3 parameters. Finally 
imposing time reversal symmetry Ai = A2 = A, we get 

iV(r, r) = (r^ + + 2(A + T^y + 2(A - T^y + T^ + 2AT^ (211) 
which implies that Q{a) = Q{l3) and since this is just a quadratic equation in 

(^^{a)^ = ^iR^ + A-T^ + {y + A + T^) - A2)i) (212) 



r2 



and 



f{r) = 2tt -2Tr-^%a) (213) 



As T ^ 00 we can find the form of a, yielding 

/(r) = 2vr(l-(l + A)-|y (214) 



This has energy E = 1.855535 x 2An'^ when A = 2.6211. Actually the true 
minimum occurs for ^ 84.6 and A = 2.6427, with E = 1.855529 x 247r^ 



which has been shown numerically (see |31|). 
The time centered 't Hooft potential 

A+A A-A 
P-^ + r2 + ^r-Ty + r2 + (T + r)2 ^^^^^ 
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gives an approximate potential 

/ A + A\"^ / A-A\"^ , , 

/(r) = 2vr-vr 1 + — ^ - vr 1 + — ^ (216) 



which is the product of two 5 = 1 hedgehogs of scale parameters A + A and 
A- A. The minimum occurs at A = 0, A = 2.6211 with E = 1.855536 x 247?^. 
The energy is symmetric in A so it is a reasonable conjecture that the minimum 
actually occurs at A = 0. 

The spherically symmetric B = 2 hedgehog is stable against perturbations 
preserving that symmetry. Its unstable modes violates this symmetry and 
span a 6 dimensional vector space. Under 0(3) rotations these decompose 
into two 3-dimensional irreducible sub-spaces, 0(3) is vectorial for one and 
axial for the other. 

In terms of the product ansatz we lower the energy if two coincident Skyrmions 
are displaced - along 3 independent axes (axial mode) or rotated relatively 
in iso-space about 3 independent axes (vectorial mode). There are actually 3 
modes which increase the energy while preserving the hedgehog form, the 3 
independent parameters found above, Ai, A2 and T. There are 6 zero modes 
corresponding to translations and global rotations (equivalently iso-rotations) 
and including the 6 modes which decrease the energy - unstable or nega- 
tive modes, giving a total of 15 modes. Using the JNR parametrization the 
minimum energy hedgehog has 

1 + !^ + ^ (217) 



where A' = ^ and T = (T^ + 2A)^ and with = 84.6, A = 2.6427. 

This configuration corresponds to a line of cubics with coUinear roots. The 
perturbations which break the coUinearity gives the vector instability while 
the perturbations that rotate the line gives the axial one. These perturbations 
of the instanton configuration give rise to exactly the same perturbations of the 
corresponding Skyrmions which reduce the energy according to the analysis of 
Bang and Wirzba[^. We will return to these unstable modes in section 4.3. 

The most attractive channel instantons and hence Skyrmions, in the Hartshorne 
description, are distinguished by concentric circles and ellipses. For high ec- 
centricities we have two well separated Skyrmions with relative iso-rotation of 
180°. The minimum energy configuration appears when the ellipse degener- 
ates to a circle of radius -| (see figure 0). Evidently the porism of triangles is 
given by a set of equilateral triangles, obtained from one another by a simple 
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^3 



Fig. 10. Reflection symmetry: xi — > —xi and X2 —X2, and 0(2) symmetry about 
the xs axis. 




Fig. 11. Reflection symmetry about the three orthgonal ctXGSI X'l ^ X'l- 

rotation. This Hartshorne configuration exudes 0(2) symmetry, which indeed 
the instanton configuration, and the subsequent Skyrmion field also exhibit. 

If the Hartshorne ellipse is only concentric with the circle then the 0(2) sym- 
metry reduces to reflection symmetry with respect to the three axes. The 
Poncelet condition requires that a + b = R, where a and b are the semi-major 
and semi-minor axes respectively while R is the radius of the circle (see figure 



ni) . When the ellipse becomes extremely eccentric (but always remains con- 



centric with the circle), it is easy to verify that the configuration corresponds 
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Fig. 12. Hartshorne diagram corresponding to two well separated Skyrmions in the 
most attractive channel. 

to well separated Skyrmions with a relative iso- rotation of 180°. Indeed the tri- 
angle degenerates to a right isoceles triangle with hypotenuse approximately a 
diameter (see figure |T2|). The configuration corresponds to Skyrmions localized 



at positions X2 and X3 with iso-rotation factors {X2 — Xi)-r and (X3 — Xi)-r 
respectively ((X2 — Xi) is a unit vector in the X2 — Xi direction). The con- 
figuration is 

U= (^i(X2^Xi) ■ T^Uix - X2)i(X2^Xi) • f 

X (i(X3^Xi) ■ r) V(x - X3)i(X3^Xi) ■ r (218) 
where U{x — X2) is a Skyrmion at X2, hence the relative iso-orientation is 

i? = i(X2^Xi) ■ f(-i)(X3^Xi) • f 

= if(X2^Xi) X (X3^Xi)) ■ r (219) 



for the Skyrmion located at X3. However the vector 

(X2^Xi) X (X3^Xi) (220) 

is a unit vector (because the angle at Xi is 90°), orthogonal to the plane 
defined by Xi, X2 and X3, and hence also orthogonal to the separation axis 
X2 — X3. Furthermore since R is of the form i n ■ r it effects an iso-rotation by 
exactly 180° as desired, about the axis (X2 — Xi) x (X3 — Xi). The discrete 
symmetries of the concentric ellipse expand to the continuous 0(2) symmetry 
of the concentric circles and implies that this configuration should correspond 
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to the minimum energy configuration which is known to have toroidal sym- 
metry p2| , ^ ,^|. After minimization with respect to the (1) free parameter R 
one actually finds the minimal torus to within a few percents. 



The manifold of attractive channel Skyrme fields forms an 11 dimensional sub- 
manifold of the full 15 dimensional manifold of instanton generated Skyrme 
fields. The ellipse simply has 2 degrees of freedom, its eccentricity and its 
orientation, instead of 4, and the orientation of the 2-plane must be orthog- 
onal to the time axis implying only the action of 0(4)/(0(2) x 0(2)) yield- 
ing 6 — (1 + 1) = 4 parameters instead of 6, reducing by another 2. Hence 
15 goes to 11 parameters. As 9 parameters correspond to the action of the 
global symmetry group of translations, rotations and iso-rotations, we must 
find the gradient flow curves in a two dimensional subspace, parametrized by 
a, h and R subject to the constraint a + b = R. These gradient flow curves 
would start approximately at the asymptotic critical point of the two infinitely 
separated, minimal energy, isolated Skyrmions, and arriving at the minimal 
energy toroidal configuration. Since the product ansatz tells us that asymp- 
totically this is a 10 dimensional manifold, the gradient flow must also yield 
a 10 dimensional manifold. This is intuitively reasonable, minimizing in a two 
dimensional manifold will typically yield a one dimensional "valley" or "path" 
of steepest descent linking together the critical points. Hence we do indeed ob- 
tain a 9 -|- 1 = 10 dimensional manifold of most attractive channel instanton 
generated Skyrme fields. 



The gradient flow has not to date been calculated. Hosaka et al 55 1 have 



exhibited a qualitatively similar manifold of constrained minima. It is obtained 
simply by letting the set of triangles vary from that of the right isoceles triangle 
corresponding to widely separated Skyrmions, to the equilateral triangle of 
the toroid through intermediate symmetric isoceles triangles. The Hartshorne 
ellipse starts at very high eccentricity and varies until it degenerates to the 
circle. The energy is minimized for each intermediate triangle, fixing the value 
of R. The minimal energy for fixed eccentricity decreases monotonically until 
the circle is reached at a/b = 1. The constrained energies are always within 
1% or 2% of the similar but fully numerical computations of Verbaarschot et 



al 56 and of Walhout 57 



Geometrically the 10 dimensional most attractive channel manifold consists of 
the direct product of a 6 dimensional (global) manifold, which is generated by 

3 independent spatial translations and 3 independent isospin rotations, with a 

4 dimensional (relative) manifold. The 4 dimensional relative manifold consists 
of "centered" Skyrme fields, which can be acted upon by the group of spatial 
rotations. The manifold is parametrized by a coordinate depicting separation 
and 3 angular coordinates taken without loss of generality to be the Euler 
angles specifying a frame of unoriented Cartesian axes. The generic SO (3) 
orbit is actually only S0(3)/V^, where V is the group of 180° rotations about 



53 



the 3 axes (and the identity), which of course leave an unoriented Cartesian 
frame invariant. At minimal separation the orbit degenerates to IRP2 which is 
the same as the sphere S"^ with antipodal points identified, physically it is the 
orbit of the symmetry axis of the toroid, taking into account that this axis is 
unoriented. 

The Atiyah-Hitchin manifoldj^, corresponding to the moduli space of cen- 



tered 2 BPS-monopole configurations has exactly the same orbit structure. 
Hence the Atiyah-Hitchin manifold is a good candidate describing centered, 
attractive channel Skyrmions. The Atiyah-Hitchin manifold has an implicit 
metric, the Atiyah-Hitchin metric. This metric, however, is not appropriate 
for Skyrmions. The Atiyah-Hitchin encodes in it velocity-dependent Coulomb 
interactions between monopoles, which are absent for Skyrmions. The true 
metric must be calculated using the Skyrme energy functional. This has been 



done by Leese et al |59 



3.3 Rational maps and multi-haryon number states of the Skyrme model 



The study of high baryon number solutions of the Skyrme model has always 
been a very difficult problem to tackle. What is lacking is a good ansatz which 
captures the symmetries and simplifies the equations of motion. 



A first glance at higher baryon number was given by Braaten et al |^D| almost 
10 years ago. Using relaxation methods on a highly (in those days) powerful 
Cray super-computer, they isolated states for the sectors 5 = 2 to 5 = 6 and 
computed their energies. Quite surprisingly, the configurations for i? = 3, 4, 5 
and 6 took very geometrical shapes, tetrahedral (5 = 3), octahedral (5 = 4), 
and less symmetrical ones for B = 5 and 6. Obviously, these look nothing like 
nuclei as described by shell or droplet models from traditional nuclear physics. 
Of course, the doughnut shaped B = 2 state, the Skyrme model's deuteron 
state, looks very little like a pair of weakly interacting nucleons. However, as 
the classical binding energy of the doughnut (100 MeV) is much greater than 
the deuteron's real binding energy (around 2 MeV), one could hope that the 
B = 2 case is not typical for the Skyrme model and that things would settle 
down for higher baryon number. We will see that it is not the case, and one 
gets the impression from the literature that the results for B = 2 to 6 were 
both unexpected and not understood: consequently they were left pretty much 
alone for the interim. In 1996 Battye and Sutcliffe confirmed these results 
(except the B = 6 state which seems to have been misidentified) using state- 
of-the-art software and hardware, as well as found the structure of probably 
minimum energy states with B = 7,8 and 9 (see figure |T3| taken from reference 



61|). This indicates that indeed the Skyrme model favours regular geometric 



configurations instead of shell model- type lumps. These shapes in fact, are 
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Fig. 13. Various baryonic density for Skyrme states of baryon number 2 to 12. 

very similar to those found in chemistry. Only time and further investigations 
will tell if indeed and why such configurations are physically relevant. 

One must note here that all these states have been found using numerical 
relaxation methods, with the algorithm feeling its way down the slopes of con- 
figuration space in search of the lowest energy possible in each given baryon 
number sector. Of course, it is impossible to be absolutely sure that the con- 
figurations obtained so far are absolute minima of the energy in each baryon 
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number sector, and not merely local minima, without using other methods. 
Thus there is still a possibility that the spectrum of real states might be dif- 
ferent, even though the fact that the energy per baryon number is very close 
to the Bogomolny bound (about 10% or so over the bound) makes this pos- 
sibility quite remote. Oscillations and perturbations around those solutions 
have been considerd to test (at least locally) the stability of the solutions, and 
to compare the resulting spectrum with known excited states of nuclei. So far 
the results look encouraging but there is still a lot of work to be done before 
comparison with nature can be made in a serious manner . 



As impressive as these numerical results might look, they are still a long way 
ahead of actually doing simulations of Skyrmion scattering processes for sev- 
eral initial Skyrmions or for large energy. Even though it actually is possible 
to do so, it still takes strong numerical skills to do the simulation and extract 
and interpret the results. To be able to push forward the current study of the 
model, and also to avoid that the study of soliton scattering is absorbed into 
purely numerical or computational physics, one needs some kind of analytical 
handle on the problem, even an approximate one (for a flavour of the diffi- 
culties encountered in the study of soliton-sohton scattering see the article by 
Crutchfield and Bell||6^). The instanton method which we have just described 
does exactly this but it is of little help for the study of large baryon number 
solutions, except under very special conditions |3^. The gradient flow curve 
method which we will present in section 4 does not help much either since 
to apply it we need to know the manifold connecting the critical points of 
the system for low energy in a given baryon number sector. However we do 
not know this manifold, and we actually deduce its structure from numerical 
simulations. 



Help comes our way from the study of the BPS monopole system. It is not 
our intention to discuss this system here (see the review by Sutcliffe[ |65|] ). 
Pioneers of this model are Bogomolgny||6^, Prasad and Sommer field ||6 7|| , 't 
Hooftg, PolyakovH], Atiyah and Hitchin[^,0, Manton^, Gibbons|3§, 
Ward||73| and Nahm[[f^, to mention a few. The BPS monopole is a topolog- 
ical soliton of a massless Higgs type model with an SU(2) gauge symmetry, 
which saturates the Bogomolny bound of the theory [^,^. Over the years, 
the spectrum of the model has been studied, and states with each value of the 
magnetic charge (which represents the winding number associated to the soli- 
ton in this model) have been isolated. Quite surprisingly, there exists a state 
with magnetic charge N with symmetries identical to those of the Skryme 
model for baryon number A^[0] (this fact has been verified for ranging 
from 1 to 9): toroidal for N = 2, tetrahedral for A^ = 3, etc. This of course 
does not indicate that BPS fields should be used to study baryonic systems, 
but it does indicate that the mathematical arsenal used in the study of the 
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BPS system could perhaps be generalized to Skyrmion system^. The tools 
of interest here are called rational maps and were introduced in the study of 
BPS systems by Donaldson[^, and recently elaborated upon by Jarvis|]78[ . 



We will not explain in detail the method nor show how it works. Instead we 
will give a taste of how it can be applied to Skyrmion systems following the 
article of Houghton et al p4 |. 



BPS monopoles are solutions of a model possessing a symmetry breaking term 
which breaks an SU(2) symmetry via a triplet of Higgs fields to a U(l) sym- 
metry, interpreted as electromagnetism. Inside the monopole, the SU(2) sym- 
metry is intact, but it has to be broken to U(l) on the outside so as to give the 
soliton finite energy (because of a potential term in the Lagrangian density). 
Restricting the Higgs field to its broken symmetry value outside the soliton 
fixes its length, but not its direction: it can take any direction in IR^ and in 
fact describes an S"^. Infinity in flat 3-dimensional space also consists of a 2- 
sphere, so the Higgs field of the monopole is a map from S*^ to S"^. Such maps 
are divided in disjoint homotopy classes numbered by the number of times the 
first is wrapped around the second S*^. This topological winding number is 
proportional to the magnetic charge of the field. The vacuum has the Higgs 
field pointing the same way everywhere (thereby having winding number and 
magnetic charge 0), while the unit magnetic charge monopole looks like a 
hedgehog from afar. By stereographic projection one can transform the map 
from S"^ to S'^ to another which maps the complex plane in another complex 
plane (if we identify a given point of the spheres with infinity in the usual 
way). Donaldson|7^ showed that there is a one to one relation between the 



field of an monopole and rational maps R{Z) of degree N. A rational map 
of degree fromC — >(D is defined as: 

Z ^ R{Z) = ^ (221) 



where p{Z) and q{Z) are polynomials of at most degree A^, with at least 
one being of degree A^ and with no common roots. The parameters of the 
polynomials generate (much like in the case of the parameters of the instantons 
of the previous subsection) a finite dimensional manifold of configurations and 
can be chosen so as to give the soliton some desired symmetries, and fine tuned 
to lower the energy as much as possible. We refer the reader to the literature for 
further details |3^. Experience shows that once a particular set of symmetries 



^ In fact, tlie similarity between BPS and Skyrme systems could run deeper than 
it appears. During his talk at the CRM-Fields-CAP 1997 workshop "SoUtons" in 
Kingston, Canada, N.S. Manton conjectured that the moduli space of vortices is a 
submanifold of the moduli space of BPS monopoles, which is itself a submanifold 
of that of Skyrmions, itself being included in the moduli space of SU(2) instantons. 
This fascinating conjecture remains to be established. 
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Fig. 14. Stereographic projection of an 5*^ of radius 1/2 on plane {XY). It gives 
Z = X + iY = tan {9/2)e^'^, and the projection angle a = 9/2. 

has been chosen and implemented in the rational map, the method gives a 
good approximation to the exact (numerical) result. 

Houghton et al |34] generalized this method to the Skyrme model by choosing 
the following ansatz for Skyrme fields: 

U{x) = e^^('')^R(2)-" (222) 



where ^^(z) is a function of Z, or in other words of the usual spherical angles 
6 and only, /(r) depends solely on the distance r from the origin. This 
ansatz is interesting in that it singles out the distance to the origin from the 
angular coordinates. As Houghton et al noted, one can understand this ansatz 
as mapping the two-spheres centered on the origin of space onto the 2-spheres 
which correspond to latitudes in the (of SU(2)). To do this the two-spheres 
of the IR^ of space are mapped via stereographic projection onto the complex 
plane (with complex infinity identified to a single point of the sphere) as shown 
in figure |T^. Elements Z of this plane are then functions of the spherical angles 
6 and 0: 

Z = tan (^/2)e''^. (223) 



We are now on familiar ground: elements of the initial complex plane Z are 
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then mapped by the rational map R{Z) into another complex plane, itself 
obtained by stereographic projection of the latitudes of the SU(2) three-sphere: 

"x" 51 — > 53 

{U{x)} 



1R3 


R X 52 
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{^R(Z)} 
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{R{z)} 



fiRiz) is expressed as follows: 



^R(z) = M Mp|2| (2^(-^)' 2S(^)' 1 - (225) 



11 + 



(224) 



where ^{R{Z)) and Q'(-R(Z)) represent the real and imaginary parts of R{Z) 
respectively. nR(z) depends on the rational map R{Z) chosen, and gives, for 
example, the ordinary radial vector n if R{Z) = Z. This particular map re- 
produces of course the ordinary Skyrmion. 

The parametrization ( P22| ) chosen for U can seem a bit restrictive. It is espe- 
cially suited to reproduce configurations localized around the origin. It prob- 
ably does not reproduce accurately, well separated Skyrmions (although it 
works surprisingly well and gives an idea of some processes). The main advan- 
tage of this parametrization is that it decouples contributions from the radial 
{i.e. ones related to /(r)) and from the angular part {i.e. related to R{Z)). 
This way, it is possible to first impose a given symmetry to the configuration 
and then minimize the angular contribution to the energy. Subsequently by 
minimizing with respect to /(r) we obtain the minimum energy possible within 
the ansatz. We note that in certain cases, varying the value of some parame- 
ters of R{Z) gradually and minimizing with respect to /(r) gives "snapshots" 
of certain scattering processes. The case of -B = 3 is a good example. 

Choosing the parameters in the polynomials p{Z) and q{Z) so as to get solu- 
tions with the right symmetries is a technical but important point which we 
will illustrate by the following example. Let us compute the expression for the 
rational map R{Z) which gives the Skyrmion field U{x) = Uq + i u ■ t with 
lowest energy in the B = 2 sector. As we will see in section 4.3 the following 
three symmetries are characteristic of a pair of Skyrmions converging together 
at the origin to form a tightly bound minimum energy toroidal solution of the 
equations of motion: 
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The most general expression for the rational map R of degree 2 is given by: 



R{Z) 



AZ2 + ^Z + u 



(229) 



where a, /3, 7, A, and v are constants (real or imaginary) to be determined 
by imposing the symmetries Ji, I2 and J3. To do this, we need to parametrize 
them using Z and R{Z). Using (p23|) and 



u = sin/(r)nR(z) 



(230) 



we see that the transformation x — ^ —x of Ji is equivalent to Z ^ —Z after 
projection on the complex plane. Similarly, y —y and z —>■ —z just translate 
to Z ^ Z and Z l/Z, respectively. For our needs we will only consider 
the new transformations Xi = /1/2 and X2 = I2I3 (the product J3/1 being 
discarded since it brings no new constraint on R). The symmetry Xi is given 
then by: 



Xi : Z 



n 



n 



n 
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n 



(231) 



The field nR(z) being then invariant under the transformation, this imposes 
the constraint 



R{-Z) = R{Z) 



(232) 



on R. Similarly, I2 is written as: 



Xo : Z 



l/Z: 



n 



n 



n 



(233) 
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A short calculation shows that n transforms under this transformation like 
'^R(z) '^i/R(z)- So R also has to satisfy 



R{l/Z) = l/R{Z). 



(234) 



We now apply the conditions (232) and (234) to the general form of i?, (|229|) 
to fix the constants a, /3, 7, A, /x and z/. Equation ( |232D sets /? and /x to 0, 
while (|234|) gives the following constraints on the remaining parameters: 



a7 = u\ 



and 



2 I 2 \ 2 I 2 



(235) 



We can set z/ = 1 by scaling each variable then we find: 



R{Z) 



or 



7 



(236) 



-7^2 + 1 7^2 + 1 

The second possibility corresponds to a rotation by 90°, Z \Z . One then 



replaces this expression in the energy of the ansatz (222) given by 
E = 47r 1 {r^f + 2iV(/'2 + 1) sin^ / + j dr 



(237) 



where 
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(238) 



As seen earlier, the contribution to the energy coming from the angular part 
of the ansatz decouples from the radial part, and X can be minimized as a 
function of the parameter 7. One finds the extremal value X = vr + 8/3 for 
7 = 0. This simply corresponds to the map RiZ^ = Z^ which possesses the 
cylindrical symmetry 

i?(e'") = e2*"i?(Z), (239) 



which contains the discrete symmetries Ji and l2- This is the toroidal config- 
uration which we will discuss more when we describe the low energy manifold 
for the sector B = 2. Replacing the minimum value of X in the expression of 
the energy and minimizing further relative to the radial function /, one finally 
obtains a value only 3% greater than the mass of the torus obtained by a fully 
numerical computations. Rational maps with more complicated symmetries 
generalise the previous discussion and are presented in reference . 
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We will close this section with the following comments about the number of 
"holes" (regions with zero baryon number density) which are present in a given 
configuration. dR/ dz is zero if the Wronskian is zero 

W{z) = p'{Z)q{Z) - p{Z)q'{Z) = (240) 



which is the numerator of dR/ dZ. The baryon number for the ansatz (|222|) 
is given by 
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(241) 



The baryon number density is then proportional to dR/ dZ and vanishes where 
the Wronskian is zero. This means that it will be zero along rays pointing to 
infinity from the origin, and whose directions are given by the roots of W{Z). 
R{Z) being of degree N, W{Z) is generically of degree 2N — 2 (the naively 
leading power of Z cancels in ( p40|) ), and so a configuration of this baryon 
number should have 2N — 2 holes in it. For the case N = 2, q{Z) = 1 so the 
Wronskian is in fact linear in Z and only has one zero, which is consistent 
with the trivial fact that a torus only has one hole. 



3.4 Search for a sphaleron in the Skyrme model: Morse theory 



Morse theory [[79| , p0| , pl| ,p^ , otherwise known as global variational analysis, re- 
lates the topology of a manifold to the number of and types of critical points 
of a function defined on the manifold. For application in field theory, the 
manifold in question is the (infinite dimensional) space of all field configura- 
tions and the "function" defined on this manifold is generally a functional, 
typically the energy functional or perhaps the action functional. The general- 
ization of Morse theory to the infinite dimensional arena goes under the name 
of Ljusternik-Snirelman theory [p^ . 



The classic, illustrative example of the application of Morse theory is furnished 
by a function defined on a torus. We take the outer radius to be R and the 
inner radius to be r < i?, with the symmetry axis pointing along the x- 
axis. The function in question should be a "Morse function" , the definition of 
which we will address below. We will consider the function defined by the value 
of the z coordinate of the Cartesian coordinates of each point of the torus, 
which happens to be a Morse function. Physically, the function in question 
is the altitude from the ^ = plane for each point on the torus. Now let us 
consider the critical points of this function on the torus. It is well known that 
every function defined on a compact manifold achieves its global maximum 
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Fig. 15. The altitude fonction /(x, y, z) = z defined on an ordinary torus possesses 4 
critical points: the global maximum A, the global minimum D, and two minimax's 
B and C . Morse theory identifies the presence of the minimaxima with the existence 
on the torus of two types of non-contractible loops C\ and Ci- 

and its global minimum somewhere on the manifold. Clearly the function 
will be critical at these two points, hence, in fact the compactness of the 
manifold, which is is a topological characterization, has implied the existence 
of two critical points. There are however, even more. These are predicted 
by Morse theory, due to the non-simple connectedness of the torus. The torus 
admits two different non-contractible closed loops (see figure ^). Morse theory 
implies that there are at least two other critical points, which are minimax's. 



Physically we can just see them: the points B and C in figure |T5[ If we consider 
the intersection of the torus with successive planes of constant the two planes 
for which the topological nature of this intersection changes correspond to 
the positions of the minimax's. Successively, as we sweep the plane through 
the torus the intersection (see figure |1^) will commence as a point (a very 
degenerate "circle" or loop), then a normal loop, but then at one point the 
loop will pinch in and touch itself and then break up into two disjoint loops. 
This is the first minimax or saddle point. Then the two loops will separate 
a little but again come back and touch. When they touch is another saddle 
point, and further they will separate to form a single loop, which will finally 
terminate by degenerating to a point. 

To prove the existence of an additional critical point due to non-trivial topol- 
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Fig. 16. Intersection of the torus with planes {a.) z = R + r, (b) z = R, (c) z = R — r, 
(d) z = 0, (e) z = -R + r, (f) z = -R and (g) z = -R-r. 

ogy, consider the somewhat more general example of a smooth function f{x) 
defined on a compact manifold with a non-contractible loop. We wish to prove 
that there exist at least three critical points for the function. It is immedi- 
ate that the global minimum /(xmin) and the global maximum /(xmax) exist. 
Without loss of generality /(a;max) > /(a^min), since if /(xmax) = /(a^min), the 
function is a constant. Also without loss of generality we may assume that the 
maximum and minimum are achieved as unique, individual, isolated points. 
This means that there are already two critical points hence if the minimum or 
maximum is achieved elsewhere there would be at least three critical points 
and we would be done. 

Now consider the set of non-contractible loops in the manifold which go 
through the minimum. For each loop find the point at which the function 
is maximum, the set {x}. If this occurs at several points for any one loop, 
choose the point x where |V/p is minimal (since we are looking for a critical 
point where |V/| = 0). The maximum f{x) is necessarily greater than f{xmm) 
since if it were not then f{x) would have to be critical at x, and then would 
already admit a third critical point at x. Now we find the minimum of the set 
{/(x)} by varying the loop. The point where this occurs must exist and will 
correspond to a critical point of the function, a saddle point. It must exist 
because a bounded, monotone, sequence in a compact manifold always admits 
a limit point. We consider the sequence of points obtained by finding x for a 
given curve and then varying successively the curve to yield xi, X2, etc. such 
that /(xi) > /(xi+i) (see figure 0). Let Xc be defined by Xc G {x} such that 
f{xc) = min(/(x)) and we call the curve passing through Xc the critical curve, 
Xc(t). V/(xc) must vanish. If |V/(xc)| > 6 {6 > 0) then consider the curve 
obtained by deforming the points of the critical curve in a neighborhood of Xc 
along the direction opposite to the gradient: 



Xc(t) ^ Xd(t) = Xe(t) - e(t) V/(fe(t)), 



(242) 
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Fig. 17. Drawing of the successive paths passing through Xmin and the local maxima 
X[ for i = 1, 2, 3, ... of the curves which converge to the real minimax Xc- 

where e(t) is non-zero only around Xc- For e(t) small enough, the deformed 
curve is still a non-contractible curve, however the maximum of the function 
along Xd{t) is clearly less than the maximum at Xc because the gradient — V/ 
points along the direction of decreasing /. This is a contradiction since Xc is 
the minimum of {x}. Hence Vf{xc) = 0. 

Morse theory goes on to give a set of inequalities relating the number of crit- 
ical points to the changes in topology of the manifold. We refer the reader to 



the literature p5|j36|j80|j81| , |83|j8^ for more details. We note that nothing de- 
pended critically on the fact that we were considering a non-contractible loop, 
it could well have been any non-contractible compact manifold, for example 
a sphere 5^. Hence non-trivial homotopy groups Hm can imply the existence 
of non-minimal critical points. The crucial ingredient for the success of the 
minimax procedure was that the minimal critical point was non-degenerate. 
In fact, a function which only admits non-degenerate critical points is called 
a Morse function. To apply Morse theory however we do not truly require 
Morse functions. The function should simply be non-degenerate along every 
non-contractible loop. 

The application of Morse theory to infinite dimensional manifolds was ana- 
lyzed by Ljusternik-Snirelman||83|| . A very readable account of the use of this 
theory was done by Taubes for the case of magnetic monopole[^. The idea 



is to first establish the existence of non-contractible loops in the configura- 
tion space Mq, the space with net monopole number zero, i.e. Hi(Mo) 7^ 0. 
Secondly consider the configurations of an infinitely separated monopole-anti- 
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monopole pair. This configuration is in Mq, and the non-contractible loop 
corresponds to rotating the monopole relative to the anti-monopole in iso- 
space by one complete revolution. The minimal critical point corresponds to 
the situation where the monopole and anti-monopole have annihilated and all 
radiation has dissipated off to infinity leaving the quiescent (symmetry broken) 
vacuum behind. This critical point of zero energy is non-degenerate along ev- 
ery non-contractible loop. Taubes searched for the existence of a non-minimal 
critical point in the configurations bounded above in energy by the asymp- 
totic critical point of the infinitely separated monopole- anti-monopole pair and 
bounded below by the energy (zero) of the vacuum. First it was shown that 
this infinite dimensional subset of the configuration space is in fact compact. 
With this information it is sufficient to exhibit a single non-contractible loop, 
with energy everywhere less than the energy of the infinitely separated pair, 
to be able to conclude that the minimax procedure will converge. It will of 
course converge to a different critical point than the trivial vacuum, since the 
vacuum is non-degenerate. The existence of this loop was shown proving 
the existence of a non- minimal critical point. 

Bagger et al ||8^ attempted to mimic this procedure for the Skyrme model. 
The configuration space for the Skyrme model corresponds to maps from IR^ 
into the group SU(2): 

C ■ {U{x) : + cx) ^ SU(2)} (243) 

R'^ + oo = S*^, hence C separates into disjoint sectors labeled by the homotopy 
classes 113(811(2)) = 7L, the integer being of course the baryon number. Hence 

C = Co + Ci + C2 + ■ ■ ■ + C_i + C_2 + • • • • (244) 

The existence of non-contractible loops in C 

ni(C)^ (245) 

implies that 

ni(Cn)7^ 0, (246) 

and that they are actually all equal. Each sector actually contains each other 
sector, simply by constructing the requisite number of baryon-anti-baryon 
pairs, and moving, say, the anti-baryons as far away as is required. 

Hi (Co) can be seen to be exactly the same as 114(811(2)). We can see this 
through the following construction. R^ -|- 00 can be thought of as the end 
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< = 1, E,^ +00, [/ = 1 



U = 1 




i e [0,1] 



0, + oo, [7 = 1 



Fig. 18. Values of U defined on the hypercube of space. 

cap of a 4-cube. The parameter along the fourth dimension corresponds to 
the parameter along the loop (see figure |T^). If we start at t = with U = 1, 
impose that U = 1 along the vertical faces (because these faces also correspond 
to the point at infinity of IR^), and = 1 at t = 1, we define a loop in Cq. 
However a 4-cube with the surfaces identified is the same as S"^, R'^ + oo = S^, 
hence we equally well define a mapping 



Uix,t) : -> SU(2), 



(247) 



i.e. an element of 114(811(2)). It is known that||87 
n4(SU(2)) = ^2- 



(248) 



Hence there exist non-contractible loops in the space of configurations of the 
Skyrme model. These loops can most physically be realized as correspond- 
ing to a Skyrmion-anti-Skyrmion pair, as in the monopole situation, the two 
solitons are rotated relative to one another through one complete revolution. 
The difference between the Skyrme model and the monopole situation is that 
while loops involving further complete relative revolutions are distinct non- 
contractible loops for monopoles, for Skyrmions all loops with an odd number 
of complete revolutions are equivalent to each other and non-contractible while 
an even number of complete revolutions yields a contractible loop. Bagger et 



al [§6| found, at least to first order in perturbation theory for large separation 
of the two solitons, that there were no non-contractible loops where the en- 
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ergy was everywhere less than the energy of the asymptotic critical point of 
the infinitely separated Skyrmion-anti-Skyrmion pair. Augmenting the Skyrme 
model with an electromagnetic interaction, however, gave a sufficiently attrac- 
tive Skyrmion-anti-Skyrmion potential for them to conclude the existence of 
a non-minimal critical point in this somewhat modified theory. 

The problem was analyzed in greater detail by Isler et al |^, for the B = 
and B = 2 situation together. Let Us{x) be the field of a Skyrmion. Then for 
well separated Skyrmions or Skyrmion-anti-Skyrmion pair, the product ansatz 
suffices. Let 



UB=2{x)=R{t)Us{x-Xi)R^{t) Us{x-X2) 
UB=oix)=R{t)Usix- xi)R\t) Ul{x-X2) 



(249) 
(250) 



where R{t) is an SU(2) matrix that introduces a relative iso-rotation. If R{t) 
varies from any point in SU(2) to its antipodal point, i.e. 



t=i 



-R{t) 



t e [0, 1] 



(251) 



then t parametrizes a non-contractible loop in Co or C2. This is the usual 
topology which demonstrates SU(2) as the simply connected (double) cover 
of SO (3). It is a reasonably straightforward computation to find the energy to 
lowest order in the separation d = — a;2||, 



„2 2(l-cos^)(3(n-rf)2-l) 



Eb=o = 2Es- 47Tfy^ '-^ (252) 

F o p I . ,2 2 (l-cosg)(3(n-rf>-l) 

for 

R = e'^^-^/\ (254) 



Eg is the energy of a Skyrmion and k is the coefficient of the fall off 
in the Skyrmion profile function /(r), /(r) ~ /t/r^ + 0(l/r®) (see equation 
([T3|)). The potential V serves to separate the (reduced) configuration space 
consisting of relative iso-rotations into two disjoint parts. Indeed this reduced 
space of configurations is isomorphic to a 3-ball of radius 2tt modulo one 
identification. 6 plays the role of the radius, n the unit vector giving the 
direction. Furthermore antipodal points are identified, (6',n) = (27r — 6, —n). 
This identification is particularly reductive for 9 = 27r, the whole sphere at 
9 = 271 corresponds to R = —1 which is identified with the origin where R = 1. 
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3{h -df -1^0 



Fig. 19. Diagram of the attractive region for the Skyrmion-Skyrmion B = 2 poten- 
tial. 

The factor 1 — cos 9 is positive semi-definite, and equal to zero only at ^ = 
or 27r. The function 3(ri ■ c?)^ — 1 varies from 2 when n ■ = 1 to — 1 when 
n ■ (i = 0. Hence when the direction of the relative rotation is chosen such 
that h - d = ±1/a/3, this defines a double cone that passes through the origin, 
demarcating the boundary between regions of relative attraction and relative 
repulsion (see figure |l^). The regions of attraction have less energy than 2 Es, 
while the regions of repulsion have more energy than 2 Eg. 

If the attractive region is in the "time-like" region of these cones, which is 
the case for the case B = 0, there is no curve which can pass from a point 
within this region to its antipode without either crossing over to the region of 
repulsion or touching the cone at its vertex. This is evident since a point and 
its antipode find themselves in opposite sides of the "forward" or "backward" 
light cone. All paths linking them must pass through the vertex or enter into 
the repulsive region which is not desired since the energy here is greater than 
that of two infinitely separated Skyrmions (or Skyrmion-anti-Skyrmion pair). 
Hence we cannot conclude the existence of a non-minimal critical point, we 
need to find a non-contractible curve where the energy is everywhere less than 
2 E^. Even though our curves are never greater than 2 E^ (when they pass 
through the vertex) we cannot be assured that the minimax procedure will just 
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converge to the asymptotic critical point (of energy 2eskyrmion)- The situation 
just described is that which apphes to the case B = 0, i.e. a Skyrmion-anti- 
Skyrmion pair. Here, there is no question that the minimum energy configu- 
ration is non-degenerate along every non-contractible curve. Hence it remains 
an open question whether indeed there are non-minimal critical points in the 
sector B = 0. 

Considering the case B = 2 (see figure |T^, the attractive and repulsive regions 
exchange with respect to B = 0. Hence the attractive region is the "space-hke" 
region relative to the cone. It is evident that now there exist non-contractible 
loops, which remain everywhere in the "space-like" region. They simply skirt 
around the cone to the other side of the origin where lies the antipode. Hence 
we show the existence of non-contractible loops which are everywhere lower in 
energy than 2 Eg. However this is still not enough to conclude the existence of 
non-minimal critical points. Indeed in the sector B = 2 the energy functional 
fails to be a Morse function to a sufficient extent. The minimal critical point 
has been demonstrated to be a toroidal configuration with axial symmetry, in 
all but a rigorous, analytical mathematical proof. The axis of symmetry has no 
direction, the toroidal configuration rotated by 180°, about an axis orthogonal 
to the symmetry axis, is identical to the starting configuration. This implies 
that the minimal energy configuration is degenerate along a non-contractible 
loop. Hence it is the strongest possibility that the minimax procedure will 
converge simply to the minimum energy toroid. 

Even though our exercise with Morse theory has led to no new solutions it is a 
worthwhile analysis allowing us to understand the model in a more profound 
way. There are several open questions raised by the analysis, two evident 
ones are do there exist sphalerons in the B = sector and do the existence 
of non-trivial higher homotopy groups imply existence of sphalerons. Indeed 
Morse theory has been used in relation with rational maps in a more recent 
article g. 



4 Dynamics of the Skyrme model: Soliton-Soliton scattering 

The dynamics of solitons is an extremely interesting and complicated problem. 
There are many different modes of excitation for a single soliton itself. There 
are certain modes, the zero modes of the classical small oscillation problem 
about the soliton, which properly belong to the soliton itself. They are usually 
treated semi-classically via Bohr-Sommerfeld type quantization rules. There 
are also regular vibrational and resonant modes which correspond to excited 
states of the soliton. Additionally there are modes which correspond to the 
scattering of (non-solitonic) waves off the soliton itself, for example, pions 
scattering from nucleons. 
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Next one can consider the interaction of two solitons with each other. Here 
one can support many forms of reactions: scattering, deformation, bound 
states, annihilation among others. Sohton-anti-sohton annihilation is partic- 
ularly problematic because of the coexistence of both perturbative and non- 
perturbative regimes. A physical example is evidently nucleon-anti-nucleon 
annihilation. The potential between a nucleon and an anti-nucleon has been 
obtained by Lu and Amado[]T3[ using the product ansatz for large distances, 
and by Lu et al |jT^ for distances larger than 0.8 fm using numerical methods. 
Physically there emerges a critical distance which is of the order of 0.8- 
1.2 fm, outside of which the interaction between the particles is essentially 
repulsive. On the other hand if they attain this critical distance, they will 
quickly combine together into a lump of mesonic matter of baryon number 
zero. Numerical simulations of the classical system are essential to understand 
this complicated process [|T^]. It was shown that the reaction for transforming 
the pair of particles into a single lump of zero baryon number happens at the 
limit of causality and that the energy left by the disintegration remains local- 
ized for a relatively long time (until pion radiation waves disperses it into the 
vacuum) |]TB[ . The remaining part of the process, the emission of pion waves 
from the lump of mesonic matter, has been studied in the most detail, using 
path integral methodsfl^, or coherent state methods [p^ . The results repro- 
duce well the experimental phase shifts. We will not treat soliton-anti-soliton 
annihilation any further. 



Even more complicated situations arise as we increase the soliton number. 
We can go to the point where one has an infinite number of solitons and all 
their various phases, fluid or solid, with possible crystalline structures. In this 
review we shall only discuss the low energy interaction of solitons, specifically 
Skyrmions, corresponding to the low energy scattering of nucleons. 



Scattering of nucleons within the true microscopic theory, the standard model, 
is impossible to treat satisfactorily, essentially because of our inability to com- 
pute anything in the low energy domain. Even for very high energy scattering 
the final processes leading to hadronization are not computable from the mi- 
croscopic theory. For low energy processes we have effective field theories, such 
as the Skyrme model, which afford more tractable descriptions of the physics 
involved. But even here, the baryons are represented as solitons and an exact 
quantum description of soliton states is still lacking. The only perturbative 
expansion feasible seems to be the semi-classical approximation. 



The semi-classical approximation serves well to describe constituent properties 
of individual solitons. Essentially, the procedure is to identify the low energy, 
collective modes of the soliton and to quantize them. The interactions between 
solitons can only be treated perturbatively at large separations, by computing 
the effective interaction potential between them. 
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In the past ten years there have been substantial advances in treating the part 
of the interaction at short distances and how to describe the scattering within 
the semi-classical approximation. Even the semi-classical approximation is not 
exactly solvable: soliton-soliton scattering involves an infinite number of de- 
grees of freedom except for some very special cases of integrable models. The 
problem could be tractable if there were some way of truncating to a finite 
number of degrees of freedom. Exactly such a truncation was suggested by 
Manton|3. 



In very general terms, one is interested in finding the low energy degrees of 
freedom of the two-soliton system. Typically one finds only a finite number 
of relevant degrees of freedom. The low-energy motion can then restrict itself 
self-consistently to these degrees of freedom. 

The canonical example of such a situation was provided by the case of magnetic 
monopoles in the so-called BPS limit p8|j58[| . Here the inter-monopole force 
vanishes exactly, the magnetic Coulomb repulsion being exactly cancelled by 
an attractive force due to the existence of a massless scalar exchange. Hence 
there exist static solutions with two monopoles situated at arbitrary relative 
orientations and positions. This set of configurations corresponds to a sub- 
manifold of the set of all configurations and is called the moduli space. Indeed 
the characterization in terms of positions and orientations makes sense for 
monopoles when they are well separated, but as they come close together they 
lose their identity. What is preserved is the dimension of the space of moduli. 
For large separation, the moduli describe the position and orientation of each 
monopole: there are three degrees of freedom for the position of each monopole 
and one internal phase (related to the residual U(l) gauge symmetry) giving 
a total four degrees of freedom per monopole. As the two monopoles approach 
one another, only the dimension of the relative moduli space remains, the 
monopoles deform completely and fuse into a single entity. Clearly, since there 
are no forces, the moduli describe an equipotential surface. This surface is also 
the set of minimal energy configurations in the two monopole sector. 

We can make an intuitive analogy with the surface of the earth. If this surface 
were perfectly spherical and frictionless, the potential as a function of radius 
would be (effectively) infinite at the radius of the earth and the equal to the 
usual gravitational potential for larger radii (figure ^). The radial motion, 
although not simple harmonic motion, corresponds to highly energetic modes. 
If we start at some point on this idealized earth, with an initial velocity that 
is tangent to the surface and arbitrarily small in amplitude, it is intuitively 
clear that the motion will remain very close to the surface of the earth. It is 
easy to prove in this case that the motion will follow geodesic curves on the 
surface of the earth (we will clarify somewhat the notion of geodesic later on). 

The problem of the monopoles is analogous. If the initial conditions corre- 
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Fig. 20. Potential V{r) on a hard, frictionless sphere as a function of the radius r. 
Re ~ 0.7 is the radius of the sphere. 



spond to being on a point of the moduli space of two monopole solutions, 
and the initial velocity is tangent to this surface and arbitrarily small, it 
was proposed [pTj^ ^ that the subsequent movement follows the appropriate 



geodesic in the manifold of the moduli space. This produces a concrete exam- 
ple of the truncation of the infinite number of original degrees of freedom to 
a finite number of relevant, effective degrees of freedom. 

The general situation with solitons cannot be adequately described in terms 
of a moduli space as in the case of BPS monopoles. For example, Skyrmions 
or non-BPS limit (non-zero Higgs mass) monopoles experience inter-solitonic 
forces; hence no static solutions exist corresponding to arbitrary relative po- 
sition of the solitons. Typically there do exist several low-lying critical points 
of the energy functional which should be involved and are important in the 
low-energy dynamics. In models where the solitons are not confined, one such 
critical point corresponds to infinitely-separated solitons. In addition, if con- 
figurations of energy lower than twice the energy of one soliton exist, then the 
minimum energy configuration represents another critical point corresponding 
to a bound state of two solitons. Furthermore there could be other metastable 



solutions such as the dibaryon of the Skyrme model[^, or the sphaleron so- 
lutions of the Weinberg-Salam model p6|,|4|, among others. The low-energy 
dynamics will restrict itself to these critical points and certain paths linking 
them together, as we will see in the next subsection. 
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4-1 General formalism 



Manton|^ suggested a possible truncation of the degrees of freedom to de- 



scribe the low energy dynamics of the soliton in these more general situations. 
He suggested that the dynamics would truncate self-consistently to the union 
of all the low energy critical points and a set of curves which pass between 
these various critical points. These curves are alternatively taken to be the 
paths of steepest descent or the gradient flow curves linking the critical points 
together. The gradient flow method or steepest descent method gives rise to 
only slightly different sets of configurations if the critical points are truly low 
lying and the gradients are small, as we will now discuss. There is also another 
formalism called the valley method which also serves to give the paths which 



connect the critical points, however, we shall not discuss this method here|]89[|. 



Gradient flow curves are mathematically described as the integral curves of the 
vector field corresponding to the gradient. This gives a first order differential 
equation: 

dxJ(A)_ dViAX)) 

^''^r- dx^ ^^^^ 



where x\X) are the coordinates of a point along the curve, g^-^ is the metric 
on the space of all configurations and V^(x') is the potential defined on it. 
The initial directions taken corresponding to the unstable directions (negative 
modes) are extracted from the matrix of second derivatives of the potential at 
the position of the critical points. 

The steepest descent curves, however, are defined in Manton[3^] by the equa- 
tion 

d^xj dV , , 

= -a? <^^^) 



augmented with the boundary condition that 

lim x\\) = x'q (257) 

A— »— oo 



where the coordinates of the unstable critical point. These correspond, 

actually, to the approximate dynamical trajectories followed by the system 
starting at t = — oo at the critical point and moving in the unstable direc- 
tions while neglecting terms that are quadratic in the velocities, t corresponds 
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exactly to A. The exact dynamics would follow from the equation 



which is written in differential geometry as the geodesic equation 

9^,{x^+rlxH'^^-^ (259) 



where the Levi-Civita connection is 

^^'-2^ + dx^ dx-J- ^^^"^ 



It is evident that the approximation is valid when x^ are very small. 

If we consider the same initial conditions with the dynamics augmented by a 
damping term representing friction, the steepest descent curves will naturally 
go over to the gradient flow curves 

d^x^ dx^ dx^ dV , dx^ /^^-.n 

^« dA^ + ^ijrM^^ = -Q-- bg,^ (261) 



(6 3> measures the amount of damping in the system). Reparametrizing the 
curves with A — 6A and taking b — > +oo wc get 

dV drJ 



which is just the gradient flow equation. 

We believe that all three equations yield essentially the same manifold to 
which the dynamics should truncate for a large variety of dynamics. What 
is required is that variations between the trajectories which determine the 
sub-manifolds go to zero. This will occur if the potential in the "transverse" 
directions is very steep. All three trajectories will be pushed together because 
of energy considerations: the trajectories cannot vary too far from the gradient 
flow trajectories because it costs too much in energy. In 4.2.3 we will discuss 
a specific example wherein the three methods will be compared. 

We will show with several examples how the truncation of the dynamics comes 
about and also some examples of how it can fail. 
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4-2 Examples of truncation for a few systems 



4-2.1 Particle subject to a potential with symmetry 

First we treat the simplest, non-trivial example conceivable, which we have 
already mentioned: a particle on a spherical surface. Consider a Lagrangian 
in with Cartesian coordinates and Pythagorian metric 



where = 6ijX^x\ 9 and cf) are the usual spherical angular coordinates, A and 
are real and positive constants. For A = the space of static solutions is 
just IR'^ and the general motion is along straight lines. For any non-zero A 
the set of static solutions is the set of sub-manifolds of IR^ corresponding to 
critical points of the potential. For this potential the minimum corresponds 
to the manifold of the sphere of radius a. In addition there is an unstable 
critical point at the origin. The typical motion for any value of A depends on 
the total energy. The motion is of course bounded since the potential energy 
rises without bound for large r. The effective potential for the radial motion 
is 



where / is the (conserved) angular momentum. This potential has only one 
minimum, pushed further out from r = a (its position for / = 0) because of 
the angular momentum barrier. If the total energy of the system is fixed while 
A — >■ -l-cxD then the motion is energetically bounded to stay within a region 
that is arbitrarily close to the minimum of V{r). This minimum corresponds 
to the surface r = a. 

Indeed, the Lagrangian ( |263| ) after the substitution r = a + t/\/\ where e = 




(263) 




(264) 



e{t) and a S> e/ V^, b 



ecomes 




A2 Ax/A 



4ae3 
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" '"2 I „• 2 /) j2\ I t /I „2 ,2 



+ sin^ ^0^) + -- - - 4 (265) 
2 2 A 

after applying the condition a ^ t/\f\. We see that the radial motion and the 
angular motion completely decouple as A — +oo. To leading order the energy 
in the radial oscillations is 

E, = ]^ 8a^ el (266) 



where eo is the amplitude of the solution e{t) = eo sin (2v^aA t) to the equation 
of motion of e, and which should be taken to be independent of A. r then 
oscillates between a — eo/VX and a + eo/y/X, while the energy in the angular 
motion is 

Ee,^ = ^a\e^ + sm^e^^). (267) 



The movement of the system with arbitrarily small initial velocity tangent 
to the sphere will be governed exclusively by the kinetic term of the La- 
grangian. The dynamical constraint of being forced to live on the sphere of 
radius a will be expressed by the appearance of the metric induced on the 
spherical sub-manifold of IR^ by the Pythagorian metric already existing in 
the ambient space. This is obtained by choosing appropriate coordinates on 
the sub-manifold, = 6, v? = 0, the standard spherical coordinates, and to 
replace for x' = x^{u^) in the metric of the ambient space 



ds =5\\ dx^ dx-' 



dx' (9xj . u . 1 
d^^^a^sin^^d^^. 



(268) 



So the effective Lagrangian will then be 



1 du" (kj) 



2 ' dt dt 

= i(a'^' + a'sin2^02). (269) 

The orbits of such a Lagrangian are great circles around the sphere, which are 
in fact exactly the geodesies of the metric (p6^ , and are traced out by the 



motion at essentially constant velocity. We compare in figure ^ the "exact" 
numerical solution with the approximate solution r = a and find a good 
agreement even for several revolutions around the sphere. 
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Fig. 21. Comparison of the "exact" (left) and approximate solution given by the 
path of steepest descent (right) for a particle subject to potential ( p63| ) . They are 
both very close to the 2-sphere 5^ (a = 1, A = 3, vq = a, vq = and = 0-5 here). 

4-2.2 Particle on logarithmically deformed two-dimensional space 

A second instructive example is analogous to a scattering problem. This ex- 
ample shows how simple changes in the metric can induce radical changes in 
the motion. 

Consider again the motion of a particle in with Pythagorian metric and a 
potential corresponding to a Lagrangian 

= ^[x^ + f + z^]-Xg\x,y,z). (270) 

We consider for the present purposes a function g{x^) which admits solutions 
to the equation 

g{x') = 0. (271) 

Then as A ^ +00 the low energy motion will be restricted to the sub-manifold 
described by ( p71|) . Take the example 

(y'(x') = 2; — Inr = 0, where r = \J x"^ + y'^. (272) 

The sub-manifold is the surface of revolution obtained by turning the function 
z = ln(x) about the z axis (essentially it is a horn with a singular vertex at 
X = y = and z = —00). A good set of coordinates to use are the usual 
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polar coordinates in the (xy) plane. As A — > +00 we can neglect the motion 
transverse to the surface, in which case 2; = Inr and the potential becomes 
zero while the kinetic term gives 



2 I ^2 I „2n2 



_ 1 

~ 2 



1 ( -2 



(273) 



Hence the dynamics corresponds to geodetic motion on the 2 dimensional 
manifold with metric 



d5^ 



l + i_)dr2 + r2d^2 



(274) 



It is clear that as r — > +00 the manifold becomes flat. The equations for the 
motion are 



(275) 
(276) 



This is just a "central force problem" with an r-dependent mass. The second 
equation corresponds to the conservation of angular momentum, /, 



(277) 



Conservation of energy yields a second integral of motion 
1 



(278) 



The orbit equation is 



dr _|_ 

d^~ i^T^n 



V2er2 - Z2 



(279) 
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Fig. 22. Plot of the radial effective potential U{r) for a particle restricted to move 
on the logarithmic surface z = lnr (vq = 1 and b = 0.5). 

where e = Vq/2, I = bvo, vq is the initial velocity and b is actually the distance 
of closest approach. We see that r >b. We can define the effective potential |9( 

f/(r) = V{r) + f/2r'^ by 



-P{e-U{r)) 



P. 



(280) 



This yields 

f/(r) = 



1 / /2 

r2 + 1 V 2 



vl l + b^ 
2 1 + r2 



(281) 



which gives the graph in figure 



There is no bound orbit; the turning point of the scattering occurs for r = b 
when U{b) = Vq/2, i.e. there is no more kinetic energy in the radial motion. 
The orbit equation is integral in terms of elliptic integrals and yields 



e{r) 



E{ib, 1/r) 



(282) 



where E is the elliptic integral of the second kind. This gives a mild attractive 
potential. Typical geodesies bend slightly towards the singularity at r = for 
large impact parameter b, however, for b near zero the particle makes several 
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Fig. 23. Trajectories for the particle trapped on the logarithmic surface for different 
values of the impact parameter b: (a) b = 1.5, (b) b = 0.5, (c) b = 0.35 and (d) 
b = 0.135 (The particle is coming from the left and its velocity is irrelevant). 



revolutions about the singularity before returning to infinity (see figure ^31) . 
The 6 = geodesic is of course singular and the particle falls into the hole at 
r = 0. 



4-2.3 Systems without a continuous set of static solutions 



Returning to the first example of the sphere, we augment it with a potential 
in the z direction. Physically this could correspond to a uniform gravitational 
field. 
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The Lagrangian that we consider is 

L = ^{r^ + r^{9^ + sm'^ 9(p'^)) - \{r^ - a^)^ ~ Ar cos9. (283) 

with A > and A ^ A. Now the critical points of the potential correspond 
only to the absolute minimum located near the south pole of the sphere {r = a 
and 6 = 7t), and a minimax near the north pole at r = a and ^ = 0. The 
conditions for the critical points are 

4A(r^ - a^)r - A cos 6^ = 
Ar sin ^ = 

which give 6 = 0, n and 

A /An 
^ = "±8A^ + ^H- 

This case falls directly into the category for which Manton proposes that the 
motion truncates to the unstable manifold formed from the union of the curves 
of steepest descent or the gradient flow curves connecting critical points. The 
proposal corresponds to the following equations: 

the gradient flow curve (|255| ): 



r = -4Ar (r^ - a^) - A cos 9 (287) 

r2^ = Arsin^ (288) 

sin^ 9(p = (289) 

or the paths of steepest descent ( |256| ) 



f = -4Ar (r^ - a^) - A cos 9 (290) 
r26'' = Arsin^ (291) 
sin^ 94> = (292) 



or the complete equations of motion ( ^5^ ) 

f = r{9'^ + sin^ 9 cj?) - 4Ar(r^ - a^) - A cos 9 
2rr + r'^9 = Ar sin 9 + sin 9 cos 

-^[r^ sin ^01 = 0. 
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(284) 
(285) 

(286) 



(293) 
(294) 

(295) 



As A/A ^ +00, all three methods simply reduce to the sphere r = a, and the 
motion is constrained to this sphere with an effective potential 

V;ff(^) = aA cos^. (296) 



The difference in the three method gives rise to surfaces which vary by a 
thickness of the order of A/ A. This can be important if this ratio is not small. 

Numerical studies show that the motion is well approximated by motion on the 
truncated sub- manifold, with an important caveat: if the particle approaches 
the saddle point too closely the approximation can be misleading. To see this 
consider the case where the particle starts at the unstable critical point with a 
finite initial velocity. The movement when restricted to the unstable manifold 
will evidently correspond to revolutions about vertical great circles. Conserva- 
tion of energy implies that the particle will always have enough energy to rise 
up to the saddle (point) and pass over the top. The actual motion however will 
necessarily excite the radial degree of freedom. If the dynamics conspire such 
that as the particle approaches the saddle point enough energy has been trans- 
ferred to the radial degree of freedom to energetically prohibit the particle to 
pass over the saddle, the subsequent motion will fall back down the sphere 
on the same side in complete disaccord with the prediction of the truncated 
dynamics. We have actually numerically discovered this kind of deviation from 
the expected behaviour from the truncated dynamics (see figure pi]) . It is clear 
that the motion is very sensitive to the initial conditions and to the coupling 
between the low energy degrees of freedom and the high energy modes which 
govern the transfer of energy between these modes. 

Another situation where the motion on the truncated manifold is very suscep- 
tible to small transverse oscillations occurs when the truncated set of configu- 
rations ceases to be a manifold. For example if the potential in ( |270| ) is of the 
form 

V{x') = X{g{x')f{x')f, (297) 

as A ^ -|-oo the solution of 

g{x') = (298) 

and 

fix') = (299) 

will give in general two different sub-manifolds for the truncated dynamics. 
This is fine if the two manifolds are disjoint; then they will be separated by 
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Fig. 24. Trajectories of the particle energetically constrained near a circle, and also 
subject to a linear, vertical gravitational potential (left) and a plot of the polar 
angle ^ as a function of time t during the process (right). For certain values of the 
parameters, the particle excites the radial mode sufficiently and does not quite make 
it over the top (even though it has more than enough energy to do so), in complete 
contradiction with the truncation method. Figures shown are obtained for A = 3 
and 6 = 1 (top), and for A = 1 and 6=1 (bottom). 
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large energy barriers and the dynamics will behave independently in each sub- 
manifold. However if they are tangent or even cross, the dynamics becomes 
very sensitive to transverse oscillations. 



Take for example the two dimensional Lagrangian 



n 2 



(300) 



As A — s> +00 the low energy dynamics will truncate to the curves 

y = (301) 



or 



x^+[y-]^' = 1. (302) 



These curves describe the x axis and a circle of radius 1 which is just tangent 
to the X axis at the origin (see figure ^). A low energy particle, in general, 
will effect tiny oscillations transverse to these curves and translate along these 
curves at roughly constant velocity. Suppose we start at a point on the circle 
with velocity tangent to the circle. Depending on the phase and amplitude of 
the radial oscillation that will necessarily be excited, at the moment that the 
particle passes the point of contact, the particle can easily be trapped by the 
valley along the x axis and move on to x — +00. Conversely, a particle moving 
along the x axis, having a slight transverse oscillation, can be trapped in the 
circle. It is clear that the large scale low energy dynamics is not independent 
of the excitation of the high energy, modes however small they may be (see 
figure |2B) . 



4-3 Application of M anion's method to the B = 2 sector of the Skyrme model 



We now apply Manton's formalism of gradient flow curves to the Skyrme 
model, for which the method was first put forward ||37[|. The first step is to single 
out the relevant critical points of the energy functional (in fact, the critical 
points of the potential energy of the model). Unlike the simple examples of 
the previous subsection, this task is quite hard, since a field theory is a system 
with an infinite number of degrees of freedom living (in this case) in a 3 
dimensional space. Furthermore, non-linearity complicates the problem, and 
it is not surprising that coordinated numerical and analytical methods have 
to be used to find solutions to the equations of motion and to study their 
stability. We then describe how the global and relative degrees of freedom 
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Fig. 25. Three-dimensional (top) and contour plot (bottom) of the potential of the 
Lagrangian ( |30CD for A = 20 and a = 1. 



factorize. Finally we present the construction of the full unstable manifold 
linking together the critical points of this sector. 
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Fig. 26. Trajectories of a particle moving on a space made of a straightline with 
a circle tangent to (left) and a plot of its position along the x axis as a function 
of time (right). The final state is extremely sensitive to initial conditions, with the 
particle either going through the intersection area at the first try (top) or getting 
caught in the loop and retracing its steps on the x axis (bottom). 

4-3.1 Critical points of the B = 2 sector 

The first solution of the Skyrme Lagrangian in the B = 2 sector was aheady 
found by Skyrmel^] in the 60's. It is the product ansatz defined as follows: 



Upa{x) = AlUs{x-Ri)A, AlUs{x-R2)A, 



(303) 



where Us is the Skyrmion field, Ai and A2 are SU(2) matrices defining the 
orientation of Skyrmion 1 and 2 respectively, and Ri and R2 represent their 
positions. The product ansatz is a solution of the equations of motion only in 
the limit where the distance between the particles R = \ \Ri — i?2|| is infinite, 
since then, in the neighborhood of each Skyrmion, the other field is the identity 
and hence the product satisfies the equations of motion (and vice versa in the 
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neighbourhood of the other Skyrmion). When R is finite, the product ansatz 
is not a solution of the theory, but it still represents a good approximation 
to the configurations along gradient fiow curves for large separations. This 
is true even though "Skyrme matter", like typical soliton fields, is quite soft 
and the mutual interactions of the fields deform the Skyrmions at a distance: 
they are no longer exact hedgehogs. However the ansatz is good enough to 
obtain an approximation to the potential between the solitons. This potential 
depends strongly on the relative orientation ^41^42 of the Skyrmions, and this 
allows it to be repulsive or attractive. Hence the solution at i? ^ +00 is a 
saddle point of the energy functional and necessarily a state of lower energy 
of most likely different structure exists. We will return to the potential and 
the product ansatz later since they are useful to understand the structure of 
the manifold of gradient flow curves of the system. 

One expects to flnd in the B = 2 sector a state which represents the deuteron, 
and indeed it would be another great triumph of the Skyrme model if such a 
state was found, with the right quantum numbers, binding energy, etc. The 
flrst localized state with baryon number two, which did not give the deuteron, 
is an immediate generalisation of the B = 1 Skyrmion and was probably 
already contemplated by Skyrme himself in the 60 's, although the result was 
first pubhshed by Jackson||88| in the 80's. The B = 2 hedgehog is defined (as 
we have seen in section 3.2.3) as 

U{x) = e'f^^'')'" (304) 



where the profile function /2(r) has to be computed numerically like in the 
B = 1 sector, but obeys the following modified boundary conditions: /2(0) = 
27r and /2(+oo) = 0. This solution is a straightfoward generalization of the 
B = 1 hedgehog, and we will call it the dibaryon in what follows to avoid any 
confusion with the ordinary hedgehog. The mass of the dibaryon happens to 
be about 1.855 x 247r^ in natural units, which translates roughly to 2.5 to 4.4 
GeV (the value depends on the ones chosen for the parameters of the model 
and e) or about 3 times the mass of a single Skyrmion. This state can 
be shown to be unstable. It has much more energy than the product ansatz 
itself, and one expects the dibaryon to disintegrate at least into a pair of 
widely separated Skyrmions. We will see that exactly this has been observed 
in numerical simulations ||91||. Hence the dibaryon is also a saddle point. Thus 
we have, so far, two types of saddle points from which to construct gradient 
fiow curves. We are still missing, however, the minimum energy B = 2 state 
to where these curves lead. 



The state which is generally accepted as having the lowest energy in the B = 2 
sector was discovered independently by Verbaarschot [p^ ] and Kopeliovich||54 
using numerical methods, and actually proposed by Manton||53| by indirect 
methods related to dipole-dipole dynamics and symmetries. Its energy is about 
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Fig. 27. Numerical simulation of the scattering of two Skyrmions via the B = 2 
torus. 

1.18 X 247r^ in natural units which is about 4% less than a pair of free 
Skyrmions. Quite intriguingly, the minimum energy state is toroidal. The 
deuteron, being the least bound of all nuclei, comprises of a proton and neu- 
tron that are quite spatially separated and distinct most of the time. The 
toroidal form well represents this spatial separation, however, it fails to repro- 
duce the distinct character of the nucleons since the Skyrmions are completely 
deformed and have lost their separate identities. On the other hand, this state 
is expected to take an active part in the scattering processes of two Skyrmions, 
especially those with zero impact parameter and with the relative orientation 
which generates maximum attraction. In such a situation, the Skyrmions move 
toward each other, deform and come together in the toroidal state before di- 
viding into two Skyrmions which depart along a direction perpendicular to the 
initial one (the scattering plane is fixed by the initial (iso) orientations of the 
Skyrmions) Il53| , p2| (see figure ^ taken from reference |9^). This 90° scatter- 



ing is typical of head on collisions of identical solitons[^. The interpretation 
of this configuration as the deuteron is still controversial given its large bind- 
ing energy (of the order of 100 MeV instead of the 2 MeV or so observed in 
nature), small radius and "strange shape" .[^ Although it is not certain that 



^ One way to resolve this issue might be to compare the cross section for nucleon- 
nucleon scattering computed from the Skyrme model with actual experimental mea- 
surements, and look for a signature predicted by the Skyrme model for Skyrmions 
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this state is the absolute lowest energy state of the B = 2 sector, it does seem 
reasonable that it is so. No other states with lower energy have been found so 
far, numerical simulations of Skyrmion-Skyrmion scattering have not uncov- 
ered another intermediate state |9^, and numerous modes of oscillation around 
the toroid have been considered without finding any negative modes[^. We 
will assume here that it is in fact the state with lowest energy and that the 
manifold of low energy dynamics of the B = 2 sector consists of path joining 
the minimum and the two different families of saddle points. 



To visualize the unstable manifold of the i? = 2 we need to know how the 
different configurations presented above are linked together by the dynamics. 
The number of degrees of freedom in the unstable and zero modes of each state 
will be most helpful. A solitonic system has an infinite number of degrees of 
freedom or modes that it can excite. Modes can be categorized according to 
their role in the stability of a particular configuration: negative modes gen- 
erate paths followed by the system when it moves to another configuration 
of lower energy, positive modes are modes associated with oscillations around 
soliton configurations and zero modes correspond to rigid motions of the soli- 
tons arising from symmetries of the energy functional which are not respected 
by a given solution. If the system has low enough energy, only a finite number 
of these modes will be excited. In fact the system will first explore only the 
zero and negative modes (if there are any); higher energy excites the positive 
modes starting with those with the lowest frequencies. The parameters corre- 
sponding to zero modes are particularly important, and are called collective 
coordinates of the system. They form what is called the moduli space for the 
case where there are no negative modes. For example, in any BPS system, 
the collective coordinates corresponding to the zero modes are sufficient to 
describe the low energy motion of the system, which then corresponds to the 
geodesies on the moduli space. In cases where there is a weak potential acting 
between the solitons (like the Skyrme model), one has to also include the neg- 
ative modes and the union of the corresponding gradient flow curves, to have 
enough freedom to let the system evolve in a satisfactory fashion. The num- 
ber of degrees of freedom rises with the baryon number, since new collective 
coordinates enter to describe the relative motion of the particles, however the 
number can decrease if new symmetries arise via the dynamical evolution of 
the system. 



going through the toroidal configuration. This is still an open problem since no 
general cross section for nucleon-nucleon scattering, neither classical nor quantum 
mechanical, has been computed so far (We present in the last subsections of this re- 
view our own semi-classical computations of nucleon-nucleon scattering angles with 
the Skyrme model under certain conditions.) 
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4-3.2 Factorization of global and relative coordinates 



Let U (x) be a completely general configuration of the system obeying the con- 
ditions that U{x) goes to the identity at infinity fast enough so that the total 
energy associated to the configuration is finite, localized in space, and the 
baryon number associated to the configuration is an integer. Such a configu- 
ration possesses generic collective coordinates corresponding to the invariance 
of the energy functional under the following 3 transformations: 



x^x + Rt (305) 
x^D{B)-x (306) 
U{x)^AU{x)A'^. (307) 

In general a configuration is not invariant under these transformations, giving 
rise to nine zero modes in the fluctuation spectrum. i?T is the position of 
the center of mass of the system, and is related to invariance of the energy 
under translation, the rotation matrix D{B) associated to the SU(2) matrix 
B parametrizes the rotational invariance of the system under global rotation, 
and finally the SU(2) matrix A parametrizes the invariance of the system 
under global isospin rotation. Because of the absence of special symmetries of 
U (x), the 9 parameters are completely independent, and if one gives them time 
dependence and quantizes the system treating them as dynamical variables, 
they give, respectively, the conserved total momentum P, the conserved total 
angular momentum L and the conserved total isospin T. The 9 global collective 
coordinates do not participate in the interesting dynamics of the system. 

We will now apply these ideas to the B = 2 sector of the Skyrme model in 
turn considering the product ansatz, the dibaryon and the torus. As we saw 
earlier, the product ansatz 

Upa(x) = AiUs{x-Ri)Al A2Us{x-R2)Al (308) 



is a solution of the theory only in the limit where R = \ \Ri — -R2II is infinite; 
we will consider it to be arbitrarily large. The product ansatz is a solution 
described by 12 collective coordinates, or twice what is needed to describe 
a single Skyrmion: 3 parameters to define the position of the Skyrmion, and 
another 3 for its orientation or iso-orientation. This last fact comes from the 
spherical symmetry of the Skyrmion: 



AUs{x)A^ ^ cos f{r) + [sm f{r) At A^ ■ x 

= cos /(r) + i sin f{r)T^D^y,{A)x'' (309) 



91 



hence 



AUsix)A^ = UsiD{A) ■ x) (310) 



where D{A) is the rotation matrix defined by (|T^ associated to the SU(2) 
matrix A. So iso-rotating the Skyrmion by a certain amount is exactly equiva- 
lent to rotating it in space by the same amount. Thus one requires only three 
parameters instead of six to describe the orientation (and iso-orientation). 
Another way to see this is to consider the following expression: 

U{x;A,B)=AUs{D{B)-^-x)A^ (311) 



where the iso-orientation A and the orientation B are considered independent. 
This state is invariant under the following transformation of A and B: 

A^AC (312) 
B^BC (313) 

since 

f/(f; AC, BC)=ACUs{D{BC)-^ ■ x)C^A^ 
= AUs{D{B)-^ ■ x)A^ 

= U{x;A,B) (314) 

by equation ( |310| ). So U{x;A,B) has a continuous redundancy parametrized 
by 3 angles, which arises from the spherical symmetry of the Skyrmion. The 
remaining 3 angles and 3 coordinates of the center of mass parametrize a 
manifold which can be written x IR^. The product ansatz therefore has 6 
independent parameters to describe each of the two particles, for a total of 12. 

We will now show how to re-express these 12 parameters Ai, A2, Ri and R2 in 
terms of global and relative ones. To do this we will reproduce the discussion 



of Walhout and Wambach||2CI||, who use the analogy with the treatment of a 
rigid body: the global coordinates will represent the transformation from the 
laboratory frame to the body fixed frame, and the relative coordinates will 
describe the system in the body fixed frame. Walhout and Wambach[20| chose 



the Skyrmions to be separated along the z axis in the body fixed frame by a 
distance R, and the Skyrmion at the position x = —{R/2) 63 to be isorotated 
by C relative to the other: 

AY = 1 (315) 
A'Y = C (316) 
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Ri = (317) 
R2 = ~es. (318) 

In the lab frame, we will use Rt to denote the position of the center of mass 
of the system, and we will position the Skyrmions in space using the inverse 
of the rotation matrix corresponding to the SU(2) matrix B, and a matrix A 
will rotate the whole system in isospace. This gives us the product ansatz 



Ulf{x)=A,Us{x-R,)A{A2Us{x-R2)Al 

= AUs[D{B) ■ (f - i?T - D{B)-' ■ |e3)] x 

CUs[D{B) ■ (f - i?T + D{B)-' ■ ^es)]C^A^ (319) 



with 



Ai = AB (320) 
A2 = ACB (321) 

R, = RT + D{B)-'-^e3 (322) 
R^ = R^-D{B)-^-^e^ (323) 

which is clearly of the form: 

U{x) = AU{D{B) ■ {x - Rt))A^ (324) 

where the global degrees of freedom are singled out. If we count the parameters 
we find that there are 13 of them instead of just the 12 which we started with. 
Walhout and Wambach||20|| (inspired by some work by Verbaarschot|9^) show 



that one angle, which we shall define shortly, parametrizes a symmetry of the 
ansatz ( |319D and is therefore redundant. In the body fixed system, the product 
ansatz (|319|) obeys the following symmetry: 

U^l{x;C) = C(ir-n)U^^(D(ir-n)-x;C)(-ir-n)C^ (325) 



where n is any vector perpendicular to the z axis, and C is defined as follows: 
C{f-n) = {T-h)Cl (326) 

If one writes C as e^'^"^, then this implies that n = 7 x 63. Let us check that 
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( |325|) is indeed a symmetry of the product ansatz. Applying the transformation 
to the product ansatz in the body fixed frame one gets: 



U^iix;C) = Usix- |e3)Cf/s(f + |e3)Ct 

-^CUs{x+ |e3)Ctf/s(x - |e3) (327) 

since D{iT ■ h) is a rotation by vr around an axis perpendicular to the z axis, 
and using ( ^26| ). We then see that the Skyrmions change places, but in the 
limit where R — > +00, which is a necessary condition for the product ansatz 
to be a solution of the equations of motion, the order of the Skyrmion matrices 



is irrelevant since the two commute. Then ( |325| ) is really a symmetry of the 
solution within the product ansatz. It is also a symmetry of the Lagrangian 
(since it only consists of global rotation and isorotation) and of the subse- 



quent evolution. The symmetry (|325|) is valid for any vector n as long as it is 



perpendicular to the z axis: the angle which defines the orientation of h in the 
(xy) plane actually spans a whole family of discrete symmetries of the ansatz. 
Since n = 7 x 63, rotating n in the (xy) plane by the angle 6 also rotates 7 by 
the same amount in the same plane, which in turn is equivalent to acting on 
the matrix C as follows: 

C ^ e^^'Ce-'^'. (328) 

It is easily checked that for the product ansatz this symmetry takes the fol- 
lowing form: 

U^p{x- C) = e-'^'f/^fp(D(e'^') ■ f; e'^'Ce-'^')e'^' (329) 



or that, in the laboratory reference frame, the product ansatz is invariant 
under the combined transformations: 



A^Ae-'-^^ (330) 
B^e'^^'B (331) 
C^e'^'Ce-'^'. (332) 

This shows that the angle of rotation 6 around the z axis implicitly contained 
in C is redundant and does not play a role in the describing the configuration. 
C is then only parametrized by 2 angles: one for the orientation of 7 in the 
(xz) plane, and the magnitude I7I =7, by which the second Skyrmion is 
rotated relative to the other. We then have, as we should, 9 global collective 
coordinates {A, B and ^t) and 3 relative ones (C and R) giving a total 
of 12 instead of 13. One should carefully note that equation (|329D by no 
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means implies that the product ansatz possesses a continuous symmetry: it 
only does if one writes it using too many parameters as collective coordinates 
(exactly like the case for the ordinary Skyrmion described above by equation 
(|314|) ). Once the correct number of collective coordinates has been used, the 
only symmetries that the product ansatz possesses are 3 reflection symmetries 
relative to the [xy), {yz) and {xz) planes (all those discrete symmetries are 
exact only in the limit where R —>■ +00). 



4-3.3 The construction of the unstable manifold 

Let us suppose that the Skyrmions always keep the same hedgehog shape 
no matter how strongly they interact with each other. Then one could ob- 
tain the Lagrangian of the system by giving time dependence to the collec- 
tive coordinates, and replacing the product ansatz in the Skyrme Lagrangian 
and computing the requisite integration over all space. These 12 degrees of 
freedom then describe a manifold, which we shall note M12, with a metric 
and a potential induced by the kinetic and potential parts of the Lagrangian, 
respectively [p5|]. The Skyrmion-Skyrmion or nucleon-nucleon dynamics is then 
obtained by doing classical or quantum mechanics in the curved space of the 
manifold M12. Unfortunately, this program is not correct since the Skyrmions 
deform and do not stay in the form of a hedgehog when they interact (obvi- 
ously since they deform enough to merge into a torus), and M12 has no hope 
of describing the full dynamics in the B = 2 sector. But one expects that the 
unstable manifold {i.e. that obtained by gradient flow curves) of the B = 2 
sector, which we shall call Aiu, and M12 will become approximately equal for 
R sufficiently large. 

Numerical studies have shown that the product ansatz give a surprisingly good 
approximation to A^i2 up to a separation of just a few fm (depending on the 
initial relative orientations of the Skyrmions). This enables one to extract use- 
ful information on the Skyrmion-Skyrmion dynamics relatively easily. One can 
for instance compute an approximation to the Skyrmion-Skyrmion potential 
even for moderate values of the separation R. The separation is not a zero 
mode of the energy any longer, but it is a relatively "soft" mode such that the 
potential energy does not vary by large values as a function of R, as long as 
R is large. In low energy scattering, the R mode can be considered a "slow" 
mode compared to the angular variables representing the rotation of particles 
which are "fast" variables. We will describe exactly such an approach in the 
next subsection. 

The first calculation of the Skyrmion-Skyrmion potential was obtained by 
Skyrme himself and is defined by: 

V'ss = ^B=2 - 2Ms (333) 
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where Eb=2 is the static energy of the system. This involves integrating Skyrmion 
profiles over all space, which is best done using an expansion in inverse powers 
of the separation R. To leading order one finds 



VssiR, C) = -Anfy ' ^" ^3 ^ ^'^ (334) 



which only depends on the relative coordinates: 2 angles and the distance R. 

Subsequently, there has been a great deal of work on the extraction of the 
nucleon-nucleon potential from the Skyrme model and its comparison with 
traditional nuclear potentials Early computations, mainly using 
the product ansatz, quickly showed that the general tensorial form of the po- 
tential so obtained was in good agreement with what is known of the nuclear 
force (1-2 pion exchange, repulsive core, etc.). However, a particularly dis- 
turbing defect haunted the problem for many years: the absence of a central 
attraction at medium range. This was all the more disturbing since it is pre- 
cisely this part of the nuclear force which binds nucleons together into nuclei. 
Many ansatze were tried generalizing the product ansatz: the modified prod- 
uct ansatz (the Skyrmions were given the freedom to change their radii as 
they interacted, but with little gain in the central attraction), the sjTiimetrized 
product ansatz]^ (a step in the right direction since it was symmetric un- 
der the exchange of the particles, and provided some central attraction but 
at the expense of the long range behaviour), and more recently the instanton 
methodp5| (which combines the advantages of all the previous ansatze), but 
the crucial central attraction still eluded all efforts. 

It was in two[^6|,0 exact, numerical analyses that a central attraction was 
found. The two approaches bifurcated in the way that they included physi- 
cal, quantum nucleonic states. In the approach of Verbaarschot et al ||56[ the 
potential is computed using numerical relaxation methods on a 20 x 20 x 40 
regular lattice for fixed separation and imposing only the discrete symmetries 
of the most attractive channel ( see below). The separation was defined to be 
the distance between the respective centers of baryon number density. They 
found a reasonable attraction (roughly 70 MeV) in the central channel. In 
their approach, the physical nucleonic states where incorporated by treating 
the angular and iso-angular coordinates as rigid quantum rotors, and the ra- 
dial motion was treated semi-classically using the WKB method. They found 
that the kinetic energy of the angular and iso-angular motion was sufficient to 
destabilize any possible bound state. However the numerical analysis was re- 
cently repeated using a more refined method P5[], the finite element approach, 
which permits an irregular lattice. In this way, for the same amount of com- 
puting time, the algorithm can sample the critical regions with a finer grid, 
and consequently improve the accuracy. This work shows that indeed the cen- 
tral channel is sufficiently attractive to support a bound state, the deuteron. 
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within the WKB method for the radial degree of freedom. On the other hand 



in the work of Walhout and Wambach ||97|| , although the potential was again 



computed for static Skyrmionic configurations, the nucleon-nucleon potential 
was extracted by projection onto asymptotic quantum nucleon states a la 
Jackson, Jackson and Pasquier|p. Here also a central attraction was found, 
however the further analysis necessary to establish a bound state was not pre- 
sented. Walet and Amado [PDf] further improved the results by including the 
A resonance as intermediate state, as well as some gluonic corrections. With 
all these corrections, the thus obtained Nucleon-Nucleon potential is indeed 
quite close to traditional (phenomenological) nuclear potentials, even though 
the arbitrariness in the definition of the distance between a pair of particles 
when they are very deformed (at small R) makes it hard to judge. Here we 
will only be interested in the long range part of the potential which takes the 
form of the one-pion exchange potential. 



The tensorial form of the potential reveals how Skyrmions will react to their 
relative orientations. The different situations can be catalogued in three cases 
or channels. In the hedgehog-hedgehog channel (HH), C = 1 and the potential 
is zero (no one-pion exchange potential) . As R gets smaller then other interac- 
tions take over and build a repulsive core. In this channel, the Skyrmions stay 
roughly hedgehog-like until R is of the order of 1 fm. If the second Skyrmion 
is rotated by vr around the axis of separation, C = ir^, then a strong repulsive 
interaction keeps the Skyrmions apart and the potential energy rises quickly 
as R decreases. In this channel (REP), Skyrmions only stay roughly spheri- 
cally symmetric for R greater than about 1.5 fm. The most interesting case 
is the so-called most attractive channel (MAC), where one Skyrmion is ro- 
tated relative to the other by vr around a direction perpendicular to the axis 
of separation. In this channel, the Skyrmions always attract and come closer 
and closer together until they fuse into the toroidal configuration. This is the 
most studied and phenomenologically interesting case since it is believed to 
contain the nucleon-nucleon binding dynamics. Of course, in a realistic gen- 
eral collision of nucleons, one does not stay in any one of these channels (spins 
are generally not correlated during real scattering) and situations should be 
hybrids of the ones described so far. 

By freezing the relative iso-orientations, the two collective coordinates rep- 
resented by C disappear, and the configurations then span a 10-dimensional 
manifold: Mf^, Mf^^ and MfJ^'^. Mf^ and Mfo^ are actually not very 
interesting since they would not represent very accurately the dynamics of the 
system. This is because the HH and REP channels lead to relatively energetic 
states: the gradient on the manifold is large and the probability of exciting 
modes perpendicular to Aif^ and A^^f^ is not negligeable. It seems clear to 
us that more than 10 collective coordinates are needed to describe the system 
in these channels. On the other hand, it seems reasonable that if the system 
has low enough energy, it will all by itself steer clear of these regions and will 
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not "feel" the existence of these higher energy states. It also seems physically 
sound that the system will spend most of its time near Ai^^'-^. There the 
energy gradients are much smaller (there is about a 4% difference in energy 
between the torus state and a pair of infinitely separated hedgehogs). The 
manifold A4^^'^ has been extensively studied by Leese et al ||5^ in connection 
with the study of the deuteron system. 

The product ansatz in the most attractive channel possesses the following 
discrete symmetries at infinite separation (for two well separated Skyrmions 
on the z axis, equidistant from the origin with one rotated relative to the other 
by 180° about the y axis): 



7T{—x,y,z) = 10 Tf{x,y,z) (335) 




1 0~ 

7i{x,-y,z) = \ -1 0\7i{x,y,z) (336) 

1 

-1 0~ 

Ti{x,y,-z)=\ 1 0\7i{x,y,z). (337) 

1 



The last symmetry is actually valid for any separation while the first two are 
approximate and are exact only in the R ^ +oo limit. These discrete sym- 
metries have been shown to be exact for all configurations of A4^^^ including 
the torus, because they are conserved by the dynamics starting from R = +oo 
in the most attractive channel. 

As the Skyrmions merge into the torus, the separation coordinate R ap- 
proaches a minimum value which depends on the definition for the separation 
chosen (which is always arbitrary to a certain point, since when the Skyrmions 
deform and come close together, the definition of separation becomes blurred). 
Another collective coordinate is effectively eliminated by the appearance of the 
axial symmetry of the torus [p.OO|] : 

U,oUS) = e^^-Vto.us(/^(e-^^') ■ f)e-^-^ (338) 



This axial symmetry is created by the dynamics of the system and generalises 
the discrete symmetries of equations ( p35| ), ( |336| ) and ( ^37| ). This implies that 
the torus is described by only 8 collective coordinates which parametrize an 
8-dimensional manifold named A4s- Manton shows that one can describe the 
topology of A^8 somewhat more precisely. Let us factor out global translations 
and global isorotations from Aig, which are represented by a factor of IR^ x 
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S^. We are then left with 2 Euler angles, which parametrize a 2-sphere. The 
appearance of such an S"^ on a manifold is called a Bolt, in the theory of 



gravity [p6| , |3T| , |l 1|] . A Bolt is a kind of "soft" singularity on a manifold. Atiyah 
and Manton further argue that because of the reflection symmetry of the torus 
relative to the plane perpendicular to the axis of the torus (the plane which 
"slices" the torus in two identical halves modulo a reflexion), the 2-sphere is 
really an IRP2: a sphere with its antipodal points identified. This means that 
the two Euler angles over-define the orientation of the axis of the torus. This 
is topologically similar to the BPS case where a torus state also exists in the 
winding number (magnetic charge) 2 sector. 

The dibaryon is at the center of the most exotic dynamics of the B = 2 sector. 
Because of its hedgehog structure, see equation ( p04| ), it has the same 6 zero 
modes as the B = 1 Skyrmion. Like all solutions of the theory, it also has 
positive modes corresponding to small oscillations about the critical point. 
However, it also possesses negative modes since it is a saddle point of the the- 
ory. It has long been suspected that the dibaryon is able to disintegrate into a 
pair of Skyrmions or into the torus itself; this was numerically checked only re- 



cently by Waindzoch and Wambach [pi|] . In their articles they use a discretized 



version of the dibaryon as the initial configuration and then perturbed it in or- 
der to study how the dibaryon is connected to the rest of the B = 2 manifold. 
The negative modes of the dibaryon had already been investigated analytically 



by Bang and Wirzba||51[] on a 3-sphere of radius L. The limit L — > -|-oo of their 
findings is in good agreement with those of the numerical simulations. They 
solve a Schrodinger-type equation for the perturbation of the dibaryon field in 
order to find the perturbations 6U which give the maximum negative energy 
gradients. This gave them three magnetic modes and three electric modes, 
whose names refer to the properties of transformation of the perturbation un- 
der rotations. We have already discussed these negative modes in section 3.2.3. 



The magnetic mode along the direction z is parametrized as follows|91]: 



SU^ : 57r° = (339) 
Sn^ = -g{r)^ (340) 

Srt^ = g(r)- (341) 
r 

671^ = (342) 

where g{r) is a function computed numerically, the modes along the x and 
y axes being obtained by cyclic permutations. This mode possesses a simple 
interpretation when one recognizes that the dibaryon configuration is very 
close to a pair of Skyrmions in the product ansatz, placed one on top of the 
other. There are two ways to lower the energy of the resulting configuration: 
by translating or by rotating one Skyrmion relative to the other. The former 
represents the magnetic mode; there are three independent magnetic modes 
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corresponding to the three orthogonal translation directions present in R^. 
The latter represents the electric mode which is parametrized by a somewhat 
more complicated form||91||: 



5f/f :5vr° = -3a(r)- (343) 

(344) 

..^^(cM-M)i| (345) 
S^'-3Ur)-^-^y^ + lb(r) (346) 

for a perturbation along the z axis (the profile functions a(r), b{r) and c(r) are 
obtained numerically). The other two can be obtained by cyclic permutations, 
6U^ and 61/^ corresponding to rotations along the x and y axis respectively. 

Ordinary time evolution of the dibaryon subject to these modes exhibits a fis- 
sion process taking place. It can be numerically seen that for the perturbation 
6U^, the dibaryon disintegrates in two Skyrmions moving along the z axis. 
This negative mode is associated with a large energy gradient: the Skyrmions 
separate relatively quickly and oscillate strongly. The electric mode SU^ is 
slightly less energetic and the Skyrmions although moving away also along 
the z axis, do so in a "twisting" fashion. Large oscillations are also present 
here. Waindzoch and Wambach also studied these same disintegrations using 
the gradient flow method. The oscillations, which are transverse to the low 
energy unstable manifold, are damped (see equations (|261| ) and ( |262| )) and it 
appears then that the magnetic mode is tangent to the path leading to the 
product ansatz in the HH channel, while the electric mode leads to the REP 
channel. This is still a delicate issue since the numerical computations are very 
hard to do and the precision required to confirm these results is quite high||91|]. 



Confirmation by analytic methods would be most welcome. Waindzoch and 
Wambach also show that the right superposition of perpendicular magnetic 
and electric perturbations can direct the system along the path leading straight 
to the toroidal configuration. 

We are now ready to put the various parts together. Figure ^ shows 
schematically what the configuration space for the B = 2 sector looks like. 
The 12-dimensional manifold A4i2 is the minimum required to describe the 
dynamics of a pair of Skyrmions or nucleons. As we saw earlier, for Skyrmions 
infinitely far from each other, these 12 collective coordinates can be identified 
as the individual positions and iso-orientations of the separate solitons. These 
states can be well represented by the product ansatz saddle point. The product 
ansatz comes in three different iso-orientations, if one chooses to freeze such 



100 



Dibaryon, E ~ SMg 




Deuteron, E ~ 1.92Ms 



Fig. 28. Diagram of the B = 2 manifold showing paths in configuration space 
connecting the dibaryon to the product ansatz in the HH and REP channels with 
paths tangent to the magnetic (M) and electric (E) modes respectively. The deuteron 
and the product ansatz MAC are connected to the dibaryon by paths comprising 
of an orthogonal superposition of E and M modes. 

a degree of freedom: hedgehog-hedgehog channel, repulsive channel and the 
most attractive channel whose collective coordinates parametrize the mani- 
folds AA^^I^'^ respectively. The product ansatz in the most 
attractive channel is linked to the toroidal configuration and its manifold M.^ 
by gradient flow curves, while the HH and REP channels are connected to 
the dibaryon by paths whose directions near the dibaryon are tangent to its 
magnetic and electric modes respectively. The number of collective coordi- 
nates along those paths changes from 10 to 6 because of the appearance of 
the spherical symmetry of the dibaryon as well the freezing of the separa- 
tion between the particles, removing 4 collective coordinates. The dibaryon is 
also directly connected to the torus as mentioned earlier by paths with initial 
directions spanned by the superposition of orthogonal magnetic and electric 
modes. Intermediate paths are generated by different combinations of these 6 
modes. At the current rate of progress in this area, the relationship between 
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the initial conditions and the final states attained should soon be clear. In the 
adiabatic approximation for Skyrmion dynamics, the paths can be followed in 
any direction (there is no pion wave emission here) but this is not physically 
plausible when energies differences of more than 1 GeV are concerned. 



4-4 Low energy nucleon-nudeon scattering 

We now turn to the problem of nucleon-nucleon scattering as reported in 



Gisiger and Paranj ape p9|j40|| following the methods proposed by Manton|37 



Scattering of Skyrmions at low energies should truncate to a dynamics taking 
place on an unstable manifold linking together all of the low energy critical 
points. We consider the scattering in a particularly simplifying approximation. 
One observes that as the Skyrmions are separated to infinity, their low-energy 
behaviour corresponds to two independent, free Skyrmions, each with 3 trans- 
lational degrees of freedom and 3 rotational or isorotational degrees of freedom. 
The moduli space is simply two copies of IR'^ x S^/Z2 and the induced metric 
is also the natural metric on these manifolds. As the separation between the 
Skyrmions decreases, interactions develop between the Skyrmions, both in the 
potential and in the kinetic energy {,i.e. , the metric). These interactions can 
be written in a systematic perturbative expansion in inverse powers of the 
separation. In what follows, we will only consider scattering at large impact 
parameter, for which the separation between the Skyrmions is always large. 
The perturbative expansion can therefore be truncated, and indeed our main 
approximation is to take only the leading term in this expansion. Surprisingly, 
this term comes from the expansion of the kinetic part of the Lagrangian and 
not the potential, and it has actually been overlooked by previous investi- 
gations. A secondary assumption that has been made and which we want to 
underline is that the pion mass is taken to be zero. This could perhaps explain 
why the interaction that we find has been missed: it comes from the two-pion 
exchange part of the interaction which is generally neglected when considering 
massive pions. We complete our analysis by using the Bohr-Sommerfeld quan- 
tization method and a modification of the method of the variation of constants 
to calculate the scattering trajectories of the nucleons. 



4-4- i Lagrangian of the Skyrmion- Skyrmion system 

Written using the left-invariant one form already introduced previously (equa- 
tion (0)), the Skyrme Lagrangian takes the form: 



f 1 
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n{u)^-n{uyn{u),-n{u)^ 



(347) 



which, as we saw earher, separates into a kinetic and a potential part: 

7^o ■ 7^i 7^o ■ 7^i - 7^o • 7^o n-, ■ n-. 



P 1 



/•2 1 

V = -—TZ- ■ TZ 

^ 2 ' ' 4e2 



7^i ■ 7^i 7^j • 7^j - 7^i • 7^j 7^j ■ 7^i 



(348) 
(349) 



The potential part has been extensively studied in the literature. It pro- 
duces the following leading order term in the separation d between the 
Skyrmions0,|,|5]: 



V = 2M + 47r/> 



2 2(l-COS^)(3(n-rf)2-l) 



rf3 



(350) 



where d = d/d. We are interested by the kinetic part of the Lagrangian which 
is less understood. 

As noted earlier, to accurately describe the motion of the pair of Skyrmions 
we need the metric on the unstable manifold of the baryon number 2 sector, 
as well as the potential defined on it. Since both are not well known one 
must resort to approximating the pair of Skyrmions by a parametrization or 
ansatz. We chose the simplest parametrization for the system, the product 
ansatz, which maintains both solitons rigidly in the Skyrmion configuration 
at all time. Of course, since the particles deform when they come close to each 
other, our parametrization is only valid when the Skyrmions are far from each 
other. This is compatible with the low energy assumption we made and which 
is necessary to ensure that the degrees of freedom of the system can indeed 
be truncated to a finite (small) number. Then one takes 

U{x) = Ul U2 = A^U{x - Ri)A B^U{x - R2)B. (351) 



A and B are time dependent SU(2) matrices representing the orientation of 
the Skyrmions and Ri and R2 are their respective positions. A symmetrized 
product ansatz is possible, which respects the symmetry under exchange of the 
two particles[0, however such an elaboration does not affect the first order 
term in the mutual interaction. 

Replacing U by U\ U2 in the kinetic part of the Skyrme Lagrangian (|T^) and 
working out the computations using the identity 

ni{U,U2) = mU,) + Dab(t/i) Km (352) 
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where DahiU) is \ ti[T^UT^U^], allows us to isolate the interaction term 
from contributions describing free individual particles. The kinetic part of the 
Lagrangian then writes 

T = Ti + T2 + Tint (353) 



where 



1 



+ 



2e2 



■ 7^f ■ 7^f + 1 4-^ 2) + 2{nl ■ + 1 ^ 2)n\ -D-n 



+ 2{n\ ■ 7^f + 1 ^ 2)nl ■ d ■ nl + atzI ■ d ■ nl nl ■ d ■ nf 
- {nl-D- nff - {Til ■ D ■ nlf -2nl- tz\ nl ■ nf 



2nl-D- 7^f n\-D- nl 



-2{nl -n] + 1 ^2)(nl- D -nf + nl ■ D ■ k 



(354) 



where the dot implies contraction over isospin, TZj^ = 7l^{Ui), 7^^ = 71^{U2) 
and D = D^]^{Ui) and Ti and 7^ are similar to (|I^) for the fields Ui and U2 
respectively. 

To obtain the full Lagrangian of the system we need to integrate the La- 
grangian density over all space: 



int- 



(355) 



Since the Skyrmion profile function /(r) is only known numerically, even 
though it falls to zero like for large r, it is not possible to compute 
the integral analytically. An expansion in inverse powers of the separation 
d = — i?2|| between the Skyrmions is then used {d was assumed to be 
large from the start in order for the product ansatz to be accurate). 

To do this we first write Ui and U2 as functions of the matrices A and B and 
the position vectors Ri and R2 using ( |351| ) and the following identity: 



^o(f^i) = (^'' - D^AUix - R^)A^)^ 7^^(A) - 

D,^{A)R\n^{Uix~R,)) (356) 

and correspondingly for U2, B and i?2- We are interested in the leading con- 
tributions in inverse powers of d of these tensors. Since the Skyrmion profile 
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function /(r) behaves like /t/r^ at large separation, we get after a short cal- 
culation for large values of |x — -Ri|: 



_abc 



-2k- 



\x 



Ri\ 



+ 01 



\x 



(357) 
(358) 



and 



\x 



\x 



Ri\ 



(359) 



where fi = (x — Ri)/\x — Ri\, and correspondingly for U2, B and R2. To go 
further we will use a standard method for this type of calculation (for a very 
complete description of this method, its subtleties and its application to the 
computation of the Skyrmion-Skyrmion potential in the product ansatz see 



In this method, we divide space in three regions I, II and C. I and II are 
regions of space close to the Skyrmions 1 and 2 respectively (see figure ^91) . 



We define the notion of closeness by stating that these are regions where the 
profile function /(r) is not well approximated by its asymptotic expression. For 
simplicity, and taking advantage of the spherical symmetry of the Skyrmion, 
we will take those regions to be spheres of some chosen radius R. Region C is 
the complementary region where both Skyrmion profiles fall off like l/{Ri — r)'^ 
(z = 1, 2). We then evaluate the integral over each region separately using 
the fact that in regions I and II, one can use the "exact" expression for the 
Skyrmion profile, while in region C both Skyrmion profiles behave according 
to their asymptotic forms. 

Let us start with the contributions from the non-linear a model, i.e. the term 
quadratic in time derivatives. By using (0) and (|356|) we easily find that: 



7^°(f/l) Dab(f/i) Kiu^) = 

-7^S(A)Pea(t/^) /^ab(t/l) D^B) R^2T^f (U (x - R2)) 

~n^{U{x-R,))R\D,,{A^) Z},b(f/i) VUU2)K{B) 

+ niiU^x - i?i))i?iDca(A^) /^ab(t/i) L'bd(5)i?i7^f ([/(x - R2)) (360) 

When integrating this expression over region I we can use the asymptotic 
expansion for the tensors of Skyrmion 2, giving 
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Region C 



\x — Ri \ > B, \x - 





Fig. 29. Diagram of the three regions I, II and C used to divide space in order to 
compute the Lagrangian of a pair of Skyrmions as an expansion in inverse powers 
of the separation between the particles 



J d^f 7^°(c/l)L'ab(c/l)7^^(c/2) 



-2k J d^x 

\x-Ri\<R 



ni{A)V,^{Ul)D^,{U{)e^-D,,{B)rint{B) 



\x-Ro\^ 



+ 



I. j 



+2K / d'' X 

\x-Ri\<R 

+0 



3 . 7^f(^7(f - 4))4/5ab(At)Z}bc(f/l)e^'^'^i^ef(5)r^7^^(5) 



\x-R2? 



|f-i?o|3 



(361) 



By translating the integration variable by and defining a new integration 
variable y — x — Ri, we easily see that this integral contributes at most terms 
of order Integration over region II yields the same result by symmetry. 

Integrating over the complementary region C is different since tensors from 
both Skyrmions take their asymptotic form and we find to leading order: 
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4^2 / e-^^DAA)r{nt{A)e^''^n^,{B)D,,,{B)r^ ^ 

. ' . |f-i?i|2|x-i?2p 

|x-_Ri|>/?, \x~R2\>R 



i+ii+c 

+ . . . (363) 
= 1^ / ^^o(-4)^o(-g)[-Pac(^)^^/^bd(A)(g- ci) - jz- d)'Dcd{A^ 3)2^5^"^] 

+0(1) (364) 

where y = x — R2 and z = y/d. To obtain ( P63|) , we used the fact that 
expanding the integration bounds from C to the whole space only adds higher 
order contributions to the integral, but does not alter the value of the leading 



order terms (see [35| for details). The last line (364) is obtained by shifting 
the integration variable by Ri and absorbing 3 factors of d in the measure. 
This is possible since we are integrating over all space. 

The fact that the main contribution comes from the faraway region might at 
first sight seem a bit strange, but it should be kept in mind that this same 
region makes a contribution in the analogous computation of the potential 
(13501) , which is as important as the contribution from the regions close to the 



Skyrmions (see [35| for more details on that point). 



Similar computations are feasible for the Skyrme term, but because of its 
structure (quartic in derivatives), it only contributes at higher order. The 
leading contribution from the kinetic part of the Lagrangian is then of order 
1/d while the potential part is of order 1/d^ rendering the latter, in principle, 
negligible compared to the contribution from the kinetic energy. Finally we 
obtain the following expansion of the Lagrangian for a pair of Skyrmions far 
from each other in the product ansatz, using the expression (|2l|) for single 
Skyrmion Lagrangians: 



L = -2 M+^Md"^ + 2A(7^^(A) 7^^(A) + 7^^(5) 7^^(5)) 

+ ^e'^^ej'^'^ C^A) C^iB) {6'^ - Sd)) D,^{A^B) + 0{l/d^) (365) 

where d = d/d, A = 2nK'^fl /:^(A) = ^A), 7^^(A) = 7^g(A) and we have 
used the relation C'iA) = D^^{A) n^{A). 

By keeping only the leading order terms in 1/d, we neglect any contribution 
from the potential part of the original Lagrangian, and the motion is then 
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completely specified by the geodesies of the metric induced by the kinetic 
term onto the unstable manifold of the baryon number two sector, as specified 
in subsection 4.3. 



The term in l/d in (|365|) , though absent from the literature of the 80's, was also 



independently obtained by Schroers||41||. He found a leading contribution which 



behaves as l/c? and even calculated sub- leading spin-orbit coupling terms. The 
only other comparable calculation to our knowledge has been done by Walhout 
and Wambach[pO| for the case of massive pions. The limit as — of their 
expression, however, does not leave a term which behaves as 1/d and hence 
does not reproduce our result. We beheve that this contribution should come 
also from their evaluation of the integral giving the induced kinetic energy 
in the faraway region (region C) and then recovering our result as m,r — ^ 0. 
Such a contribution would also be proportional to {e~"^^'^)'^, what they call 
"two pion exchange" . We also add that this 1/d term is of leading order in an 
expansion in inverse separation with respect to a scale determined by f^^ and 
e which has nothing to do with the length scale set by the pion mass. 

The Skyrmion-Skyrmion Lagrangian we obtained describes dynamics in a 12- 
dimensional moduli space via equations of motion which are highly non-linear 
and quite complex. The Lagrangian possesses several symmetries, and associ- 
ated conserved quantities such as total isospin which is connected to invariance 
under left isorotation 

A^CA and S -> CB (366) 



and the total angular momentum related to invariance under the following 
operation: 

A^AC and B ^ BC and d'' ^ D^^{C^)d'' (367) 



(where C is a constant SU{2) matrix). However they are not much help in 
simplifying the equations of motion or even solving them. In order to go fur- 
ther, and to keep numerical analysis to a minimum, we use the perturbation 



method of Lagrange pO| familiar in celestial mechanics to compute approxi- 
mations to the equations of motion. We will describe how this perturbative 
scheme works. 

We are dealing here with a system described by a Lagrangian of purely kinetic 
nature which we will note: 



: Tl + T2 + Tl 

: To + Ti (368) 
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where Tq = Ti + T2 is the free Lagrangian, Tj is the interaction part and 
Ti ^ To because of some small factor in Tj {1/d in this case). Let denote 
the free coordinates of the system, which we take to completely describe its 
state. In our case qf will represent d\ A or B. The free canonically conjugate 
variables are defined as 

Pi-^Uqf,qf). (369) 



Here are taken to be d\ 7l\A), 7l\B), C\A) and 0{B). In the free system, 
which is described by Tq, these quantities are conserved. Indeed, without the 
interaction term, the Lagrangian describes a pair of free, spherically symmet- 
rical tops moving and spinning at constant velocity. 

Adding Ti to Tq complicates things and removes those conservation laws. Even 
though the system is still described by the same coordinates = qf, the 
canonically conjugate variables are changed; we find 



= PiW,P?)=P? + Api (370) 

where Ap\ has a well defined expansion in 1/d. The aim of the method is to 
find an accurate expansion in 1/d of the equations of motion of the system, 
enabling one to only keep the dominant terms. This is done using Poisson 
brackets. Let us denote by C^{q^,p^) quantitites which are conserved in the 
free system. In our case, they will be rfj, 7^'(A), TZ\B), 0{A) and 0{B) 
but they could also represent more complicated functions of these quantities. 
C^{q,p) will no longer necessarily be conserved and they can be separated in 
a free and an interacting part as follows: 



C^(gi,pO=C^'(gi,P° + Api) 

^C^{qf,P^) + -^C\qf,pf)Ap,. (371) 

The interacting part can be written in a l/d expansion. The time derivative of 
is given by the Poisson bracket of with the Hamiltonian of the system: 

-^C^ = {C^ H} (372) 

CI X 



where H = Hq + Hi = J2i qiPi ~ {Tq + Ti). Hi being suppressed by a factor 
of 1/d relative to Hq, writing H as Hq + Hi also provides an expansion in 
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this small parameter. Finally the Poisson bracket itself can be written as an 
expansion in 1/d. Indeed, since: 



A = A and A = yfM)A = A_yf^)A (373) 

dqi dqf dpi j ^ dpi ^ dp'j dp^ j ^ dpi ^ dp^ 



d 



'dp\ d d 



then 



{A{q,p),B{q,p)} = Y, 



dA{q, p) dB{q, p) dA{q, p) dB{q, p) 



dqi dpi 



dpi dqi 



E 



OA OB OAdB ^fdAdB dAdB-^^dAp 



n 



0/ V J 



dq^ dp^ dp^ dqf ^ ^ dqf dpf dpf dqf ' ^ dpi 
^^ OAdB dAdB ^^dApy 



{qO,pO+Ap) 



^ dqf dp^ dpf dqf ' ^ dpi 



= {A5}o + {A,5}i 



(374) 



Replacing these expressions in ( |372|) , we find the following expansion in l/c? 
of the time derivative of C^: 



dt 



{^C^'(g°,/)Api,//o(g°,/)}o + 

{(:7V,p°),^^o(g°,P°)APi}o + --- 



(375) 



since {C^^H\\i is automatically of higher order in l/d and {C^^Hq\q = 
exactly. The only asumption made here is that d is large. 

For the case = d^, since Hi does not depend on d^, the last three terms in 
( |375|) vanish, thus 



dt 



d'' = {d\Hi}o 



2A 



Md2 



6'^d^ + + 6'^d^ - ?,Sd}(t X 

i'"'(^^'^C{A)C^{B)D^^{A'^B) 

using the free Poisson brackets of the free system: 



(376) 



{d\W} = 5'- 



(377) 
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{n^{A),n\A)} = -^e^^^n^{A) 
{7^"(A),i:^(A)} = o 

{n^{A),DUA)} = -^e^'"'D,M) 
{C^(A),DUA)} = ^e^''D,aA) 



(378) 

(379) 
(380) 
(381) 

(382) 



where 11^ are the conjugate momenta to d\ (Because of the symmetric nature 
of the free Hamiltonian, the same brackets are true if we replace A hj B 
everywhere. Furthermore all the mixed brackets between A and B are zero.) 

In the case where is n^{A), C^iA), n^{B) or C^{B), we get 



A 



dt 



dt 



C\B) 



2Md 

A 

2Md 

A 



2Md 



^i^c^ihdc^(^A)C'^{B){6» - d'S)t^''^Di^{A)D^^{B) 

+ ■■■ (383) 
^iacgjbd^c(^)^d(^^(^ij _ d'd?)e^''^D^{A'^)D^[B) 

+ ■■■ (384) 



A 



2Md 



e'''''C'{A)D^^{A^B) + e^' D^{A^ B)C\A) 



e'''''C\B)D^^{A^B) + e'^'"' D,,{A^ B)C\B) 



+ ■■■ (385) 



+ ■■■ (386) 



where {C^, Hi}q is exhibited and the dots represent the remaining very com- 
plicated terms which are non zero and actually are not negligible, being of the 
same order in l/d as {C^, Hi\o- 

Our approximation is reliable as long as the separation d between the particles 
is large enough for the conjugate momenta to stay close to their free values. 
As we have already worked with the undeformed product ansatz approxima- 
tion, and neglected the potential, which are both valid for large d, we feel 
confident that we have not lost any meaningful information by making this 
further approximation. If d is kept large we should then find geodesies similar 
(qualitatively at least) to those given by the exact equations of motion. 



The system of equations ( P76|) , ( P83|) , ( P84| ) , (|385| ) and (|386|) is still quite com- 
plicated and we will treat it with one final approximation method, namely the 
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method of variation of constants. This perturbation scheme consists of replac- 
ing the variables in the right hand side of these "conservation equations" by 
their free trajectories. This gives rise to a "variation" of the previously con- 
served "constants" . The procedure can be iterated indefinitely to give higher 
order corrections. One should however always maintain consistency with the 
first (Lagrange) approximation. The range of validity of this further approxi- 
mation is rather hard to define, but it is clear that only slowly varying trajec- 
tories in phase space with large d can be considered. The method of "variation 
of constants" is only useful with respect to the equation for d, where it gives 
the scattering trajectory. The change in the spin or the isospin governed by 
equations (|383|) , (|384|) , ( p85| ) and ( p86| ) cannot be treated with this approxi- 
mation method because of the long range nature of the interaction. The results 
give an infinite change in these angular momenta, which is not reliable. We 
will get back to this problem in the last subsection of this article. 



4-4-2 Quantization of the Skyrmion spin/isospin states 

We now use the semi-classical Bohr-Sommerfeld quantization rules applied to 
the classical motion of the free Skyrmion to construct quantum states corre- 
sponding to nucleons. We take such an unusual route because we need quantum 
states which are described by specific classical trajectories and consquently 
suitable for the approximation of variation of constants. In contrast, the or- 
dinary quantization method does not fit into our scheme since it provides 
nucleonic states as quantum wave functions. It is nevertheless interesting to 
consider it in the light of our previous developments and we will briefly de- 
scribe it in the following paragraph. 

Adkins et al quantized the Skyrmion and constructed spin and isospin 1/2 
states. They considered the Lagrangian for a single spinning Skyrmion (|2ll) 
which can be written as 



L = -M + A tr[iU] 

3 

= -M + 2Aj2 a- (387) 

i=0 

where ai parametrizes A = + ia ■ r with ag + = 1. The time depen- 
dent SU(2) matrix A defines the rotational characteristics of the Skyrmion. 
By doing so they make the low energy hypothesis, describing the system ap- 
proximately using only its zero modes (all radial oscillations, deformations, 
etc. are neglected). This is exactly like the BPS case: geodesic motion on the 
minimum energy manifold. One must solve |3^ the Schrodinger equation on 
the manifold of low energy dynamics (in this case, static solutions), finding 
the wave functions for each state of the system in a highly non-Cartesian {i.e. 
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curved) moduli space. The one Skyrmion system is a very simple and elegant 
illustration of this method. 



By considering the truncation of the system ( P8?| ), the configuration of the 



system is just represented by a point moving on the 3-sphere of the SU(2) 
group defined by the familiar constraint: 

AA^ = l = al + a^. (388) 



This is a straightforward generalization of the system studied in the first exam- 
ple of subsection 4.2.1 where the system was a particle following great circles 
on a 2-sphere. By computing vfa (the variables canonically conjugate to the a) 

r)T 



substituting in the Lagrangian (|387| ), and performing the usual Legendre 



transformation, we obtain the familiar expression for the Hamiltonian of a 
spinning rigid body of mass M and moment of inertia A: 

H = M + 1^7fl (390) 



The usual canonical quantization procedure gives: 
1 ^ 52 

H = M - (391) 
2A ^ daf ^ ^ 



which is the Laplacian over the 3-sphere of SU(2): the Schrodinger equation 
we have to solve on the low-energy manifold comes naturally. Solutions to this 
equation are well known: symmetrical, traceless polynomials in Oi. Statistical 
considerations impose the order of the polynomial to be odd for the Skyrmion 
to be a fermion, and even for a boson. Adkins et al |^ of course chose the 
former. They found the following states four spin/isospin states: 



|p T) = -(ai + ia2) (392) 
vr 

\p i) = --(ao-ia3) (393) 

TT 

|n T) = -(ao + ias) (394) 

TT 

\n [) = (ai - 102) (395) 

TT 
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These wave functions are only suitable for a completely quantum mechanical 
treatment of the Skyrmion-Skyrmion system. As stated earlier we follow the 
semi-classical route and will obtain these same states via the Bohr-Sommerfeld 



method Ml . 



We also start from the classical Lagrangian of equation ( 1387] ). Using the ex- 
pression of the matrix A as a function of the usual Euler angles: 

^ = g-i"r3/2 g-i/3r2/2 ^-l^r;l2 ^ggg^ 

where a, 7 G [0, 1tx\ and /3 G [0, tt], or for the 



ao = cos( 2 ) cos( 2 ) 


(397) 


oi = sm( 2 ) sm( 2 ) 


(398) 


02= sm(2)cos( 2 ) 


(399) 


03= cos(2)sm( 2 ), 


(400) 



we find in another form the Lagrangian for a rotating rigid body: 

L = -M+^A[d2 + + 7^ + 207 cos (401) 

The angles a, /3 and 7 obey the following equations of motion: 

« + 7 cos /? — 7/3 sin /? = 

/3 + d7 sin /5 = (402) 
7 + a cos [3 — a(3 sin /? = 

and refer, in connection with the rigid body system, to the motion of the body 
fixed axes relative to the laboratory fixed axes. In our case however, they have 
a different interpretation. Indeed, contrary to an ordinary rigid body which 
has only one conserved vector quantity in the laboratory reference frame, 
(namely the total angular momentum) the quantized Skyrmion will have two 
such conserved quantities: the spin and the isospin. The Euler angles will 
then specify those two quantities. In fact we can say that the angle /? roughly 
will fix the relative orientation of these two vector quantities, while the time 
derivatives of a and 7 will fix their magnitudes. For simplicity without loss of 
generality we will choose the axis of polarization of spin and isospin parallel 
to the axis 3 of isospace. Following the convention used by Adkins et al 
we will chose the following expressions for the spin and isospin generators: 
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/3 = -2A.nl{A) = -2A7^^ (403) 
^3 = 2ACl{A) = 2AC\ (404) 

Following our choice of quantization axes, we take (3 to be either or tt and 
look for solutions to the system of equations ([402|) which are now given by: 



ci ± 7 = (405) 
7 ± a = (406) 

The variables a and 7 actually decouple now and the spin and isospin gener- 
ators are given by (if P is constant) 



J3 = A(d + cos Pj) 

= A(7 ± a) (407) 
5*3 = —A{a cos /5 + 7) 

= -A(7±d), (408) 

depending on whether /3 = ("+" sign) or n ("— " sign). It is then easy to 
choose (3 so that spin and isospin are parallel or antiparallel, and a and 7 so 
they are positive or negative. There only remains to apply the quantization 
rules of Bohr-Sommerfeld to the system. In the original problem of the hydro- 
gen atom, this method was used to compute the allowed radii of the electron 
orbits. In the case of the Skyrmion, the radius of the orbits is already defined 
since the system moves at constant angular velocity along trajectories which 
are the great circles of SU(2), therefore having radius 1. The quantization 
condition will fix the angular velocity so that the spin and isospin have the 
right value, namely 1/2 for nucleons, 3/2 for the first nucleon resonance and 
so on. The Bohr-Sommerfeld quantization condition, derived by path integral 
methods, is ||103|] 



W = J2Ji = (n + Oh, n = 0,±l,±2,--- (409) 



where the Ji are the action angle variables and is a correction factor arising 
from the functional integral over Gaussian fluctuations about the classical 
trajectory, which we will neglect. The action-angle variables Ji are defined by 

J, = jp,dq, (410) 



where pi is the momentum conjugate to the coordinate q\, and the integral 
is taken along a closed path followed by the system during one period in the 
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plane {qi,Pi) of phase space. In our case qi represents the angles a, j3 and 7. 
Due to the cyclic nature of the angles a and 7, and are readily computed: 

Jq, = da = 27rpQ, = 27rA[d + 7 cos/5] (411) 
•A = 'j'Pi^i = 27rp7 = 27rA[7 + dcos/5]. (412) 



and by fixing (3 to either or tt, Jp is 0. Then we find, applying the quantization 
rules, 

Ja + J^ = nh, n = 0,±l,±2,---. (413) 



The spin and isospin generators can also be expressed as functions of the 
action angle variables using (|407D , (pSD , (M ) and (|4T2|) : 



4 = ^ (414) 
ft=-i. (415) 



By symmetry, the Skyrmion has equal magnitude of spin and isospin, |/; 
I J3I, so for n = 1 in ([413|) we have \Ja\ = | J7I = 1/2. 



3| 



In the case where (3 = 0, then using (|407D and 



/3 = A(d + 7) = 2Aw = ±1/2 (416) 
^3 = -A(d + 7) = -2ACU = Tl/2 (417) 

with (a + 7)/2 = 0(t) = cut + 0o, a solution of ( TO ) and ( pej) . This type 
of angular motion produces Skyrmions with spin and isospin antiparallel, and 
the states \p |> and \n t>. The proton state corresponds to isospin +1/2 
along the 3 axis in isospace while the neutron corresponds to isospin —1/2. 



Replacing this solution in equation (|396|) gives 



A = cos 0(t) -isin0(t)r3 (418) 

where a; > corresponds to the state \p l> and u < to \n 
In the case where P = n, 
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/3 = A(d-7) = 2Acu = ±1/2 (420) 
^3 = _A(7 - d) = 2Kuj = ±1 /2 (421) 

with (a — 7)/2 = '?/'(t) = ut + ipQ and the corresponding matrix 



^ = i[sinV^(t)r^ - cos^(t)r2] (422) 

= _ie-iV'(t)rV2^2gi^{t)r3/2 ^^23) 

represents the state |p t> with u > and J,> with u < 0. We can see the 
similarity of the forms for the matrices A corresponding to nucleon states and 
the wave functions obtained in the directly quantum version of Adkins et al 
( |395|) . Finally there remains to compute the magnitude of uo so as to have spin 



and isospin 1/2. This is done by solving the following equation: 

l/gl = - = 2k\uj\ ^uj = \r^ 50-100 MeV, (424) 



using the value of A from reference POU . The uncertainty is caused by the values 
of /tt and e (via the moment of inertia A) which are subject to variations 
depending on which observables are chosen to be best reproduced by the 
model. 



4^.4^.3 Nucleon-nucleon scattering 

We are now ready to compute nucleon-nucleon scattering trajectories. As men- 
tioned earlier, we use the approximation method of "variation of constants" 
and simply replace into the right hand side of equations (|376| ), ( P83|) , ( |384|) 



(|385|) and (|386|) the semi-classically quantized trajectories found in the pre- 
vious section for spin and isospin, as well as taking d{t) corresponding to its 
free trajectory, a straight line at constant velocity. To calculate the change 
induced in the previously constant quantities we integrate the equations from 
t = — oo to t = +00 over one free trajectory. This will enable us to obtain 
scattering angles for the trajectories. 

As we mentioned earlier, this computation scheme does not work for the time 
evolution of the spins and isospins. Equations ( |383D , ( p84| ), p85| ) and p86| ) 
all have on the right hand side a factor 1/d which behaves like 1/t for large 
values of t since d{t) = v t + 'y. Then 

_ 1 1 (425) 



and similarly for TZ^. Then TZ and C change by an infinite amount between 
t = —oo and t = +oo (because the right hand side of the equations integrate 
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to a log which varies very slowly). This means that our approximation scheme 
is too crude or that our treatment is not valid for zero pion mass. Indeed, if 
7^ then the usual Yukawa cut-off factors e~"^'^'^ will arise on the right 
hand side of the equations of motion, and will render the changes in C and 
TZ finite. There does not seem to be any easy solution to the problem so we 
will not discuss spin/isospin changes further. We will just say that results are 
(roughly) compatible with the exchange of charge carrying (pions) and spin 
carrying (vector mesons) as intermediate particles. 

By contrast, equation ( |376| ) shows that d/ dt behaves like and the 

approximation method works well. We present below the results for the scat- 
tering of nucleons for some particular cases of the initial polarizations, using 
our semi-classical formalism, namely scattering of particles whose spin and 
isospin are polarized along the z axis 

A{t),B(t) = cos0(t) — isin0(t)r'^ for the states |p |> and |n ]'> 

or i[sm-ip{t)T^ — cos-?/'(t)r^] for the states |p |> and |n |> (426) 

and the relative motion initially given by 

d{t)=vt + ^ (427) 



where 7 is the impact parameter vector if v and 7 are chosen orthogonal and 
the time of closest approach is at t = 0. These give simple two-dimensional or- 
bits which can be obtained analytically. The tensorial nature of the interaction 
implies that the forces depend on the angle between the axis of separation and 
the spin polarization. If we choose the spin polarization along an axis tilted 
with respect to the normal to the initial scattering plane we get complicated, 
three dimensional scattering trajectories. 

It is important to observe that in our formalism an additional parameter 
arises which describes the initial state of two incoming, polarized nucleons. 
This parameter, along with the impact parameter, the initial velocity and the 
direction of polarization, actually selects the particular scattering trajectory 
followed by the nucleons. The parameter describes the relative orientation of 
the Skyrmions at a fixed (initial) time. It plays in fact a role similar to a hidden 
variable. An incoming pair of physical nucleons, in our formalism, has a fixed 
value for this parameter, which is only "measured" after the scattering takes 
place. In a physical experiment consisting of incoming beams of nucleons giving 
rise to collisions or scattering of pairs of nucleons, the value of this parameter 
will be uniformly distributed. A similar parameter arises in the case of the 
scattering of BPS monopoles. This parameter enters the computations via S 
and e which are defined by 
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(428) 
(429) 



in self evident notation. 



There is an immediate separation of the scattering into two cases, depending 
on whether Dab(^^-B) is time independent or not. When it depends on time for 
large values of the ratio (uj/v) there is an exponential suppression of the scat- 
tering where v is the relative velocity and 7 the impact parameter. This is quite 
evident: for slowly translating Skyrmions, the prescribed rotations imposed by 
selecting semi-classically quantized nucleon states have the effect of averaging 
the interaction to zero. Interactions which depend on the relative orientation 
of rapidly spinning bodies are common in the classical or semi-classical treat- 
ment of soliton systems. A similar numerical example is presented by Piette 
et al ||104|| in the case of "baby-Skyrmions" (Skyrmions in a two dimensional 
space), where classical trajectories and scattering of spinning solitons is stud- 
ied. They too exhibit an "oscillatory" interaction which nonetheless generates 
a net force between the particles and non-trivial scattering. 



First case: D.^\^{A^B) time independent We first present the expression 

for time variations of the previously constant relative momentum p = {M/2)d 
for the scattering of protons or of neutrons with various spin polarisations. 



P T P T 

— / = --— cos(25)ci'^ 
pip I ^ 



(430) 



{zz) pUI ^ (431) 

= (— sin(5), cos((5), 0) 

4pk = _ rjk ^ 4^k^ . J _ gjk/^ . J) 21 

(m) T ^ i ' (432) 

= (sin((5), - cos(5), 0) 



The right hand sides can be interpreted, at this level of our approximation, as 
coming from a spin-spin channel and a tensor channel interaction. We stress 
that this is only a correspondence: the true effect of the kinetic term is to 
supply a non-trivial connection in the geodesic equations on the low energy 
sub-manifold and not to modify the potential. 
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To find the actual change in p and hence the scattering angle, wc integrate 
these equations from t = — cxd to t = +00. For the cases of scattering of 
protons or neutrons on each other respectively we find that the scattering 
angle depends on the variable S which corresponds to the phase lag between 
the rotation of A{t) and B{t). Straightforward integration gives 



= -Aa;2 cos(25) 



' ■dt + —v^ 



2^2 + ^2)3/2 2 



(433) 



^2(^2^2 + ^2)1/2 ^2(^2^2 + ^2)1/2 



(434) 



This yields 



p^(+oo) = -2Au;2^os25^ + 

j''v 2 



(435) 



from which we calculate the cosine of the scattering angle 



cos^ : 



p{+oo) ■ p{-oo) 
|p(+oo)||p(-oo) 



4A 



^ 16A2 



+ LU^ cos2 2S) 



1/2- 



(436) 
(437) 



Second case: Dah{A^B) time dependent For the cases of collisions be- 
tween protons and neutrons, the expression for p''{+oo) is more complicated, 
because of the time dependence of Dab(^^-6). The expressions for the scat- 
tering of p on n each contain a time dependent A^B. When integrated these 
yield an exponentially suppressed variation in the dimensionless group {cuj/v) 

~ e-^^^l (438) 



Thus in the limit v ^ we get negligible scattering in these cases. These 
exponential suppression factors appear in the solutions via G functions which 
can be expressed from Bessel functions by 

GnA^, y) = -^x'^K^ix, y) (439) 
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where K^{x) is the Bessel function of the second kind. We get 



pint 



-j-p^ = —pr cos 4a;t + 2e (f 
dt 



(440) 



Which gives after integration 



Px(+00) 



Mjv 
2Au}' 



sin2e G'i,i(4a;,7/^;) + ^ 



Py(+oo) = f^cos2e Go,i{^u;n/v). 



(441) 



pint 



f'^ = ( - sin(2a;i + e), cos(2u;t + e), 0) 



(442) 



similarly integrates to 



Px(+00) 

Py( + 00) 



2Au>- 



■sin2e[5Go,2(4a;,7/^) + '^G^,2{^,^/v) 



+f^G2,2i^,^/v) + ^Gs,2i^n/v) 



I Mv 



2Aa;2 



cos2e[^Go,2(4a;,7^) + ^16-1,2(40;, 7/^) 

+ ^G2,2(4a;,7^) + :^G3,2(4a;,7^y 



(443) 



and 



{Hi) p in i 



P = (- sin(26,;i + e), - cos(2wi + e), 0) 



(444) 



to 



Px(+00) 

Py( + 00) 



2Aa;2 



2Aa 



sin2e[f?Go,2(4a;,7A) - f Gi,2(4u;, 7^) 



+ i|G2,2(4o;,7^) - ^G3,2(4a;,7H 



+ 



Mv 



cos2e[-^Go,2(4a;,7/^^) + i|Gi,2(4a;, 7^) 
-^G2,2(4^,7^) + :|G3,2(4a;,7/ 



(445) 



This is to our knowledge the first analytical calculation of nuclcon-nucleon 
scattering from essentially first principles, without recourse to ad hoc models 
or potentials. To calculate the classical scattering cross-section we need to 
compute the scattering for all different polarizations relative to the initial 
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scattering plane. This would comprise a different project which would probably 
be best achieved by numerical methods. Therefore we are unable at this point 
to make a direct comparison with experiment. 

Let us now make a few remarks on our results. In the limit that the initial 
velocity vanishes, for fixed uj and 7, we recover 90° scattering. This is, how- 
ever, not surprising as it is a property also shared by the Coulomb and many 
other interactions treated within our approximation. 90° scattering is hardly 
remarkable except at zero impact parameter, where of course, it is impossible 
to avoid the region of close proximity of the nucleons and it is important that 
the configurations pass through the minimal, toroidal configuration. 

We have made several approximations in our treatment, which deserve some 
discussion. First we want consider the method of variation of constants. To 
check its accuracy, we observe that equation ( [430| ) is a simple Kepler problem 
(for this particular equation). We solve it directly to find the exact value of 
the scattering angle and then compare wiith the result obtained by variation 
of constants. The exact scattering angle 9 is: 

cos e = 1 - ^^.^,.,4 • (446) 



Defining a parameter x 
cos2 25A2cu^ 



and choosing v and 7 so that x is much smaller than 1 we get from equations 
(pTp and (TO: 



cos e = -j=^== ^ 1 - 8x + O(x') (448) 
1 — 4v 

cos e = — — ^ ~ 1 - 8x + O(x'). (449) 
1 + 4x 

With the previously defined values and definitions we get 

X ~ 10"^° MeV-2 X < 1. (450) 

7^w* 



This imposes very loose restrictions to v and 7. For these conditions, the 
approximation seems to work well and respect the other assumptions. 

The second main approximation is the expansion in inverse powers of d, ne- 
glecting all terms beyond the dominant contribution (from the kinetic term). 
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We have from the start restricted ourselves to large separations between the 
Skyrmions. In this regime, the Skyrmions are well described to leading order by 
the product ansatz and the manifold of collective coordinates is parametrized 
by the variables of the product ansatz. The induced metric, calculated to 
leading order behaves as 1/d while the induced potential behaves as In 
principle there can be a region where the contribution of the induced metric 
dominates and we can neglect the potential. We find that the metric induces 
an interaction which can be interpreted, within our approximation method, as 
a spin-spin and a tensor interaction. Unfortunately it seems that the domina- 
tion by the metric term is not physically realized. The induced kinetic term is 
multiplied by essentially the frequencies of angular rotation of the Skyrmions: 



Tint = ^^e^-ejbd ^a^^) ^d^^) ^^ij _ ^i^j) D^b(At5) + 0(l/rf2) 



d 



(451) 



while the potential term has two extra powers of the separation in the denom- 
inator: 



d^ 



cose){3{n-df -1) 



(452) 



So 



-^^int — ^int + Mnt 
-2 f2, ,2 



implies that Tint ^ ^int if 

d > — . (454) 

UJ 



Thus for the kinetic term to dominate, the frequencies should be much larger 
than the separation. This corresponds to a region of validity for a separation 
of about 3 fm and greater. However there is much latitude available since the 
values of /^r and e which go into determining cu are fixed only by choosing two 
experimental inputs, f^^, e can vary as much as 10-30% thus we do not feel 
overly concerned with exceeding the regime of validity. Our approximation 
would of course be better justified for the case of A — A scattering where 
u; ~ 300 MeV corresponding to a separation of 1 fm. In any case, we do not 



123 



believe that it is physically reasonable to consider the scattering of nucleons 
with the metric term alone and we expect a contribution from the potential 
term which is of the same order of magnitude. We do not expect, however, 
any great, qualitative modification of the scattering upon inclusion of the 
potential term, it is of similar strength but actually contributes only in the 
tensor channel for the case of massless pions. 

We have shown how to formulate the nucleon states within the semi-classical 
approximation. We have treated the scattering and computed the scattering 
angles in a systematic perturbative approximation. Future work should in- 
clude consideration of a non-zero pion mass, which leads to a central channel 
interaction, a better control of the perturbative method, a departure from 
the product ansatz, and a proper treatment of the region of close proxim- 
ity, to test the validity of our formalism in the phenomenology of low energy 
nucleon-nucleon scattering and of the static quantum states in this sector. 
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